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FOREWORD 


This document is the final report for Contract NAS5-11106 to 
perform the post-launch and prime data reduction and data analysis 
of the AS&E grazing incidence x-ray telescope experiment on the 
pointed section of the OSO-IV satellite, and the AS&E rotating 
wheel experiment. The contract covers the period 19 November 
1967 to 30 April 1970. The spacecraft was launched on 18 October 
1967 and the experiment has operated successfully throughout its 
lifetime. At the present time it continues to return useful informa- 
tion on the location of sources of solar x-ray activity, although 
there has been some loss of sensitivity. 

The scientific objectives of the pointed x-ray telescope experiment 
are to study the structure of solar x-ray sources in several wave- 
bands over an extended period, and to study the relationship of 
these sources to problems of solar physics. The objectives of 
the rotating wheel experiment are to locate sources of soft x-rays 
in the celestial sphere and to determine their intensities and 
spectral properties within several wavebands. 

The reduction and analysis of the OSO-IV data has taken place 
in several phases. Prior to the launch of the spacecraft, computer 
programs for data reduction were developed and tested. This 
phase covered the period 3 0 June 1966 to 31 March 1967. 

The work reported in this document covers the second and third 
phases. Post-launch data reduction commenced 19 November 1967, 
one month after the launch of the spacecraft, and continued until 
3 0 June 1968. Preliminary data analysis commenced during this 
phase. The character and quality of the scientific results to be 
gained from the experiment were defined during this phase of the 
analysis and preliminary results were reported to the scientific 
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community. The prime data reduction and analysis phase began 
1 July 1968. During this phase, routine data reduction was 
completed for all data received in the form of computer readable 
experimental data and spacecraft attitude magnetic tapes 
(essentially the period from experiment turn on, 25 October 1967, 
to spacecraft tape recorder failure, 12 May 1968). This data 
has been examined, significant events have been selected, 
and the data for this period has been analyzed in detail. The 
results of these analyses have been reported to the scientific 
community as they became available. 

The Principal Investigator on this contract was Dr. Riccardo 
Giacconi. The Project Scientist for the data reduction and 
analysis of the OSO-IV AS&E pointed x-ray experiment, covered 
in Sections 1. 0 through 6. 0 of this report, was Dr. Giuseppe S. 
Vaiana. The Project Scientist for the data reduction of the 
AS&E x-ray stellar rotating wheel experiment on OSO-IV 
(Section 7. 0) was Dr. Herbert Gursky. 
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1.0 INTRODUCTION 


The grazing incidence x-ray telescope on the OSO-IV spacecraft 
obtained solar images in the 2. 5 to 12 X wavelength region, 
nearly continuously from 27 October 1967 to 12 May 1968. The 
spatial resolution of this instrument is on the order of one arc- 
minute. The time resolution can be as short as 5. 12 minutes, 
but is usually 20. 5 minutes. 

The instrument has sufficient spatial resolution to locate and 
identify the sources of the solar x-ray emission and to provide 
an estimate of the source size. The instrument can measure 
the x-ray emission from a particular region even while the 
majority of the solar x-ray flux is coming from other regions. 
Consequently, the x-ray telescope is able to observe time 
variations in the x-ray flux emitted by an individual active 
region and to follow the evolution of this activity for several 
solar rotations. 

The body of this report is devoted to the results of the analysis 
of the OSO-IV data. This analysis was conducted concurrently 
with higher resolution rocket investigations as part of the AS&E 
Solar X-ray Astronomy Program. 

1. 1 AS&E Solar X-ray Astronomy Program 

The AS&E program of x-ray observations of the quiet and active 
solar corona, which began ten years ago, has been directed 
towards the simultaneous acquisition of data with high spatial, 
spectral and temporal resolution by means of grazing incidence 
x-ray telescopes. During this period, the grazing incidence 
optics have evolved from the first crude versions to high 
precision telescopes such as the S-054 instrument to be flown 
on ATM (better than 2 arc-second resolution). Figure 1-1 
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15 MARCH 1965 
( AS&E / GSFC ) 


20 MAY 1966 

( GSFC ) 




8 JUNE 1968 

( AS&E ) 


Ha 

8 JUNE 1968 

( COURTESY OF ESSA ) 


Figure 1-1 Three x-ray photographs of the sun showing the evolution in 
x-ray telescope performance over the past few years. Also 
shown is an Ha picture taken at the time of the 8 June 1968 
rocket flight. The flare is the brightest feature close to the 
center of the disk (overexposed in this particular exposure). 
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illustrates the progress of x-ray imaging experiments during the 
last several years. 

The OSO-IV experiment is part of a systematic study of the 
quiescent and active x-ray corona by rocket and satellite 
observations. The following list summarizes these observations: 

First use of grazing incidence optics (Aerobee 150 
flown in November 1963) 

1 2 

X-ray photography of the sun at solar minimum ’ 

(Aerobee 150 flown in March 1965) 

Long-term monitoring of x-ray emission by active 
regions with arc-minute resolution^’ 5, 6 
(OSO-IV satellite launched in October 1967 and 
still operating) 

First attempt at flare x-ray photography (only partially 
successful because of pointing control failure) 

(Aerobee 150 flown in March 1968) 

Observation of an x-ray flare with few arc-seconds 
resolution. First use of slitless spectrometer 
technique 7, 8 » 9 (Aerobee 150 flown in June 1968) 

Observation of quiescent coronal structures on the 
disk with few arc-seconds resolution and limb 
observations to more than 0. 3 solar radii in height 
(Aerobee 150 flown in April 1969) 

Flare x-ray cinematography to observe the spatial 
development of an x-ray flare and to associate the 
microwave and hard x-ray bursts with centers of 
emission inside the flare region (Aerobee 150 flown 
in November 1969) 

Observation of the x-ray corona near the time of a 
solar eclipse to unambiguously associate large radius 
eclipse observations with disk features. Observation 
of the lunar de-occultation of an active region to 
determine the altitude dependence of the x-ray emission 
with greater precision. (Aerobee 150 flown in March 
1970) 
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With the exception of the AS&E grazing incidence telescope 
experiment on OSO-IV, described in this paper, all of these 
observations were made from sounding rockets. Rocket flights 
have a typical duration of less than five minutes. Thus, 
rocket-borne detectors cannot observe the time variations in 
the behavior of x-ray active regions. In order to observe 
the time behavior of the x-ray activity of individual emitting 
regions, a satellite-borne instrument is required with sufficient 
spatial resolution to distinguish and locate the emitting regions. 

This long-term monitoring with moderate spatial resolution is the 
unique characteristic of the x-ray telescope experiment on 
OSO-IV. The still photographs from recent rocket flights, 
however, show details of x-ray emitting regions with a spatial 
resolution orders of magnitude better than available with 
the OSO-IV experiment. Therefore, in order to relate the 
results of the OSO-IV analysis to the higher resolution photo- 
graphs, it is appropriate to summarize some of our rocket 
results. 

Our three most recent rocket flights have revealed the existence, 

8 9 

within x-ray emitting regions, of structures ' with dimensions 
on the order of a few arc-seconds. The structure of an emitting 
region may closely resemble that of the corresponding Ho- plage 
(Figure 1-1 and 1-2) when seen on the disk. 

Examination of the spatial distribution of the x-ray and Ho emission 
near the limb (Figure 1-3) shows that the x-ray emission extends 
to considerable heights (100, 000 to 150, 000 km) in the corona 
above the active region. This three-dimensional structure 
may take the form of loops connecting portions of the same active 
region or of nearby active regions. The configuration of the x-ray 
emitting plasma appears to be governed by the magnetic field. 


1-4 



Figure 1-2 An x-ray exposure of the sun taken 8 June 1968 at 1745 UT 
while an importance In flare was in progress. The exposure 
was made with the AS&E 34 cm 2 x-ray telescope. 
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Figure 1-3 Comparison between x-ray and Ha emission shows the three 
dimensional structure of active regions. Three loops can be 
seen to connect active regions in the left set. 
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The structure is quite complex and does not correspond to a 
simple bi-polar field configuration within one active region. 

The x-radiation of the flare of 8 June 1968 derives neither from 
a single bright point, nor from a diffuse blob of unconfined 
plasma. Instead there is a general correspondence between 
the regions of x-ray and Ho emission. (A detailed comparison 
of the spatial distribution in Ho- and x-rays is shown in Figure 
1-4.) A significant exception to this is seen where part of the 
brightest x-ray emission bridges a dark lane between two Ho- 
emission centers located on opposite sides of the zero longitu- 
dinal field line. Roughly 50% of the total x-ray emission of 
the flare comes from a region of about 15 arc-seconds dimension 
about this point. Figure 1-5 presents a comparison between the 
photospheric magnetic structure, the x-ray structure and the Hey 
structure. The three rectangles correspond to the same portion 
of the flare region and are at the same scale. We note the co- 
incidence between the bright x-ray emission and the regions 
of higher longitudinal field gradient. 

Exposures taken in the slitless spectrometer configuration have 
provided spectra from 3 to 15° of individual active regions. 

The spectra obtained confirm the marked difference between 
the flare region and other active regions (Figure 1-6). The 
spectrum of the flare region is substantially harder than that of 
the active region, as expected. Significant differences, 
however, are also present between the spectra of non-flaring 
active regions. 

On the rocket flight of 8 April 1969, the structure of the x-ray 
emitting regions of the general corona was observed on the 
disk for the first time (Figure 1-7). In addition to the bright 
limb, familiar from earlier x-ray photographs, the disk is seen 
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8 June 1968 , 1743 UT 




Figure 1-4 Comparison of the flare region in Ha and x-rays. 
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Figure 1^5 Comparison of the flare in magnetic field, x-rays, and in Ha. 
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6 Spectra obtained by the slitless spectrometer in the 8 June 1968 
rocket flight. The flare spectrum is substantially harder than 
either of the other two active regions shown. The spectrum of 
region 9445 is significantly harder than of region 9443. 

The ordinate is left blank pending absolute energy calibrations. 
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Figure 1-7 An x-ray exposure of the sun taken 8 April 1969 by AS&E's 
34 cm 2 , rocket- borne, x-ray telescope. In addition to the 
active regions and bright limb, both point-like and diffuse 
structures appear on the disk. 
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to display inhomogeneous structure. Small, concentrated, point- 
like features and more diffuse structures are both noticeable. 

A preliminary examination indicates that the latter structures 
are more diffuse than the chromospheric structure evident in 
CaK spectroheliograms. The point-like structures correspond 
to particularly bright features of the chromospheric network 
with bi-polar fields. The nature of these inhomogeneities in 
the low corona is not yet clear. The relationship between these 
structures and the magnetic sector structure and longitudinal 
velocity structure of the solar wind requires further investigation. 

1.2 Scientific Objectives 

The scientific objectives of the post-launch and prime data 
analysis phases of the OSO-IV pointed experiment were to 
analyze and interpret the data obtained from the OSO-IV x-ray 
telescope experiment. In particular, (1) to use the spatial 
resolution of the instrument to distinguish individual active 
regions in order to study their location and size during both 
flaring and quiescent periods; (2) to use the extended duration 
of the experiment to study the temporal behavior of the x-ray 
flux from active regions; and (3) to investigate the relationship 
between the soft x-ray emission of active regions and their 
manifestations at other wavelengths (e.g. , Ho and microwave 
radio emission). The specific tasks to be performed encompass 
four areas; routine data reduction, x-ray plage analysis, 
x-ray flare analysis, and scientific collaboration with other 
experimenters. 

1.2.1 Routine Data Reduction 

Upon receipt of GSFC supplied experiment data on magnetic tape, 
AS&E reduced, computed and printed out the data in raster format. 
Duplicate tapes were submitted to NSSDC as required. 
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1.2.2 X-ray Plage Analysis 

Computations were performed to arrive at spectral hardness 
indices, temperature, and electron densities over selected 
active plages. Correlation analysis was conducted of plage 
x-ray brightness with temperature, magnetic field intensities 
and gradients, Ho and other related solar parameters. The plage 
evolution was studied, and the findings were interpreted using 
specific plage models. 

1.2.3 X-ray Flare Analysis 

Graphs were constructed of x-ray flare brightness and x-ray 
flare hardness index as functions of time. Computations were 
made of peak irradiance of flares summed over the entire 
plage subfield, of x-ray flare equivalent temperature and of 
total energy output of flares. Correlation analysis was performed 
of x-ray flare brightness and Ho peak brightness, flare sizes 
and Ho sizes, flare brightness, and other flare parameters. 
Results were interpreted using specific flare models, and an 
attempt to improve flare forecasting was made. 

1.2.4 Scientific Collaboration 

AS&E has established collaborative relationships with other 
OSO-IV experimenters and experimenters of other spacecrafts. 
During the period of performance of this contract, AS&E met 
and exchanged information with these individual scientists 
and with interested government agencies. The contractor made 
available the data and findings, at the discretion of the NASA 
technical monitor, by publishing it in technical journals and 
presenting it at technical meetings. 
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2.0 DESCRIPTION OF EXPERIMENT 


2. 1 Instrumentation 

The grazing incidence x-ray telescope is a two-mirror, image 
forming system of optics. Incident x-rays reflect off the first 
mirror, which is a paraboloid of revolution, onto the second 
mirror, which is a confocal hyperboloid of revolution. The 
radiation reflected off the second mirror subsequently forms 
a true image of the source in the focal plane of the system. 

For x-radiation of wavelengths greater than 2 $, angles of 
incidence are less than the critical angle necessary for total 
external reflection; reflection efficiency is therefore high. 

A detailed description of grazing incidence x-ray telescope 
design and application has been given by Giacconi et al. 111 

The telescope on OSO-IV (shown mounted into the instrument 

in Figure 2-11 is fabricated of electroformed nickel; it is 

approximately 7. 6 cm in diameter by 16 cm long; its collecting 
2 

area is 2. 01 cm ; and its focal length is 83. 6 cm. Its inherent 
resolution is 20 arc-seconds. The grazing angle of reflection 
for rays incident parallel to the axis is 40 arc-minutes for both 
reflections. 

The reflection efficiency of the x-ray telescope was measured 
at the 8. 3 X aluminum line. The reflection efficiency on axis 
was found to be 28%. The theoretical reflection efficiency 1 * 
at 8. 3 8 for two reflections from nickel surfaces pitched at 
equivalent angles to the x-ray beam (that is, 40 arc-minutes 
for both reflections) is approximately 65%. This difference 
between the theoretical and actual efficiencies is considered 
typical of x-ray telescopes fabricated by the electroforming 
technique. The theoretical and measured reflection efficiencies 
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are shown in Figure 2-2. 


Transmission as a function of wavelength for all filters was 

calculated from the values of the mass absorption coefficients 

12 

tabulated by Henke et al. and measured areal densities 

-2 

(gm-cm ) of the filter materials used. Their transmissions to 

8.3 S Al -Kg- x-radiation were then experimentally measured and 

were found to be consistent with calculations. Besides the 

wheel- mounted filters, there is a fixed filter in the optical 

path of 0. 00015 inch Mylar with 2000 i? evaporated Al, adding 

2 

primarily 0. 53 mg/cm of C. 

The detector is a "photoemission-scintillation" detector. 

13 

Lincke and Wilkerson have described its use in the extreme 

ultraviolet. The detector (Figure 2-3) consists of an external 

photocathode and a scintillation detector. The photocathode 

is a conical shell of nickel substrate, the interior of which is 

coated with Csl. This substance has been found, by Lukirskii 
14 

et al. to have an extremely high photoelectric quantum 
efficiency for x-radiation. The photocathode is held at a 
potential of -10 kV. The scintillation detector uses a plastic 
scintillator mounted on a lucite light pipe, and viewed by a 
14-stage photomultiplier. The scintillator is coated with a 
film of vacuum-deposited aluminum; the aluminum film is at 
ground potential. 

Incident x-radiation strikes the interior of the conical cathode 
and causes, with high probability, the emission of one or more 
photoelectrons. The conical shape of the photocathode sub- 
strate structure and the planar shape of the aluminum filter over 
the scintillator form an electrostatic lens which directs the 
photoelectrons into the scintillator. This produces enough 
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Figure 2-3. Photocathode detector chamber. 
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scintillation light to be detected by the photomultiplier with 
reasonable efficiency. 

2.2 Experimental Procedure 

2.2.1 X-ray Telescope 

A schematic representation of the AS&E grazing incidence x-ray 
telescope experiment is shown in Figure 2-4. The instrument 
uses a grazing incidence system of mirrors (an x-ray telescope) 
which form a true x-ray image of the sun on the Aperture Wheel 
in the focal plane of the telescope. The wheel contains four 
apertures of two different sizes (one arc-minute and four arc- 
minutes) which determine the spatial resolution of the optical 
system. The Aperture Wheel was programmed to step auto- 
matically between complete scans of the sun, or to stay at any 
desired position. As the OSO-IV pointed section scans the sun, 
its image moves over the Aperture Wheel. Only x-radiation 
from that region of the sun at which the instrument is pointed 
passes through the aperture for further analysis. 

2.2.2 X-ray Detector Package 

The x-radiation from a single picture element subsequently passes 
through one of three filters in the Filter Wheel located im- 
mediately behind the Aperture Wheel. The filters, together 
with a fixed filter in the optical path of the x-ray telescope 
and the reflection properties of the telescope itself, permit 
analysis of x-radiation in three wavelength bands. These are 
approximately: (a) 2.5 to 12 R, (b) 2.5 to 11 R, and (c) 2. 5 to 
9 $. The Filter Wheel can be programmed to step automatically 
between complete scans of the sun. It cycles through four 
positions:, three filter positions and a calibration position. 

It can also be programmed to stay at any desired fixed position. 
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Beyond the filters', transmitted x-ray quanta strike the photo- 
cathode where photoelectrons are emitted and are then 
accelerated into the scintillator viewed by a photomultiplier 
tube. The counts accumulated between data sample times 
are thus proportional to the number of x-rays emitted from the 
corresponding picture element of the sun. A total of 1920 such 
sets of accumulated counts comprise the data from one complete 
scan of the sun. 

2.2.3 Data Transmission 

The pulses from the photomultiplier are subsequently counted 
and the result telemetered to the ground. These sets may be 
reconstructed on the ground into an array of 48 by 40 numbers 
which constitute an x-ray map of the sun. The position of a 
given number in the array corresponds to a unique position on 
the sun and the magnitude of the number corresponds to its 
brightness in x-radiation. Data accumulation time between 
these numbers is 0. 14 seconds. Due to readout delay time, 
the actual time interval required to accumulate the data for 
one x-ray filtergram is 5. 12 minutes. Thus, if the Filter Wheel 
is cycling (as is usually the case), 20. 5 minutes elapses 
between consecutive filtergrams through the same filter. 

Table I lists the instrumental characteristics of the OSO-IV 
telescope experiment. 
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TABLE I 

OSO-IV INSTRUMENTAL CHARACTERISTICS 


SPECTRAL RANGE: 

2.5 - 12 £ 

(2. 65 

2 

At gm/cm Beryllium filter) 


2.5 - 11 X - 

(1.4 

2 

(x gm/cm Mylar filter) 


2.5-9 8 

(9.4 

2 

fj. gm/cm Beryllium filter) 


SPATIAL RESOLUTION: 1 arc-minute 

4 arc -minute 

TIME RESOLUTION: 5.12 minutes 

20. 48 minutes 

4 2 

1.4x10 photons/cm -sec 

3.5 x 10 ^ ergs/cm^-sec (10^°K) 


MINIMUM SENSITIVITY: 



3.0 CALIBRATION PROCEDURES 


3. 1 Ground Processing of Telemetry Information 
The initial output "raw raster" consists of solar x-rays, 
background counts that are due to spacecraft generated noise 
(a function of raster position), and energetic charged particle 
counts. The latter counts depend upon spacecraft orbital 
position and time, but are independent of raster position. 

To determine the number of true counts due to solar x-rays, 

N, as a function of raster position, one must subtract from 
the raw raster an "average background raster, " which contains 
the spacecraft generated noise and the average charged particle 
counting rate.' 

The average background raster is generated from raster scans 
55 

utilizing the Fe calibration source, instead of the x-ray 
filters, in the optical path of the telescope. Rasters obtained 
when the spacecraft was passing through the South Atlantic 
Anomaly are removed entirely from data processing. 

The number of true counts, N, accumulated during the 0. 14 
second live time interval between data words is related to the 
spectral distribution function, <£, at the spacecraft by the ex- 
pression: 

N (a, b, t, T) = 0. 14 aJf (b, X) <t> (a, t, T, X ) d X (3. 1) 

where a = Aperture Wheel index, indicating the angular 
resolution mode of the experiment; 

b = Filter Wheel index, indicating the transmission filter 
in position, and thus specifying an effective wave- 
length passband; 

t = time; 

T = spectral parameter indicating the nature of the source 
spectrum; 
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A = sensitive collection area of the telescope aperture; 

X = independent variable, wavelength, in £; 

<t> = spectral distribution function of radiation at the 
spacecraft (units are photons/cm 2 -sec-S); and 

F = instrument detection efficiency function. 

The function of F (b, X ) is expressed by 

F (b, X ) = r ( X ) e (t) q ( X ) x (b, X ) . (3.2) 

where r = telescope reflectance; 

£ = scalar part of the photocathode detection efficiency, 
calculated from the in-flight calibration mode; 

q = photocathode quantum efficiency indicating the 
probability that at least one photoelectron will be 
emitted upon irradiation by a single photon of 
wavelength X; and 

x = filter transmission function. 

The functions F (b, x ) for the three filters used in this experi- 
ment appear in Figure 3-1. A more detailed discussion of these 

15 

relations is available elsewhere. 


3. 2 Determination of the Calibration Counting Rate 

The average background raster is constructed in the following 

manner. The constant, N is the average counting rate due 
55 cal, 

to the Fe calibration source, and is proportional to the x-ray 
detector efficiency. The quantity Njq sun is generated by 
averaging the counting rate of 40 raw raster positions, where the 
effect of solar x-rays and spacecraft generated noise are 
minimized. In practice these positions are usually the top 
2 0 and bottom 2 0 lines in the first half of a raster. 80 lines 
per raster were used for data from later in the spacecraft's 
lifetime, to improve statistical precision due to falloff of the 
experiment detector efficiency. The least sensitive mode of 
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Figure 3-1. Instrument detection efficiency function for 
the three filters used in the experiment. 


3-3 


is generated by 


the experiment is used. The quantity N^q 

averaging the counting rate in the corresponding 40 positions 

of the calibration rasters taken during the same time period. 

N , is then calculated by subtraction of N. rt from the 
cal _ 40 sun 

corresponding ^ 

N , can now be subtracted point by point over the entire 
cal 

"average calibration raster, " obtained by summation of all 
calibration rasters during the specified time interval. The 
average background raster for this interval is thus produced. 

An x-ray filtergram of the sun, containing only solar x-rays, is 
then obtained by point by point subtraction of the appropriate 
average background raster from the raw raster. 

Time variation in the sensitivity of the detector is monitored 

in the calibration mode. Figure 3-2 shows the time variation 

of the detector sensitivity for the first six months of the 

satellite’s lifetime. The detector sensitivity. N . has-been 

55 cal 

corrected for the natural decay of the Fe source, whose mean 
life is 3.90 years, and normalized to its value on 2 7 October 
1967, the starting date of routine solar scanning. It is 
apparent from the graph that a one-third reduction of the 
value of N ^ occurred during the first few weeks of operation. 
The exact cause of this decay is unknown, but probably is 
associated with deterioration and/or coating of the Csl photo- 
cathode surface. The sensitivity then decreased more slowly 
thereafter. The sensitivity remained high enough to allow 
monitoring of the x-ray flux of active regions throughout the 
period of this report. 

3.3 Further Corrections to the True Counting Rate 
3.3.1 Methods of Approximation 
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Figure 3-2. Time history of the relative x-ray sensitivity of the 
instrument, N ,, normalized to 27 October 1967. 




The true counting rate per data word on the reduced x-ray 

filtergram must now be corrected for the time variation in 

detector sensitivity shown in Figure 3-2. From such a graph 

a daily multiplicative factor is computed that is simply the 

normalized rate of N . on any day to the value on 27 October 

cal 

1967. It should be emphasized that all such corrections are 
relative to the initial in-flight conditions of the spacecraft. 

To derive information on the total flux of an x-ray emitting 
region on a filtergram, ring integrals are formed about the center 
of the region as described in Section 4.4. For ease of compu- 
tation and to avoid overlapping active regions, it is convenient 
to truncate the integrals at arbitrary radii from the center of 
the emitting region. Thus it is necessary to determine correction 
factors for the one and 4 arc-minute aperture stops to estimate 
the remaining flux in the uncounted rings. 

Anticipating a result from Section 5.2.3, we find by comparing 
the observed x-ray intensity distributions of flares with the 
telescope response function that the flare brightness dis- 
tributions are on the order of one arc -minute or less. Because 
our smallest telescope aperture diameter is one arc-minute, we 
make the assumption that the brightness distribution of flares 
represents the point source scattering function of our instrument. 
We chose the brightest x-ray flare, uncontaminated by adjacent 
regions, that occurred during the period covered in this report 
for determination of this "point source correction. " This flare 
occurred at about 2200 UT on 16 November 1967. We calculated 
ring integrals for rasters near the flare's peak intensity in our 
two aperture modes. A table was produced listing the per- 
centage of the total flux contributed by each ring. Thus to make 
the correction, the flux calculated out to a specified ring is 
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divided by that ring's contribution factor to derive the total 
estimated flux. 

This method is also applicable to any emitting region observed 
on rasters obtained with the 4 arc-minute resolution mode, since 
the majority of the flux from such regions is emitted from areas 
smaller than 4 arc-minutes. 

3.3.2 Absolute Flux Calibration 

In order to convert our measured counting rates to an absolute 
flux scale, we have compared our measurements for selected 
x-ray events with those of other experiments. One of the most 
useful comparisons was with data obtained by the University of 
Iowa detector on board Explorer 33, which operates in the same 
approximate waveband as our detector. We have used these 
calibrations to calculate representative flux values for x-ray 
emitting regions on time scales of interest. Some of these 
data appear elsewhere in this report. 
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4.0 COMPUTER REDUCTION AND ANALYSIS ROUTINES 


4. 1 General 

The data analysis and reduction routines for the OSO-IV pointed 
x-ray telescope experiment are utilized on AS&E's in-house 
IBM system 360/40 Data Processing Facility. The program 
formats are mostly written in PL/I language and may be read 
from PL/I, assembly language, or COBOL. A few short, older 
programs are written in FORTRAN IV. Listings of the important 
data reduction programs are presented in the Appendix. 

Data reduction commences at AS&E when compatible magnetic 
data and spacecraft attitude tapes are received from Goddard 
Space Flight Center. When the tape information has been 
processed onto in-house tapes, the originals are returned to 
GSFC. Routine data processing has essentially been completed 
from the start time of telescope routine scanning, 25 October 
1967, to the time of the failure of OSO-IV’ s second tape recorder, 
12 May 1968. Reduction has been delayed for the periods 
28 February to 6 March 1968 (telemetry tape) and 23-28 April 
1968 and 7-13 May 1968 (attitude tapes) because corrected 
magnetic tapes have not been received from GSFC. 

4.2 Routine Conversion of GSFC Magnetic Tapes to Raster 
Format 

4.2.1 Raw Raster Generation Program 

The GSFC telemetry (“Main Frame") tapes are converted to tapes 
called Telemetry for the same time period at AS&E. GSFC 
attitude ("Subcom") tapes are converted to tapes called Aspect- 
Ephemeris (A-E) also for the same time period. The spacecraft 
attitude and telemetry data are then combined onto a single 
"Raw Raster" tape, and the information printed out with the 
use of the Raw Raster generation program. Two print sheets are 
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required per raster scan, one for essentially the attitude informa- 
tion and the other for the telemetry readout, or raster scan, 
itself. 

The following information is included on the attitude sheet: 
the date, the Universal Time (UT) of the first raster data word 
and day of the year to 3 decimals (1 January 1967 = 1), the 
filter wheel position as determined from the Analog Subcommutator 
(ASC) #5, the aperture wheel position from ASC #7, the filter 
wavelength interval, the ephemeris and aspect data at the start 
of the raster scan (angles in degrees and altitude in kilometers), 
the sidereal start time (UT) in degrees, the number of unin- 
telligible raster scans skipped since the last raster on tape, the 
number of data words recorded on the present raster (1920 for a 
nominal scan), the number of questionable and missing data 
words on the raster as set in the telemetry tape by GSFC, and the 
average readouts from the four ASC voltage monitors for five 
equally- spaced intervals across the raster. 

The raster scan consists of two-digit "words" in a hexidecimal 
code. The scan begins in the lower left comer with a 3-4 word 
signature. The recording time increases vertically up the first 
column, down the second, up the third, and so on until the end 
of the raster is reached in the lower right comer. Details on 
the spacecraft telemetry have been recorded elsewhere. ^ 

Figure 4-1 shows a typical raw raster scan. 

4.2.2 Background Raster Generation Program 

The Background Raster generation program uses only the Raw Raster 
tape described in the previous Section. The method for calculation 
of the average background raster is discussed in Section 3. 0. 
Average rasters are computed daily from 25 October 1967 to 
26 January 1968, and weekly thereafter. The time interval was 
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RAW RASTER SCAN STARTING AT 
01 1/2 MIL BERYLLIUM 
11 40 MIL APERTURE 

UT NOV 
LAMBDA 

lit 
1 =t 

1967 

2-00 

UT TIME 
LAMBDA 

2 

2= 

2 11.48 
13.00 

DAY OF YEAR 315 
KN = 1296. 00/ EP SI LON 

1 5 10 

15 


20 

25 


30 

35 40 


48 C303030301C40 1010 1010 1C 10002030 1040203040400000502020102020302030204040202020306 *8 

OOC3COOOOCCOO 101030402050201030400040201 0202070304040302060004020305010206060805 
02020CCC040 1020 3020202C2030 302030 1000203020 1020304030202040406020404030 504030305 
45 03060104060402030101020300040200030603020503020500040302060204000305050504060304 45 

0106030 1020304030200020101000402060400020305030402020401040305020302030402050604 
06030103030504050205030404060302030302040303030102000400060602060403060204010603 
05C 5050 103000202030400 0204020204060 1020400000 50 302000 102040203030503020 103070406 
02 03 030305030403070202C20204030 1030005020401 030302040303020202020 A04040703020606 
40 06020604050508060 30 105C003060405 0104020202000002040 103020 305050703030300050706 04 40 

04020405040406050104020406010301040302040402060401040504040306060605050302090808 
060)60905030508040503040404050301050603070105040303040103040406030403030508040705 
0807060E07070 70303050706050 10004030203040304000302040202020 104050407030602050801 
05C5080BOAC706050S0403060203C0030 50203060002040504020605040004030306080708060807 
35 07CB000 A 14 070C090702 04 030 20 50202 030 104 03 0603 05 07 04 05060 70 00304030 50 3080 106 05 060 A 35 

080B 101117130COF030507C605020302010105030304040403010303040202060006050905070C04 
13132424101115080407050306050503040403040202040303020203040208040506040708050608 
00 182E2831120DQFOCOB050404050004030203050408030105020200030209030503050A060B0707 
1E2448402F2F140P 0908 09 05020404030 10303030303020 1000406 0005080 5050B06020408070905 
30 131C3E52402B160E080C0305040205040304080401050303010603030409d30609020508040A0609 30 

142 546645 8 301C130C080A040509C40202040601050 1030407020 105040204050806050904070805 
192847604 134201206090 1020403020304040300090403030 1090 304050606050407060E05060804 
1 A3 13 A604830130E0706060A0906030504060505030204050 1040804020703040407070906070807 
16333F43 3 628 1813 OB 0608 C40 10504030 3030608030 106050 00307 05 040508030 50908060407 1007 
25 14IA28391F1E12091C05060405050306030203020A03060604090403060308040604050EOC060607 25 

IS 14202817 16051 106030309070804070 10702040402010302070501 0905070A0807050AOC08070C 
OE 15 1613 13 1008061 1060906050 5030508010103020105070104030EOB07090705030C080607090E 
0607140COA00020404060307C20103020304050506050A050304020508050908030DOBOB090EOD1 1 
060 A0706060A060 70704 040605040 1060304050207030203070 50 30 3050706050C0A0C0B1 5 17OC0E 
20 07030C04080705040307030506020601040204020601050204000601050E090B0905140F1 1161911 20 

01080805030605040503020803080203030204020004030505050807080300040E11 151827282B1E 
06050602060306050404020303 030306020004040404020 l 02030406050909 1008 13 16232F2E41 IF 
060209050204020409060 103C40503060 10204030 10304090303070 507080D0B060F152B46634324 
02C40 10204020 1050 10402 01020 10 10 104080406020203040609050 50207090B0C 14 163 A48404C42 
15 000304040 10402 0202060704030 504020 10404010 606060 l 0303050 503080C0C0B0E l 5 34425 A4A IF 15 

010 305020301020404010 10301020403030404060303050904040207040608070800 15 1B3A473E26 
0203040 105020402020304 030505030601 O5O5OIO3O2O9O7O3O4O4O0O6O4O4O6OAOFOF14222E2 A 19 
000202000 101070502010301060305020 1040502010402030402020 80408080508090F08101F1313 
0 1C400030 1 050 10202020 104040 1 030502060303040 10503030503070 7030605060AOF 14080E 100B 
10 02CS050203C 0020 103030 10302000 10001010301020 1050402040309 060702040707090COBODOCOC 10 

0203010103030401050501 01030702010 104Q2040205040207030302030502080908070AOCOAOD07 
01000202 020204040 00 104040 1000302050402030304040202020 1020501 02070 A0408040B 09000 A 
C 30 3010 002C20 509040403 CO 000 50202020402020400020202 0407040404040604050AO 606050404 
0003020203C004020 1010100020204040302020202010604050 ib2COOS046lo'203050A08b60FOB09 
05 OOCCCOCCO 101020 1030205050 1030 1070404040203000 10202030407050706020206030407080706 05 

FF05C4060 1030 1010 10202010002010 10202020102060301 040 1020 30208030606060A0506030606 
FFC 10 1020203020 1000001 C00200020 l 030506010303050402010603030304020205060506060703 
FF02020000C0030 1030200020201020307020204030203050301030203040606050604040406070A 
01 FFC2010303040002C100000 1030 100010 101020201030201 02000304020204030604060B06050802 01 

1 5 10 15 20 25 30 35 40 


Figure 4-1. Typical solar raw raster scan (11 November 1967; 0202-0207 UT). 

The most sensitive mode of the experiment was used for this scan. 
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changed to improve the statistical fit, because of the deteriora- 
tion in the detector efficiency. Two print sheets are again 
required for each average raster scan. The UT interval in both 
date and day number appears across the top of each sheet. 

The interval extends through one exact 24 hour day, or 7 day 
week, with no overlap. 

Rasters are rejected by the program for the following reasons. 

The first part of the first raster scanned after the spacecraft 
reacquires the sun from a night side pass is unreliable. There- 
fore, these "sunrise rasters" are rejected from both the cali- 
bration raster and sun raster samples. Also questionable and 
missing data words are excluded from the sample. 

A calibration raster is rejected if at least the first 20 columns 
are not present, if the average count per word exceeds 5 in the 
first 20 columns (this criteria excludes rasters obtained when 
the spacecraft is passing through the South Atlantic Anomaly - 
see Section 4. 8), or if there are more than 100 interspersed 
missing and questionable data words. Also if the count per 
word anywhere on the raster exceeds 2 0, that word is rejected. 

Only sun rasters obtained with the instrument in the 0. 002 inch 
Beryllium filter and one arc-minute aperture modes are used to 
calculate N. n .In addition to the calibration raster rejection 
criteria, the last 2 columns of a sun raster are removed if there 
are less than 40 complete columns (i. e., less than 1920 words). 
This is necessary because of anomalous noise at the end of such 
"short rasters". Data words are rejected if their count on sun 
rasters exceeds 10. 

The number of rejected rasters and words and the number actually 

used in the reduction are listed on the first print sheet. Also 

listed are the values and statistical errors for N.,. , N._ ,, 

40 sun 40 cal 
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and N . The raster scan itself is printed out in the same 

Cell 

format as the raw raster, with the exception that each word is 
given in decimal notation, the first digit in units and the second 
in tenths. Figure 4-2 shows a typical daily average background 
raster. 

4.2.3 Corrected Raster Generation Program 
This program combines the information contained on the Raw 
Raster, A-E, and Background Raster tapes for the appropriate 
period, to produce a "corrected raster" or final x-ray filtergram 
of the sun. Again two print sheets are required for each raster 
scan. The attitude sheet contains the same data as presented 
on the equivalent raw raster sheet, with the addition of the value 
of N ^ used for the requisite average background raster. 

The corrected raster program uses the following algorithm: 

N. = N. - B. (4. 1) 

ill 

where N. = the true count rate of a word i; 

i 

N. = the raw count rate of word i from the raw raster; 

l 

B i = the background count rate from the same position 
on the background raster. 

A word is not printed on the corrected raster if either or both 
of the following conditions hold: 

N. <1. 0 

i 

N.<2/b7 (4.2) 

Words are suppressed in this manner to make x-ray emitting 
regions stand out for ease in analysis, and to minimize noise 
counts. A printed word consists of a two digit decimal in log 

•base 2, i. e., log N, . Thus the data word 30 is equivalent 

3 0 ^ 1 

to 2 or 8. 0 decimal. 
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TIME INTERVAL FROM NOV 11.1967 0 0 


0000 TO NOV 12.1967 0 0 


0000 


DAY OF YEAR FROM 


AVERAGE BACKGROUND RASTER 

1 5 10 15 20 25 30 35 40 

48 15191418192220212123221125232126202521191925272427183228333131373231403131333233 48 

17181417222017221721191722182117172425251822252722242826292629302928293640413844 
16182018252318212418212318182025211824242225252425302627312629273632302736313440 
45 23201816201815202117193122241818242123212622203127292632342227323731223835374841 45 

17141718202015242323212019172626192427282928322736213032313332273133343935423949 
16161722181718202122172118222326302626312623292923283027243027322729373942413745 
17 18 17 1818202 1182C2522 172222 182026252827202827262732262535273831 322532383 1444241 
12162214172024242521202223302425212628263127242529313130372828373435323543344743 
40 18131823161823172018182220262125192326302722302333273128333133373740343244494743 40 

18231722202223211924262026242822231919222523282629262930373529353241383535464840 
13122122212122-312528231820233023252034303028232928272825273239313141373839404241 
17181821212122232116232626222525232529302829302927342938353235353732404643424554 
22162020242220202325252725263130312626262628273227343734253625353035333637394147 
35 16212119212926252623273024212622272826252329263030353837343230363538493846454652 35 

182720221718272128222223232427202323293026283238323734323136313536353938475O5360 
17122217212521252926232628262825282926302938303226373236293938353241384747385453 
19202224232328232624222325281622342727313724232934313329393641373739425051455454 
23192018132726272026202325232032272730253336293631322831272934424135544047515857 
30 17151619252927222222252723262428272531332931282431323130323433423740435155545266 30 

22131719192322192627231825252827212630273032283033312831283934334444543852656164 
24122017242222222322252130262328282626323127353634303238393837403647455154595659 
17172316232525182221252930272225322728282832363334353435363137484755594552595753 
23222 5 171 92921 23 2C2 122 2323 2530252 527332727203 13236303433424 3 33354749546463626666 
25 21181920202421222630262727242530332230252724322533343032343443444857525359617771 25 

24181820172220232226252322192128253033272129263430273134444137494545535658596272 
16231721172422181618232223233132322823323028323233333640344136465744586566587772 
15201720212322202225262022273124263127273035312833263236324843505355556267656876 
15171624192323282026242422262526312732282829252934253841344144445050595954586775 
20 20201527162923272121212927252826292522332831283223313835343738464650585158677069 20 

1625I523171924222C2020 301 82422282524302728223232313036374041 52474455596267656079 
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Figure 4-2 Typical daily average background raster {11 November, 0 UT, 
to 12 November 19 67, 0* 1 UT). 
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Typical corrected raster scans appear as Figures 4-3 to 4-8. 
These 6 scans show the development of flares in 2 active regions 
for the time period included in Figure 5~4. Figure 4-7 was 
produced by combining its equivalent raw raster (Figure 4-1) 
with the daily average background raster covering the same time 
period (Figure 4-2). 

4. 3 Fulfillment of NSSDC Requirements 

Raw raster and background raster data tapes have been submitted 
to the National Space Science Data Center (NSSDC) as required 
in the Statement of Work. A card deck of the corrected raster 
program (see Appendix) was also included. The tapes are 
complete for the period 25 October 1967 to 12 May 1968 with the 
exception of the gaps noted in Section 4. 1. 

4.4 Active Region Ring Integration Routines 
4.4.1 CalComp Plotting Routines 

A major effort of the reduction program at AS&E has involved 
the generation of x-ray intensity vs time plots of individual' 
emitting regions, followed for periods of up to 60% of a single 
solar rotation. Several regions have been followed for more than 
one rotation. 32 regions have been plotted for the time period 
covered by this report. A list of the regions followed and the 
periods covered are given in Table II. Most have been analyzed 
in detail and the results appear in Section 5. 0. 

The plots are produced in the following manner. An active region 
is selected for study for various reasons, mainly based on its 
overall intensity and level of activity over an extended period 
of time. Once the region of interest is selected, its approximate 
center is located in arbitrary raster coordinates on the corrected 
rasters. The coordinates of the center are then listed as the 
region moves across the solar disc. 
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00 


TABLE II 

ACTIVE REGION TIME HISTORIES LIST 


Me Math Approximate 

Plage No . Period Followed 

(19 67) 


9034 

28 Oct. , 

9* 1 


1 Nov. , 

9 

9047 

27 Oct., 

6 h 


12 Nov. , 

19 

9062 

10 Nov. , 

°h 

- 

16 Nov. , 

21 

9073 

1 0 Nov. , 

°h 

- 

27 Nov. , 

11 

9082 

17 Nov. , 

16 h 


4 Dec. , 

0 

9088 

24 Nov. , 

15 h 

- 

4 Dec. , 

0 

9091 

24 Nov. , 

15 h 

- 

4 Dec. , 

0 

9101 

1 Dec. , 

15 h 


14 Dec. , 

14 

9107 

5 Dec. , 

n h 

- 

6 Dec. , 

21 

9108 

2 Dec. , 

n h 

- 

17 Dec. , 

0 

9110 

8 Dec. , 

16 h 


12 Dec. , 

8 

9115 

8 Dec. , 

16 h 

- 

23 Dec. , 

23 

9118 

12 Dec. , 

5 h 

- 

17 Dec. , 

0 

9128 

25 Dec. , 

5 h 


31 Dec. , 

6 


h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

h 


McMath Approximate 

Plage No, Period Followed 

(1968) 


9132/3 

2 Jan., 23^ 

- 

6 Jan. , 6? 

9146 

2 Jan., 23^ 

- 

19 Jan., 12? 

9184 

26 Jan., 6* 

- 

1 Feb., 20^ 

9188 

29 Jan., 23? 

- 

11 Feb., 5? 

9204 

8 Feb., 11? 

- 

22 Feb., 23^ 

9222 

18 Feb., 8^ 
21 Feb., 0\! 

- 

28 Feb., 8^ 

9224 

- 

28 Feb., 9? 

9225 

22 Feb. , 18^ 


28 Feb., 9^ 

9267 

10 Mar. , 21? 

- 

22 Mar. , 6^ 

9273 

23 Mar. , 3^ 


1 Apr. , 6? 

9281 

21 Mar. , 10? 

- 

27 Mar. , 0* 

9285 

19 Mar., 2jJ 

- 

3 Apr. , 15? 

9286 

24 Mar. , O^ 

- 

6 Apr., 3? 

9313 

8 Apr. , 8? 

- 

23 Apr. , 2 " 

9337 

24 Apr. , 13? 

- 

30 Apr., lOjJ 

9358 

30 Apr. , 21? 

- 

9 May, 20^ 

9364 

27 Apr. , 10? 

- 

12 May, 23* 

9372 

3 May, 1 


10 May, li n 
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Figure 4-3 Corrected raster, or solar x-ray filtergram, for 10 November 1967, 
1557-1602 UT. Note the two active limb regions: McMath 
plages 9 073 (left) and 9047 (right). 
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CORRECTED RASTER SCAN STARTING AT UT NOV 10» 1967 UT TIME 19 8 31.08 
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Figure 4-4 Corrected raster for 10 November 1967, 1908-1913 UT. A flare 
in region 9 047 (right) is near its peak on this raster. 


4-10 



CORRECTED RASTER SCAN STARTING AT UT NOV 10, 1967 UT TIME 21 4 44.16 
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Figure 4-5 Corrected raster for 10 November 1967, 2105-2110 UT. The 
bright x-ray flare in region 9073 (left), evident in Figures 
5-3 and 5-4, is near its peak on this raster. 
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CORRECTED RASTER SCAN STARTING AT UT NOV 11, 1967 UT TIME 0 5 58.49 
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Figure 4-6 Corrected raster for 11 November 1967, 0006-0011 UT. 
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CORRECTED RASTER SCAN STARTING AT UT NOV 11, 
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Figure 4-7 Corrected raster for 11 November 1967, 02 02-0207 UT. 

Compare with Figure 4-1, the equivalent raw raster scan. 
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CORRECTED RASTER SCAN STARTING AT UT NOV 11, 1967 UT TIME 4 48 4.25 

01 1/2 MIL BERYLLIUM LAMBDA 1 = 2.00 LAMBDA 2= 13.00 KN = 8311.46 

11 40 MIL APERTURE EPSILON * .15 N BAR CAL * .75 
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Figure 4-8 Corrected raster for 11 November 19 67, 0448-0453 UT. 
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A region can generally be followed for as long as desired until 
it rotates around the limb of the sun. However, several factors 
limit our detection ability of these regions. The primary problem 
is the decay in detector sensitivity. Nearly half of the regions 
were studied during the first two months of the spacecraft's 
lifetime, when the sensitivity was highest, and thus statistical 
fluctuations in the count rate were lowest. The instrument 
aperture mode is an important detection factor, because the 
4 arc-minute aperture mode has 16 times the flux gathering 
power of the one arc-minute mode. The filter mode is of lesser 
importance for this process, since only 20.5 minutes will elapse 
between filtergrams taken through the most sensitive filter. 
Rasters obtained with high noise counts (see Section 4. 8) limit 
the number of rasters available for analysis. The greatest 
number of "noise rasters" are those acquired in the South Atlantic 
Anomaly; typically these rasters account for 15% of the total 
number of rasters obtained daily. 

The list of center coordinates and the duration for each co- 
ordinate for a given emitting region, are entered into the 
particular ring integration program desired. Three tapes are 
used with these programs for the period of interest: the Raw 
Raster and Background Raster tapes, and an empty tape on which 
is printed the program results for use on AS&E's CalComp plotter. 
Raw raster rejection criteria are the same as discussed in Section 
4. 2.2. 

The ring integrals are calculated in the following manner. The 
coordinate of the center word is defined as Ring 1. Ring 2 
consists of the 8 data words surrounding the center word. Ring 3 
the 16 data words in the next square ring outwards, and so forth. 
The total uncorrected ring integral is then formed by summing 
the contributions from each ring over as many rings as desired. 
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Generally the integral is taken to include 10 rings. Ring sums 
are calculated about the same central coordinate independently 
for the raw raster of interest and its corresponding daily average 
background paster. Then the background sum is subtracted from 
its corresponding raw raster sum for each ring sequentially 
outward from the center, and the remainder integrated. If at any 
time the remainder for any ring is < 0, the integration is truncated 
at the previous ring. The program is required to integrate to at 
least 2 rings, since occasionally the central word (Ring 1) 
remainder is < 0 but contributions can still exist in outer rings. 

Finally corrections are applied to the integral to account for the 
decay in detector sensitivity and the point source response 
function (Section 3.3.1). The results are printed out in a format 
that lists for each raster scan the date, the UT and day number of 
the beginning of the scan, the filter and aperture modes, the 
integral with the number of rings counted and its statistical error, 
and the corrected value with its error. 

The final output is generated on the CalComp plotter. These 
plots present the time history of an individual emitting region. 

The corrected counting rate integrals are plotted on a vertical 
log scale vs a linear time scale. Corrected integral values less 
than 100 counts/0. 14 seconds are not plotted, because of their 
low statistical significance. 

Three different programs determine which filter and aperture mode 
combinations are to be displayed on these plots. For all emitting 
region time histories, plots are produced which display all filter 
and aperture mode combinations with error bars. The upper half 
of Figure 5-4 is an example of this type of plot. A second 
program plots only those integrals obtained when the instrument 
is in its most sensitive filter mode (0. 0005 inch Beryllium). 
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Both aperture modes are used, but error bars are not printed. 

This program uses a linear time scale of 0. 5 days/inch, whereas 
the first program uses 1. 1 hours/inch. This scale reduction 
allows one to easily detect systematic long-term variations in 
the emitting region. 

Typical examples of this compressed time scale type of plot 
appear as Figures 4-9, 4-10, and 5-3. Typical error bars are 
shown on the right side of Figures 4-9 and 5-3. Figures 4-11 
through 4-26 are tracings of the important active regions 
followed during the period covered by this report. All of these 
tracings are drawn in the 0. 0005 inch Beryllium filter and 4 arc- 
minute mode. 

The 0. 0005 inch Beryllium filter plotting program was not used 
for emitting regions studied after the end of January, 1968. 

Instead another program was used, which plots all integrals 
acquired with the 0. 0005 inch Beryllium and 0. 00025 inch Mylar 
filters in the 4 arc- minute aperture mode only. This change was 
dictated by the reduction in detector sensitivity, which decreased 
the detectability of non-flaring emitting regions through the one 
arc-minute aperture. The Mylar filter transmits at least 70% of 
the incident flux of the thin Beryllium filter, and is therefore 
useful for filling in time gaps on the plots created by the 
elimination of the one arc-minute integrals. 

4. 4. 2 Raster Active Regions Summing Program 
This program performs ring integrations of up to 7 selected 
emitting regions on a single raster. The appropriate raw and 
background raster tapes are used witn the following input data: 
the day number and UT of the desired raster, the number of 
regions on that raster for which ring integrals are needed, and 
the center coordinates for each region. Any number of rasters 
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Figure 4-9. Time history of the 2.5-12$ x-ray emission of Me Math plage region 
9034 in late October 1967. The aperture mode is one arc-minute. 
Typical error bars appear on the right side. 
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Figure 4-11. Tracing of the time history of the 2.5 - 12 £ x-ray emission 
of McMath plage region 9101. 
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Figure 4-14. Tracing of the time history of the 2.5 - 12 $ x-ray emission of McMath plage 
region 9115. 
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Figure 4-16A. Tracing of the time history of the 2.5 - 12 % x-ray emission of McMath plag 
region 9146. 
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Figure 4-1 7A. Tracing of time history of t 
plage region 92 04. 
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Figure 4-20. Tracing of the time history of the 2.5 - 12 & x-ray emission of McMath 
plage region 9225. 
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Figure 4-2 3A. Tracing of the time history of t 
region 92 85. 
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Tracing of the time history of the 2.5 - 12 A x-ray emission of 
McMath plage region 9313. 
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Figure 4-26. Tracing of the time history of the 2.5 - 12 % x-ray emission of Me Math plage 
region 9364. 



can be used within the time period of the raw raster tape. The 
ring integrations are performed in the same manner as discussed 
in the previous Section. The program prints out the input data 
and the same information as in the previous Section, with the 
exception that each raster scan will have a number of ring 
integrals listed under it. 

This program is used to calculate flux values that can be directly- 

compared with solar spectroheliograms published daily by the 

17 

Radio Astronomy Institute of Stanford University. The scans, 
taken between 2000 and 2100 UT, give radio brightness tempera- 
tures at about the same resolution as our data, and at a wave- 
length of 9.1 cm (see Figure 5-1). This analysis has been 
completed up to 12 February 1968 and the results are discussed 
in Section 5. 0. 

4. 5 Average Raster Program 

This program is useful for delineating the boundaries of emitting 
regions and improving the statistical precision of the data for 
faint active regions. The raw rasters obtained with a specific 
filter and aperture mode are summed over the particular period 
of interest. The counting rates for this sample are averaged 
and the corresponding background raster is subtracted point by 
point across the scan. The same rejection criteria as used in the 
background raster generation program (see Section 4.2.2), are 
utilized in this procedure. Rasters obtained in the South Atlantic 
Anomaly are listed with the input data and systematically removed 
from the sample by the program. 

The print-out format is similar to that of the corrected raster 
generation program (see Section 4. 2. 3), except that 1. 0 is added 
to the value of the average corrected counting rate to allow 
printing of fractional average count rates in logarithmic format. 
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Thus the output printed word represents 

log 2 (N +1.0) (4.3) 

The information generated by this program is stored on tape 
for future reference. The analysis is complete from 25 October 
to 26 November 1967. 

4. 6 Raster Orientation Program 

This program calculates the projected angles, a and 6, which 
are used to determine the orientation of the solar x-ray disc 
on a raster, a is the projected angle between the solar north pole 
and the pole of the ecliptic plane, the projection plane being that 
of the raster. For any given date the magnitude of a can be 
determined from the formula: 

cos a = 005 (4.4) 

cos B 0 

where i = the constant angle, 7° 15', between the solar and ecliptic 
- poles, and 

B q = the heliographic latitude. 

The sign of a is taken to be + between about 8 September to 
6 March for any year, and - the rest of the year. 

6 is the angle relating the pole of the spacecraft’ s spin axis 
with the solar north pole, as projected on the raster plane. 

The spin axis pole is always horizontally left of the center 
of the solar image on a raster. Angle 8 is determined by combining 
angle a with the spacecraft roll angle at any particular time. 

This result must be subtracted from 270°, the angle measured 
counterclockwise between the spacecraft spin axis pole and 
the north point of the raster plane. Figure 4-27 shows these 
relations clearly. 

The necessary quantities for these calculations are taken from 
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Figure 4-27. Raster format showing the relation of the 
orientation angles. 
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the appropriate A-E tape by the program, and the following A-E 
data is printed out: the spacecraft orbit number, the day number, 
the orbit start and end times (UT), the times of spacecraft 
entrance into and exit from sunlight, the spacecraft pitch 
and toll angles, and the computed angles a and 8. This program 
is essentially complete for the period covered by this report. 

Thus for any raster scan, the position of the solar north pole 
can be quickly located, knowing the angle 8 and the exact 
center of the solar image on the raster. The exact center of the 
solar image is displaced 0. 6 arc-minutes to the right and 0. 4 
arc-minutes up from the raster center. 

4. 7 Contour Plotting Program 

Techniques for producing smooth contours of the OSO-IV filter- 
gram data have been pursued. A typical contour map is shown 
in the center of Figure 5-1. This plot was produced by hand 
from an averaged filtergram of 3 corrected rasters obtained on 
25 October 1967 in the 0. 0005 inch Beryllium, 4 arc-minute 
aperture mode. An outline of the solar photospheric disc is 
superimposed on the x-ray map. The contour lines result from 
smoothing counting rates from adjacent words around areas of 
similar intensity. The values associated with given contour 
intervals represent true corrected decimal counting rates. 

We expect to produce such a contour map for each day from 
25 October 1967 to 12 May 1968 by use of a computer program, 
to be written in the next phase of data reduction. The OSO-IV 
x-ray data for this period then can be directly compared with 
that of other experimenters. 

4.8 Noise Suppression Criteria 

Various forms of data rejection criteria were discussed in Section 
4.2.2. The problem of noise and the methods used for suppressing 
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it on the output filtergrams are discussed in this Section. 

The Van Allen Belt of energetic charged particles dips far into 
the ionosphere over the South Atlantic Ocean to produce the 
South Atlantic Anomaly. The spacecraft's orbit intersects this 
extensive particle cloud several times a day. During these 
passes when our instrument is operating, the detector is saturated 
regardless of its mode of operation. Thus, these rasters are 
easily identifiable by their high counting rates, and can be 
readily eliminated from the sample. Figure 4-28 shows a typical 
corrected raster scan obtained when the spacecraft is just 
entering the Anomaly. 

The following methods are used to eliminate South Atlantic 
Anomaly rasters. In the background raster program, (Section 
4. 2.2), if the average count per word in the first 2 0 columns 
of a raster exceeds 5 that scan is removed. In the raster 
averaging program (Section 4. 5), a raster is rejected if the sum 
of the scan columns 10, 20, 30, and 40 exceeds 4000. This 
same criteria is used to identify South Atlantic Anomaly rasters 
in the ring integration routines (Section 4.4), but the rasters 
are processed normally with the words "Van A. Raster" printed 
out beside the counting rate. This is done to avoid the 
possibility of accidentally excluding flares with this criteria. 

During several orbits per day, raster noise appears abnormally 
high. These noise counts, along with the more usual random 
noise counts, are eliminated by the criteria discussed in 
Section 4. 2. 2. 

Analysis has indicated that the wavelength response functions 
shown in Figure 3-1 lead to anomalously high short wavelength 
x-ray fluxes in comparison with- the results of other observers. 

The most likely explanation for this discrepancy is a decrease 
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Figure 4-28 Typical corrected raster obtained when the spacecraft was 
just entering the South Atlantic Anomaly. 
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in orbit of the long-wavelength response' function of the 
detector photocathode. Recalculations of the response functions 
have been made in order to compensate for this apparent loss 
in long- wavelength efficiency. 
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5.0 SCIENTIFIC DATA ANALYSIS 


5. 1 General 

The OSO-IV x-ray telescope experiment was the first instrument 
with the capability of identifying the sources of solar soft x-ray 
emission, and of monitoring individual sources for periods on 
the order of the lifetime of the source. As a result, the instru- 
ment provided a unique source of information on the x-ray 
activity of solar active regions during the time period covered 
by this report. 

The analysis of the OSO-IV telescope results has followed 
three main lines. First, the angular resolution characteristics 
of the instrument have been used to determine the spatial 
dimensions of x-ray emitting regions. Second, time histories 
of the x-ray activity of individual, selected active regions have 
been studied in order to determine the characteristics of the 
time behavior of these regions. Third, the data has been com- 
pared with observations made at radio and optical wavelengths. 
The results of this analysis are contained in this Section. 

5. 2 Location and Size of X-ray Emitting Regions 

5.2.1 X-ray Emitting Regions 

A typical raster scan of the sun (Figure 5-1) shows eight 
distinct regions as sources of solar x-ray emission. The 9. 1 cm 
and CaK spectroheliograms of the same date are presented for 
comparison. As was originally noted by Chubb, et al. , 
the x-ray emitting regions are associated with solar active ' 
regions. Therefore, the x-ray brightness of the coronal limb 
in the instrumental passband is less than 150 photons/cm -sec- 
(arc-min) 2 . 
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Figure 5-1. An x-ray contour map of the sun produced 25 October 1967 by the OSO-IV x-ray telescope expert 
ment. Also shown are 9.1 cm and CaK maps for the same date. 



There is a substantial range in the x-ray emission of individual 

active regions. The total photon flux observed from region 9034, 

4 

the brightest region shown in Figure 5-1, is (4. 6 + .2) x 10 
2 

photons/cm -sec. The flux observed from region 9035 is only 
3 2 

(3.8 + .8) x 10 photons/cm -sec. Thus the range in x-ray 
flux in the detector passband is at least a factor of 12. More- 
over, several of the active regions observed in CaK are below 

the threshold of sensitivity of the x-ray experiment at this time 

3 2 

(approximately 2x10 photons/cm -sec). The physical con- 
ditions in the various active regions must differ for this to 'be 

the case. The ratio of the areas observed in CaK for regions 

17 

9034 and 9035 is only 2.2 ; therefore unless the vertical struc- 
tures of these active regions are very different, either the 
temperatures or densities of the x-ray emitting portions of these 
two active regions must be different. 

5.2,2 Vertical Extent of Active Regions 

8 9 

Recent high resolution solar x-ray images ’ obtained with a 
rocket-borne telescope indicate that the x-ray emitting regions 
extend to substantial heights in the corona. The resolution of 
the OSO-IV telescope is insufficient for a direct measurement 
of the vertical extent of an active region. However, the vertical 
extent of an active region can be estimated by observing the 
time at which the region becomes completely occulted by the 
solar limb due to rotation. The vertical extent is then calculated 
geometrically. In general, the time at which the x-ray intensity 
from the active region drops below the threshold sensitivity of 
the x-ray telescope is known to + 20. 5 minutes. The un- 
certainty is never larger than 95. 8 minutes. The vertical extent 
of the emitting region is calculated from the heliographic co- 
ordinates of the center of the underlying CaK region and the solar 
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2 1 8 

rotation rate w = 13.39 - 2.7 sin X (degree s/day) , where X 
is the heliographic latitude. Closely spaced regions of similar 
heliographic latitude, but different heliographic longitudes, 
are difficult to distinguish at the limb. Consequently, we have 
eliminated such cases from consideration for this analysis. 

Using data obtained between October 1967 and 1 January 1968, 
we were able to determine the time of first appearance at the 
east limb of seven emitting active regions and the time of 
disappearance at the west limb of eight. Measured vertical 
extensions ranged from 12, 000 km to 193, 000 km. The mean 
height of the top of x-ray emitting active regions was 105, 000 km. 
The results are tabulated in Table III. The accuracy of the 
proceedure can be estimated from the vertical extents computed 
for a single region, independently at the east and west limbs. 

The difference in height observed for region 9073 at the east 
and west limbs is 12, 000 km. 

Two major sources of error are present in these estimates. 

First, the x-ray output of active regions is highly variable. 
Therefore, the possibility exists that the x-ray emission of the 
active region could drop below the threshold of detectability 
because of time variations before the point of geometrical 
occultation is reached. This would result in an underestimate 
of the height of the active region. We have attempted to 

t 

minimize the possibility of this error by considering only the 
brightest x-ray emitting regions. Second, we have assumed 
that the top of the x-ray emitting region is located directly 
above the center of the CaK plage. An uncertainty in the height 
estimate proportional to the longitudinal extent of the active 
region is thus introduced. 
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TABLE III 

VERTICAL EXTENT OF ACTIVE REGIONS 


Me Math 

Appearance 

Time 

Disappear- 

Time 

Latitude 

Appearance 

CMP to 

0 

Height 

CaK Area 

Plage No. 

Date (1967) 

(UT) 

ance Date 

(UT) 


to CMP 

Disappear- 

(Degrees) (lCr km) 

(CMP) 




(1967) 



(days) 

ance (days) 




9034 



1 Nov 

0839 

N12 


8.3 

20. 16 

43 

3700 

9041 



6 Nov 

1116 

N19 


9. 1 

29.24 

102 

2700 

9047 



13 Nov 

0258 

S21 


9. 6 

35.22 ' 

157 

(62 00) 

9048 


, 

11 Nov 

1442 

N23 


8. 6 

21.61 

53 

2500 

9066 

8 Nov 

0319 



N2 0 

7. 7 


10. 67 

12 

2500 

9067 

8 Nov 

0319 



N18 

9. 0 


28.19 

94 

3200 

9073 

10 Nov 

1735 



N14 

9.2 


31. 73 

123 

(7800) 

9073 



29 Nov 

0345 

N14 


9. 3 

33.06 

135 

(7800) 

9088 
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S19 
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2 Dec 
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2 6 Dec 
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38.35 
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7500 
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12 Dec 

0633 



N21 
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36.52 

171 

9100 

9120 

14 Dec 

0740 



S2 7 

9.4 


30.64 

114 

4300 


*Millionths of a Solar Hemisphere 



In Table III we have listed the CaK area of the active regions 
as observed near central meridian passage. If we assume that 
the outline of the base of the x-ray emitting region conforms, 
in general, to the outlines of the CaK plage {as observed by 

g 

Vaiana et al. , that the topology of all x-ray emitting regions • 
are similar, and that the region does not change in size between 
the center of the disk and the limb, then the vertical extent of 
a region should be proportional to the square root of its CaK 
area. It is evident that this relationship does not hold. It 
would therefore appear that the shapes of the x-ray emitting 
portions of active regions differ from region to region. 

5.2.3 Size of Flares 

One can estimate the size of x-ray flares observed on the disk 

by comparing the observed x-ray intensity profiles with the 

expected response of the telescope to sources of various sizes. 

We have examined the intensity profiles of 42 bright x-ray 

flares observed between 27 October 1967 and 3 December 1967. 

For each of these events the diameter of the telescope field 

stop was one arc-minute. We have compared these intensity 

profiles with those calculated from the response of the telescope 

for circular sources of various diameters. Figure 5-2 is a 

histogram of the relative frequency of occurrence of the measured 

full widths at half maximum intensity (FWHM) . The mean value 

of the FWHM is 1.51 arc-minutes, which corresponds to a 

source diameter of 1. 24 arc-minutes. The standard deviation 

of the distribution is 0. 30 arc-minutes. We have compared the 

measured FWHM of the x-ray intensity distributions with the 

corresponding values of the measured area in Ha for the 24 x-ray 

events from this sample which were listed as H a events in the 

17 

ESSA Solar Geophysical Data Bulletin. There is no apparent 
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FWHM OF FLARE INTENSITY DISTRIBUTION (arc-min. ) 

Figure 5-2. Histogram of the relative frequency of occurrence of the measured 
full widths at half maximum (FWHM). 



relationship between the measured area in Hq- and the observed 
width of the x-ray distribution. If the jitter in the OSO-IV 
pointing control (about 0. 25 arc-min) is taken into account, the 
observed x-ray intensity distributions are consistent' with source 
dimensions of one arc-minute or less. 

In a similar manner we have examined the intensity profiles of 
10 bright x-ray flares observed at the limb. The mean value of 
the observed FWHM is 1. 64 arc-minutes which corresponds to a 
source diameter of 1.35 arc-minutes. The standard deviation of 
the observed distribution of limb flare FWHM is 0. 38 arc-minutes. 
Therefore, the vertical extent of the x-ray flares observed by the 
OSO-IV telescope would also appear to be on the order of one-arc 
minute. 

From the OSO-IV results, it would appear that the majority of 

the soft x-ray flux from flares is emitted from regions whose 

characteristic dimensions are on the order of one arc-minute. 

We can therefore set an upper limit on the volume of the x-ray 

29 3 

emitting region for the observed flares of about 10 cm . If 

2 

we then assume an emission measure, N e V, on the order of 

48 -3 19 

10 cm as reported by Culhane and Phillips for the flares 

observed in the 1 to 3 $ region of the spectrum, the electron 

density in the x-ray emitting region of the flare must be greater 

9 5 —3 20 

than 10 * cm . Gabriel and Jordan have derived electron 

9 11—3 

densities between 10 and 10 cm for non-flaring active 
regions. The higher densities could be achieved if the flare 
regions were filamentary in structure, or if the flare volumes 
were substantially smaller than the OSO-IV limiting resolution. 

4 
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5. 3 Time Variations in X-ray Emission From Active Regions 

5.3.1 The Range of Variation in the X-ray Emission of an 
Active Region 

The magnitude and variety of the changes in the instantaneous 
x-ray emission observed from an individual active region can be 
noted by examining a plot of the time history of such a region 
over a period of several days (Figure 5-3). The most notable 
feature of this plot is the extreme range of variability in 
instantaneous x-ray counting rate. The total variation between 
the highest counting rate shown (15. 6 minutes after the x-ray 
peak of the importance 3b flare of 2130 UT on 16 November 1967)- 
and the lowest rate shown (0513 UT on 10 November) is 160. 

In order to extrapolate back to the peak of the x-ray flare, we 
use the fact that the x-ray intensity of the 16 November flare 
fell exponentially with a time constant of 17. 5 minutes. This 
would indicate a maximum range in the 2. 5 to 12 2 x-ray intensity 
of a factor of 400 in a seven day period. This range is typical 
of that exhibited by solar active regions during periods of flare 
activity. For the period 2 7 October 1967 to 6 January 1968, the 
maximum range of x-ray flux variation exhibited by an individual 
active region during a passage across the disk was shown by the 
region corresponding to McMath-Hulbert Plage No. 9091. This 

3 

region was below the threshold of detectability, 3. 8 x 10 
2 

photons/cm -sec until 24 November 19 67. On December 19 67, 

6 

0652 UT, we observed a flare with a peak x-ray flux of 6.5 x 10 
2 

photons/cm -sec. Thus, the range in x-ray flux was at least a 
factor of 1700. 

Even if flares are not considered, the x-ray counting rate for an 
individual active region varies considerably during a single 
passage across the disk. The non-flare x-ray counting rate of 
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Figure 5-3. Time history of the 2.5 - 12 $ x-ray emission of McMath plage region 9073 for 
7 days in November 1967. Also shown are the x-ray spectral hardness (the 
ratio of emission at 2. 5 - 9 $ to emission at 2. 5 - 12 A in percent) and the peak 
9. 1 cm radio brightness temperature of the region. 



region 9073, for example, (shown in Figure 5-3) rose by a 
factor of twenty in a period of four hours and then slowly fell 
back to approximately twice its original value over a period of 
four days. Variations of a factor of ten in x-ray counting rate over 
a disk passage are typical. 

5.3.2 Time Behavior of Active Regions 

Intense, impulsive x-ray bursts are immediately obvious on the 

plot of the time history of the x-ray intensity of an active region 

(Figure 5-3). Such events appear as a sequence of a few points 

with substantially higher x-ray flux than the mean level of 

emission from the active region immediately before or after. 

Enhancements of a factor of ten in the x-ray photon flux from an 

active region during an impulsive burst have been observed. 

Enhancements of a factor of two or three are typical. Because 

the interval between comparable measurements (20. 5 minutes) 

is of the same order of magnitude as the duration of these 

2 1 

impulsive events {see, for example, the data of Teske ), 
the probability of observing the peak of an event is relatively 
low. Consequently, the flare enhancement in x-ray emission 
observed by this experiment is, in general, less than the actual 
enhancement. 

In addition to the impulsive x-ray events, the time history of an 
individual active region reveals more gradual variations in x-ray 
intensity. The durations of these merge into the longer impulsive 
x-ray flares. These events may indeed reflect the pre-and post- 
flare state of the active region. The longer events presumably 
contribute to the slowly varying component of the solar x-ray 
emission. 

In general, the time profiles determined from the OSO-IV data 
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are affected by the presence of unresolved rapid fluctuations In 

x-ray brightness. Figure 5-4 illustrates this point. The time 

history of the x-ray emission in the three OSO-IV wavebands is 

shown for two active regions. Also shown are higher time 

resolution, 2 to 12 S total solar x-ray flux measurements for this 
22 

period . One effect of the low time resolution of the OSO-IV 
instrument is to increase the apparent duration of the impulsive 
x-ray flares. However, the x-ray images allow the identification 
with only small ambiguity of the active regions responsible for 
the impulsive bursts. Moreover, by comparison with high time 
resolution data, it is easy to distinguish those portions of the 
OSO-IV data which are affected by the impulsive x-ray bursts. 

By means of comparisons between the OSO-IV data and higher 
time resolution total solar x-ray flux measurements, we have 
been able to distinguish periods in which the time history of a 
particular active region is unaffected by unresolvable fluctations. 

For convenience in analysis, we have chosen to divide the x-ray 
event arbitrarily into three classes on the basis of their durations. 
Events of duration less than 8 hours are termed "flares". Events 
of duration greater than 24 hours are called "slow variations". 
Events with intermediate duration are termed "long enduring 
brightenings " . 

The establishment of three separate classes of events is quite 
arbitrary. It is difficult to distinguish between the shorter 
duration members of the intermediate class and impulsive events. 

It is sometimes equally difficult to distinguish between long 
enduring brightenings and slow variations. In general, however, 
gradual brightenings appear to be superposed on slow variations 
and flares on gradual brightenings (as in Figure 5-3). 
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Figure 5-4. 


Solar activity from 10 November, 


08 h UT, 


to 11 November 


1967, 09 a UT. The top half of the graph shows the x-ray 
activity from two active regions as observed on OSO-IV. 
The bottom half shows tracings of data from Explorer XXXV 
(2 - 12 £) and 2800 MHz radio flux. 
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The relationship between long enduring brightenings of the active 
region and the occurrence of large x-ray flares tends to justify 
our classification. We have selected the 12 impulsive x-ray 
events with peak-observed flux in the 2. 5 to 12 $ waveband 

g 

(0. 0005 inch Beryllium filter) greater than or equal to 5 x 10 
2 

photons/ cm -sec which occurred during the period 27 October 
to 3 December 1967. Three of these events occurred at times 
when the non- flare level of x-ray emission from the active region 
was below the threshold of sensitivity of our detector. The time 
behavior of the active region before and after the 9 remaining 
events was examined. In all 9 cases, the large flare occurred 
during the course of an x-ray brightening of duration between 
6 and 18 hours. 

We have also examined 16 easily identifiable long enduring 

x-ray brightenings. Of these, only 9 produced major x-ray 

6 2 

flares with peak flux greater than 5x10 photons/cm -sec. 

In fact, no impulsive x-ray events were observed at all in the 
case of 3 of these gradual x-ray brightenings. An examination 
of the fast time resolution records of Explorer XXXV^ i n the 
2 to 12 £ band and OSO-III in the 1 to 1. 6 band^ confirm 
these results. Therefore, we can conclude that although 
the existence of a gradual x-ray brightening seems to be 
required for the occurrence of large impulsive x-ray flares, the 
existence of a gradual x-ray brightening does not necessarily 
imply that such a flare or indeed any x-ray flares will occur. 

In Table IV, we list some of the observed characteristics of the 
three event classes. The quantities listed are the extremes of 
the range of measurements made on 49 flares, 16 gradual events, 
and 7 slow variations. The ranges of duration of the three 
classes are defined. The peak fluxes are calculated from the 
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TABLE IV 

TYPICAL CHARACTERISTICS OF THREE EVENT CLASSES: MEASURED QUANTITIES 



Flares 

Long Enduring 
Enhancements 

Slow Variations 

Duration 

< 5 min. - > 1 hr. 

8-24 hrs. 

2-5 days 

Peak Flux 

10~ 3 - io“ 2 

10" 3 -5 x 10~ 3 

CO 

1 

0 

i-H 

1 

1 

o 

i— i 


ergs/cm 2 -sec. 

ergs/cm 2 -sec. 

ergs/cm 2 -sec. 

Total Emission 
(2.5 - 12 X) 

2 

1-10 ergs/cm 

10-10 2 ergs/cm 2 

2 

10-10 ergs/cm 

Rise Time 

<20 min. 

1-5 hrs. 

1/2-4 days 

Decay Time 

<20 min. -1 hr. 

2-6 hrs. 

1/2-4 days 

Spectral Hardness 
Index* 

>50 

40 

35 


*Spectral Hardness Index is the relative x-ray intensity observed in two wavebands. 
It is a measure of the effective temperature of the emitting regions. 



peak x-ray counting rate extrapolated from the observed time 

behavior of the event. The conversion from counting rate to 

energy units is performed by assuming an exponential spectral 

shape characterized by the measured spectral hardness index 

(defined below). The total output in the detector passband 

observed at 1 a. u. is obtained by numerical integration of the 

event time profile. The range of total output values is listed 

by order of magnitude. The rise and decay times are the time 

constants of exponential fits to the event counting rate profiles 

observed through the 2. 65 jigm/cm Beryllium filter. The peak 

spectral hardness index is the ratio of the counting rate 

2 

measured through the 9.4 yqvci/cva Beryllium filter to the rate 

2 

measured through the 2. 65 y gra/cm Beryllium filter (in percent) 
at the peak of the event. The spectral hardness index increases 
as the x-ray spectrum becomes enhanced in shorter wavelength 
photons. 

Examination of Table IV shows that the total energy emitted at 
in the detector passband for each of the three event classes 
is roughly comparable. The lower peak fluxes of the longer time 
scale events are compensated by their longer durations. The 
peak spectral hardnesses for the event classes decrease as the 
duration increases, implying higher temperatures for the more 
impulsive events, but there is substantial overlap. 

Since the total output in the detector passband from each type 
of event is comparable, the majority of the energy emitted by 
an active region will be emitted during the course of the most 
numerous events. Thus it would appear that the greatest portion 
of the total observed output in the detector passband would come 
from flares. In fact, this proves to be the case. We have 
integrated over time to find the total energy emitted in the 
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detector passband by the active region corresponding to McMath 

plage 9073. Over the eighteen days that this region was 

30 

observed, it emitted a total of 1. 7 + 0. 5 x 10 ergs in the 

2.5 - 12 $ band. The importance 3b flare of 16 November 1967 

29 

alone emitted 4.2 + 0. 8 x 10 ergs of this energy. 

5.3.3 Spectral Behavior of Active Regions 

We have attempted to fit the rising and falling portions of x-ray 
events to exponential. curves. This is done independently for 
the three different x-ray filters. In general, the exponential rise 
and decay times are found to be shorter for the same event 
observed through the thicker filters than through the thinnest. 

For example, in the case of the event of 1630 UT, 10 November 1967 

o 

(Figure 5-4), the rise time as observed through the 2. 65 ^gm/cm 
Beryllium filter (2. 5 to 12 $) is 73 + 5 minutes. The rise time 
measured for the 9. 4 n gm/cm 2 Beryllium filter data (2. 5 to 9 $) 
is 51 + 7 minutes. Similarly, the decay times are 325 +_ 25 
minutes for the thinner filter and 260 _+ 20 minutes for the thicker 
filter. Thus, the x-ray spectrum becomes enhanced in higher 
energy photons as the event rises toward maximum intensity.- 
During the decay phase, the spectrum softens. This behavior 
has been noted for 9 cases in which there is a statistically 
significant difference in the characteristic times for the two 
Beryllium filters. In 13 other cases, there was no significant 
difference in the characteristic times. No events have been 
measured in which the spectrum softens towards the peak of the 
event. 

In Figure 5-5 we plot the relationship between the counting rates 
measured through the 0. 0005 inch Beryllium filter and the 
spectral hardness index for the bright active region shown in 
Figure 5-3. The relationship between spectral hardness and 
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2.5 -12 A X-RAY COUNTING RATE (counts/.l4) 



SPECTRAL HARDNESS INDEX 

Figure 5-5. Relationship between the 2.5 - 12 & 

x-ray emission and the spectral hardness 
index for McMath region 9073. 


5-18 



x-ray counting rate is monotonic within the experimental errors. 

If we assume that 2. 5 - 12 X x-rays are produced by thermal 
mechanisms,- the spectral hardness index will be related to the 
temperature of the emitting plasma. The exact form of this 
relationship depends on the shape of the x-ray spectrum as a 
function of temperature including the effect of emission lines. 

The x-ray counting rate, on the other hand, is a function not 
only of the temperature, but also of the emission measure. 

Thus the monotonic relationship between the counting rate and 
the spectral hardness index shown in Figure 5-5, implies that 
for the six day period plotted in Figure 5-3, either the emission 
measure of this active region did not change significantly, 
or changes in emission measure are proportional to temperature. 

If an exponential shape is assumed for the x-ray spectrum, 
a limit of a factor of 3’ can be put on possible changes in the 
emission measure of the region. 

Observations of impulsive x-ray flares at shorter wavelengths 

19 23 

by Culhane and Phillips and by Hudson et al. have shown 
that the emission measure remains constant during the decaying 
portion of such events. Here, over a much longer time scale 
we observe a relatively constant emission measure in a single 
active region. However, there are significant differences in 
emission measure from region to region. For example, during 
the period on 10 November 1967, shown in Figure 5-4, the 
x-ray counting rates observed from region 9047 were a factor of 
five higher than the rates observed from region 9073 at the same 
value of the spectral hardness index. 

5. 4 Relationship Between X-rays and Other Solar Emissions 
5.4.1 X-ray and Ho- Flares 

In general, individual x-ray impulsive events are associated with 
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individual H a flares. Teske and Thomas have shown that, 
with only one possible exception, all Hey flares of importance 
one or greater are accompanied by an impulsive soft x-ray event. 
Most, if not all, subflares are also accompanied by detectable 
x-ray emission. 

The reverse correlation is not quite as obvious. Many x-ray 

events occur in active regions which are at or beyond the limb 

in Hey. Thus the Hey event, if any, accompanying a particular 

x-ray event may be unobservable. We have selected a set of 

x-ray events for which there was a high probability that any 

accompanying Hey event could be observed. The center of 

intensity of the x-ray event was required to be within 75° of the 

central meridian. The x-ray intensity was required to be at 
5 2 

least 2.5 x 10 photons/cm -sec. 

Between 27 October and 3 December 1967, 35 impulsive x-ray 
events were observed which satisfied these conditions. A 
search was then made of the solar flare listings of the ESSA Solar 
Geophysical Data Bulletin. An x-ray event was considered to 
have an Hey counterpart if an Ho- flare or subflare was found in 
the active region producing the x-ray event with the time of H a 
maximum occurring within 10 minutes of the x-ray maximum. 

Hey counterparts were identified for 32 of these 35 x-ray events. 

Hey events varied in importance from f to 3b. In Figure 5-6 we 
have plotted the peak observed x-ray intensity as a function of 
Hey importance classification for those events observed with the 
2.5 to 12 $ filter. Although there may be a tendency for large 
x-ray events to be associated with Hey flares of greater importance, 
a large range of Hey importance classifications can be covered 
at moderate x-ray intensities. 
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Figure 5-6. Ho- importance clissification of solar flares as a function of peak x-ray photon 
flux. 



It would appear that virtually all x-ray events with intensity at 

1 a.u. in the detector passband greater than 2.5 x 10^ photons/ 

2 

cm -sec are accompanied by Hq- flares. However, these results 
show that the magnitudes of Hq- and x-ray emission associated 
with a particular event are not well correlated. 

5.4.2 X-ray Active Regions and the Slowly Varying Component 
of the Radio Flux 

The general correlation between solar soft x-ray emission and 

25 

microwave emission is well known. We have already 

mentioned the correspondence in location between the x-ray 

emitting regions and the sources of the 9. 1 cm radio flux. 

2 1 6 

Several observers ’ have noticed the general correlation 
between the peak radio brightness temperature, T^, of a solar 
active region and its x-ray flux. 

We have compared the 9. 1 cm peak brightness temperature of 

active regions measured by the Stanford Radio Astronomy 
17 

Institute with the value of the x-ray photon flux measured 
by the OSO-IV telescope at the same time. Figure 5-7 shows 
this comparison. Measurements made during impulsive x-ray 

flares have been omitted. There is a strong tendency for high 

/ 

x-ray flux to be accompanied by high T^; however, there is a 
significant scatter in this relationship. One could conjecture 
that this scatter is due to the variations between active regions. 
On the other hand, it is easily shown that even for a single 
active region, the relationship between 2.5 to 12 .8 photon flux 
and 9. 1 cm is not monotonic. Figure- 5-8 shows the tra- 
jectories on the N , T, plot of individual regions of high activity. 

X D 

In Figure 5-8 A the x-ray flux actually decreases as the radio 
brightness temperature increases. The trajectory of another 
active region (McMath 9091) is shown in Figure 5-8 B. Here, 
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Figure 5-7. Comparison of x-ray photon flux and 9. 1 cm peak 

brightness temperature for individual active regions. 
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9.1 CM Tb (10° °K) 


Figure 5-8. Time behavior of the x-ray photon flux 

from individual active regions as a function 
of 9. 1 cm radio brightness temperature. 

The trajectories of two active regions are 
5-24 shown. 



in general, the x-ray counting rate is proportional to the radio 
brightness temperature, however, a large increase in x-ray 
emission is accompanied by only a small change in radio 
brightness temperature. 

Although the behavior of these two regions is different, both 
types of behavior may be understood, at least qualitatively, 
in terms of a model in which both x-rays and microwaves are 
emitted thermally from the same plasma. 

During the period between 10 November 1967 and 15 November 

1967 spectral hardness data was obtained for region 9073. 

This data implies in the context of a thermal model that the 

emission integral of the region remained constant while the 

electron temperature decreased. If the emitting region is 

considered optically thin to 9. 1 cm radiation and if this 

radiation is emitted by the thermal bremsstrahlung process, 

the radio brightness temperature should increase as the x-ray 

-1/2 25 

intensity decreases, T b ~T . On the other hand, the 
generally linear relationship between and the x-ray 
intensity of region 9091 might indicate that the temperature 
of this region remained constant while the emission measure 
increased. The jump in x-ray intensity at constant radio 
brightness temperature might indicate an increase in temperature, 
and the ensuing decrease in both T^ and x-ray intensity might 
indicate a drop in emission measure at constant temperature. 
Unfortunately, spectral hardness information was not obtained 
during this period and so the temperature behavior of the region 
is unknown. 

The behavior observed for region 9073 is inconsistent with a 
model in which the region is optically thick to 10 cm radiation. 
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If the region were optically thick, the observed decrease in both 
the spectral hardness and x-ray flux would have been accompanied 
by a decrease rather than an increase in radio brightness 
temperature. 

5.4.3 X-ray Events and Microwave Bursts 
24262728 

Several observers ’ ’ * have studied the relationship 

between soft x-ray flux increases and microwave radio bursts. 
Detailed correlation for any but the largest events has been 
hindered by the lack of spatial resolution and also the con- 
comitant low signal-to-noise ratios. This problem is alleviated 
by imaging x-ray detectors, such as that of OSO-IV, in that 
the region responsible for the x-ray emission can be identified 
and analyzed independent of the x-ray flux from other active 
regions . 

As a preliminary step, we have established the statistical 

relationship between reported 10 cm radio bursts and the soft 

x-ray events observed by the OSO-IV telescope. During the 

period 27 October to 3 December 1967, 62 distinct outstanding 

events of all classes at 2 700 and 2800 MHz were reported by 

29 

Covington, Gagnon and Moore. For this analysis, we have 
considered complex events as single events. Fifteen events 
occurred while the spacecraft was either within the earth's 
shadow or passing through regions of high trapped particle 
flux. X-ray brightenings of at least two standard deviations in 
the flux from an individual active region were observed in con- 
junction with 43 of the remaining 47 outstanding events at 
10 cm. Three of the four radio bursts that were not observed 
to be x-ray events were of duration less than the 307 seconds 

necessary to complete a solar scan. None of the radio bursts 

-22 -2 

which were unobserved in x-rays exceeded 4x10 watts-m 
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-Hz * in peak, flux at 2800 MHz. Thus, of the 39 centimeter 

wave radio bursts of duration longer than five minutes, 38 

produced noticeable x-ray brightenings. Similarly, 37 of the 

-22 

38 centimeter radio bursts with peak flux greater than 3x10 
-2 -1 

watts -m -Hz produced noticeable x-ray brightenings regard- 
less of duration. We therefore consider it to be reasonably 

certain that any radio event of longer duration than 5 minutes 

-22 -2 -1 

and greater peak flux than 4 x 10 watts -m -Hz , has an 
observable x-ray counterpart. 

30 

Kawabata has pointed out the statistical correlation between 
long enduring microwave radio bursts (both gradual rise and fall 
bursts and post burst increases) and soft x-ray flux enhance- 
ments. We have examined the relationship between a specific 
long enduring microwave burst at 2800 MHz and the corres- 
ponding x-ray event. Analysis of this type is necessarily re- 
stricted to those events whose durations are long compared to 
the 20. 5 minute time resolution of the OSO-IV instrument. 

During the period 1650 UT to 2030 UT on 10 November 1967 

2 9 

(shown in Figure 5-4), Covington, Gagnon and Moore 
distinguish two separate radio events at 2800 MHz. One is a 
complex event which we associate with the impulsive x-ray 
bursts observed by the solar x-ray flux monitors. Most of these 
Impulsive bursts can be attributed to region 9 047. The second 
microwave event is a gradual rise and fall radio burst (simple 
3 A) which we associate with the x-ray brightening of region 9073. 

In Figure 5-9 we show the time histories of the 2800 MHz 
31 

simple 3A burst and the x-ray brightening of region 9073. If 
the association of the radio event with the x-ray event is correct, 
we see that the x-ray brightening begins slightly earlier than 
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Figure 5-9. Trajectory of the x-ray photon flux and 2800 MHz radio flux of 
and active region during the rise phase of a long enduring burst. 
The time (UT) at which a measurement was taken is listed next 
to the data point. Loci of constant temperature and constant 
emission measure are shown for a model in which both the x-ray 
5-28 and microwave optical depths are small. 


the radio burst and continues for much longer. The x-ray 

intensity rises to two standard deviations above the previous 

mean level at 1640 UT. The radio burst is listed as starting at 
29 

1650 UT. The duration of the radio burst is given as 220 
minutes. The duration of the x-ray event is at least 650 minutes. 
The time of the peak of the x-ray event is somewhat uncertain, 
but it is clear that the peak of the radio burst occurs 20 to 40 
minutes before the peak of the x-ray event. The presence of this 
substantial time delay would appear to argue against models 
in which the long enduring radio burst and long enduring x-ray 
burst are produced by the same mass of plasma by a thermal 
process. 
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6.0 CONCLUSIONS 


The OSO-IV pointed x-ray telescope experiment continues to 
function satisfactorily. Since the failure of the second tape ■ 
recorder on OSO-IV on 12 May 19 68, only data transmitted 
during real-time passes over the receiving stations has been 
available. Recently the amount of quick-look data received 
from GSFC has been greatly reduced, because of the telemetry 
overload on the Stadan network. At the time of tape recorder 
failure, the detector efficiency of the OSO-IV telescope was 
about 30% of its value on 2 7 October 1967. Since then it has 
stabilized at somewhat less than 2 0% of this value. 

The instrument has been operated almost entirely in its most 
sensitive mode (0. 0005 inch Beryllium filter; 4 arc-minute 
aperture) since tape recorder failure. Occasionally it is 
commanded into the calibration mode for reference. Despite 
the decrease in sensitivity, the instrument can still detect all 
major flares that occur on the sun during our real-time passes. 

As mentioned earlier in this report, information obtained with the 
OSO-IV x-ray telescope has been exchanged with other ex- 
perimenters, including those from the University of Iowa, 
University College of London, Algonquin and Dominion Astro- 
physical Radio Observatories, the Stanford Radio Astronomy 
Institute, the University of Hawaii, the University of Maryland 
and the University of California at San Diego. Results of the 
data analysis have been presented to the scientific community 
at various meetings. Abstracts or reports are available for some 
of these meetings^ 5 Others include two papers presented 

at the Midwest Cosmic Ray Conference in Iowa City, Iowa, in 
March, 1968: Paolini, F. R.; Vaiana, G. S. ; Giacconi, R. ; 
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Reidy, W. P. and Zehnpfennig, T. , " Spectroheliograms of X-ray 
Flares from OSO-IV," and Vaiana, G. S.; Paolini, F. R. ; 

Giacconi, R. ; Reidy, W. P. and Zehnpfennig, T. , "X-ray Plage 
Studies from OSO-IV. " A paper by Vaiana, G. and Krieger, A. 
entitled "Results from the OSO-IV Solar X-ray Telescope (3-13 $)" 
was given at the Fiftieth Annual Meeting of the AGU from 21-25 
April 1969 at Washington, D. C. A presentation entitled "Results 
from the OSO-IV Soft X-ray Telescope Experiment" by Giacconi, R. ; 
Vaiana, G. and Krieger, A. was made at the Orbiting Solar Ob- 
servatory Workshop at Boulder, Colorado, from 4-15 August 1969. 
A paper is in preparation reviewing the results of the OSO-IV 
telescope experiment as presented in this report. 
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7.0 OSO-IV: STELLAR X-RAY WHEEL EXPERIMENT 


The Wheel Experiment, which was designed to measure cosmic 
x-rays, malfunctioned early in the flight. Also during the short 
time that it operated, it was apparent that the sensitivity of the 
primary instrument, the photoelectric detector (PED) was well 
below that determined during the final calibration, possibly as a 
result of degradation of the electron-multiplier between the time 
of the final calibration and the launch. 

The malfunction in the form of high counting rates was observed 
in the PED beginning sometime between 17 and 18 November 
1967. This was apparently caused by a malfunction of our 
instrument, but it was never possible to isolate the nature or the 
origin of the malfunction. 

There was no indication that the high counting rate was caused 
by factors associated with the spacecraft, and there was no 
indication that the spacecraft was adversely affected by the 
malfunction. The malfunction was characterized as follows: 

1. There was a high counting rate (100 counts per second) 
induced in the PED. The onset of this high rate was not neces- 
sarily coincident with the turn-on of the instrument at night. The 
time delay varied between a few seconds to half an orbit. Follow- 
ing its onset, the high counting rate persisted with some fluctua- 
tion until turn-off. 

2. The Nal detector in its high gain state was also being 
affected, although the increase was not as large as it was in the 
PED. 

3. The remaining detectors, the anthracene and the anti- 
coincidence, were not affected. 

The various housekeeping functions incorporated into the instru- 
ment did not show any significant variation since launch. A plot 
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of the values of the monitors for a number of orbits is shown in 
Figure 7-1. The small changes that occurred were consistent 
with drifts in the analogue circuitry and in the A/D converter. 

Thus, the malfunction was not a result of a power supply failure, 
nor was it adversely affecting the power supply. Figure 7-2 is a 
plot of the onset time of the high counting rate in units of ASC 
frames following turn-on through Orbit 700. There appears to be 
significant differences in the onset time, particularly between 
Orbits 530 and 62 0. However, it was not yet clear what the 
significance was. 

Plotted in Figure 7-3, 7-4, 7-5 and 7-6 are the counting rates 
around the orbit for the several detectors for a series of four . 
consecutive orbits. Figure 7-3 shows the PED count rates. For 
Orbits 619 and 62 0 the count rate was low and apparently normal 
for about half the night. Figure 7-4 shows the count rate in the 
anthracene detector. A large peak is observed at around 0° longi- 
tude indicating passage through the South Atlantic Anomaly. Else- 
where the anthracene detector is apparently counting normally; in 
particular, no abrupt change is noted in coincidence with the 
onset of the high count rate in the PED. Figure 7-5 shows the count 
rate in the plastic anticoincidence detector. Again the South Atlan- 
tic Anomaly is apparent at 0° longitude. In addition, variations 
around the orbit are noted, which can be correlated with changes in 
the cosmic ray fluxes associated with the changing magnetic lati- 
tude. Again, no increase is observed coincident with the onset 

of the high rate in the PED. Figure 7-6 shows the count rate in the 

15 

Nal detector. As noted in our earlier report, the Nal detector ex- 
hibits a more complex counting rate profile around the orbit than any 
of the other detectors, which simply reflects the fact that the device 
is sensitive to a greater variety of radiation than the others. 

The South Atlantic Anomaly is apparently a variation dependent 
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on magnetic latitude. There is also an indication of an increase 
in count rate coincident with the increase that occurs in the PED, 
but of a much smaller magnitude. 

The conclusions from these data are that the malfunction was 
causing extraneous counts to be thrown into the PED channel and 
to a lesser extent into the high gain state of the Nal detector, 
but that these detectors were still functioning, both before and 
after the onset of the high rate. As noted, the remaining detectors 
and the lower gain states of the Nal detector were apparently not 
being affected. The implication of this malfunction, as far as the 
experiment is concerned, can be summarized as follows. Little 
or no useful data could be obtained from the PED, which is the 
primary x-ray detector during the high count rate phase. The 
signal-to-noise ratio in the high gain state of the Nal detector is 
being adversely affected by perhaps a factor of 2. Under these 
conditions, a significant number of the primary scientific objec- 
tives of the experiment could not be met. 

We attempted to understand the origin of the malfunction without 
success. The basis for considering it to be a malfunction and not 
simply an increase in count rate due to radiation is primarily the 
failure to connect the observed rates with' anything physically 
real. There was no apparent variation of the count rate in the 
PED with filter wheel position, nor are there any systematic 
variations with time or around the orbit- as might be expected to 
occur with charged particles; neither was there an up-down 
asymmetry of the radiation. In addition, the sudden onset of the 
high rate is difficult to understand as a real increase due .to 
radiation. The appearance is that of some component that was 
being overs tressed either electrically' or perhaps thermally and 
suddenly becoming very noisy. An example of this might be 
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voltage breakdown of a capacitor. No malfunction of this kind 
has been observed during the entire OSO Wheel Instrument test 
and qualification program. The only malfunctions that occurred 
have been catastrophic failures induced by high voltage discharge 
during thermal vacuum tests. On at least two occasions, we did 
experience high count rates induced by RF interference. This 
occurred only when the arms were in a down position and affected 
only the high gain state of the Nal detector. Thus, we did not 
feel that this malfunction was related to an RF problem. 
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Figure 7-2 
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APPENDIX 



DOS PL/I COMPILER 360N-PL-464 CL3-5 


RAWRAST 


04/24/70 


/* AS CE OSO-D PROGRAMS // GET RAW RASTERS */ 


I 

? 

3 

4 


5 


6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
23 
26 
29 

31 

32 

33 

34 

35 

36 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


/* ASSE OSO-O PROGRAMS // GET RAW RASTERS */ 

RWST • • PROCEDURE OPTIONS (MAIN) , • 

DECLARE ILINEPIC PICTURE *99’, ILINEOUT CHARAC TER ( 2 ) 

DEFINED ILINEPIC, COLID CHARACTER ( 80 ), . 

DECLARE RECID POINTER, . 

DECLARE OUTAPE FILE RECORD SEQUENTIAL BUFFERED OUTPUT 

ENVIRONMENT (V(3980) MED I UM ( SYS008 , 2400) NOLABEL LEAVE 
BUFFERS ( 2 ) ) ; 

DECLARE 01 RASTER EXTERNAL ALIGNED, 2(UTIME, UTIM05) FLOAT B INARY ( 53 ) t 
02 (ISIZE, ISUBS l f IBADCT, ISKPED, IMISS, IENDN, IEOF, 

* JFWEEL , JAWEEL) FIXED B INARY (31), 02 ( I RASRY ( l 920 ) , 
ISC05RU0), I SC07R ( 10 ) , ISC27RC10), ISC44R(10)> BIT<16),. 
OECLARE (DATE, SDATE , YRDATE ) FLOAT BINARY ( 53 ) , . 

DECLARE ARRAY (27) , ISUBSCR FIXED B INARY ( 3 i) , KAY(40) CHARACTER^),. 
DECLARE VOLTRAY (4) PICTURE *99*, VCR AY ( 4 ) CHAR AC TER ( 2 ) 

DEFINED VOLTRAY, • 

DECLARE FNAHE CHARACTER ( 24) , ANAME CHARACT ER ( 20 ) , . 

DECLARE MONAME CHARACTER ( 4) , SGN CHARACTER ( L ) , . 

DECLARE (IFPRNTC, I FT APEC ) CHARACTER! t ) , ( l F TAPE , IFPRNT) BIT{1>,. 
DECLARE ISBEGN BIT(l),. 

DECLARE MESSAGE CHAR ACTER ( 50 ) ,. 

DECLARE (NTRAST , NXRAST , REWEPH, GETEPH, REWFIL, RUNF1L, SDAY r RUNEPH, 
DEGOMS, FLTLUK , APTLUK, OSODLP, RUNTL ) ENTRY; 

INT2 = 2 , • 

COLID = f 1 5 10 15 20 25 30 

35 40*,. 

CALL NTRAST , . 

UTO = 39490, • 

GET EDIT ( UTBEGN , UTEND, IFPRNTC, IFTAPEC) (2F(I5,5), 2 A ( L ) } , * 
IF IFTAPEC = * • THEN I FT APE= ■ 0* B, . ELSE DO, . IFTAPE = M*B f . 

CALL REWFIL (OUTAPE),. OPEN FILE (OUTAPE),. END,. 

IF IFPRNTC = * 1 THEN IFPRNT - *0*B,. ELSE DO,. IFPRNT - 
CALL REWEPH,. END,. 

ISBEGN = f 0* B , * 

NEXTRASTER..- CALL NXRAST,. 

IF IEOF NE 0 THEN GO TO ENDF I LERAS TERS , . 

IF UTIME LT UTBEGN THEN GO TO NEXTRASTER,. 

IF UTIME GT UTEND THEN GO TO HAXT I MEXCEEDEO , . 

IF NOT ISBEGN THEN DO,. ISBEGN = *1*6,. UTLOTAPE = UTIME, .END,. 
IF IFTAPE THEN WRITE FILE (OUTAPE) FROM (RASTER),. 

AJTHIGHTAPE = UTIME, . 

IF NOT IFPRNT THEN GO TO NEXTRASTER,. 

SDATE ■= UTIME + UTO,. 

CALL SDAY (SDATE, [YEAR, YRDATE, MONTH, DATE, MONAMF ) , . 

IDATEI = DATE,. 

HRS = 24. * (DATE - IDATEI) , . 

CALL DEGDMS (HRS, ISGN, I HR ,' MIN, SECS),. 

CALL FLTLUK (JFWEEL, FNAME , QLAM1 , QLAM2, FNX, FEX),. 

► CALL APTLUK (JAWEEL, ANAME),. 

/* GET ASPECT-EPHEMERIS DATA */ 

CALL GETEPH (UTIME, ARRAY, SPRATE, ASPANG, ISDAY ) , • 
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/* A$£E OSO-D PROGRAMS // GET RAW RASTERS */ 


51 

52 

53 

54 

55 

56 

57 

58 


59 


60 
6 1 
62 

63 

64 


65 

66 


67 


68 

69 


70 

71 

72 

73 


DELT = OSOOLP (UTIME, ARRAY),. 

CALL NEWPAGE,. 

SIDTIM ^ SIDGRN (UTIME + UTO) f . 

/* COMPUTE AND PRINT EPHEMER I S QUANTITIES */ 

CELONG * ATAN ( ARRAY (2 ) , ARRAY ( 1 ) ) , . 

TRLONG = (SIDTIM - CELONG) * 57.29578,. 

TRLAT = ARRAY ( 22 ) * 57.29578,. 

CLONGD = CELONG * 57.29578,. 

PUT EDIT { ' EPHEMER I S AT START CF SCAN*, 'CELESTIAL LONGITUDE', 
CLONGD, • DEGREES 1 , 'POSITION VECTOR 1 , ARRAY ( 1 ) , ' KM*) 

{ SK IP I 2 ) , A, X < 13 ) , A, F ( 9 ,2 ) , A, XI15), A, FU0,1>, A),. 
PUT EDIT ( ' T ERRESTRI AL LONGITUDE', TRLONG, ARRAY { 2 ) , 

•GEODETIC LATITUDE', TRLAT, ARRAY ( 3 ) , 

•GEODETIC ALTITUDE', ARRAY (23)) (SKIP, X(37>, A, F(9,2), 

X ( 38 ) , F< 1 0, 1 ) , SKIP, X ( 4 1 ) , A, F{9,2), X(38), F(10,1>, 
SKIP, X ( 41 ) , At FI8 v lt I «• 

RA = ATAN J ARRAY { 20 ) , ARRAY ( 19) ) * 57.29578,. 

ROLLA = ARRAY ( 27 ) * 57.29578,. 

DEC = ASIN ( ARRAY (21 ) ) * 57.29578,. 

ALPHA = ARRAY ( 26 ) * 57.29578,. 

PUT EDIT ('ASPECT AT START OF SCAN', 'PITCH ANGLE', ALPHA, 

• DEGREES', 'SPIN AXIS', ARRAY{ 19) , 'ROLL ANGLE', ROLLA, 
'UNIT VECTOR', ARRAY ( 20 ) , 'RIGHT ASCENSION OF SPIN AXIS', 
RA, ARR AY< 21 ) , 'DECLINATION OF SPIN AXIS', DEC) (SKIP(2), A 
, X ( 24 ) , A, F ( 9, 2 ) , A, X(22), A, Fill, 7), SKIP, X{48), A, 

F ( 9 , 2 ) , X ( 2 8 ) , A, F ( 1 1 , 7 ) , SKIP, X(30), A, Ft9,2), X(39), 

F ( 1 1 , 7) , SKIP, X ( 34 ) , A, F(9,2)) f . 

SIDTIMD = SIDTIM * 57.29578,. 

PUT EDIT ('SIDEREAL TIME AT START OF SCAN', SIDTIMO, • DEGREES' 
, 'SOLAR', ARRAY ( 7) , 'UNIT VECTOR', ARR AY ( 8 ) , ARRAY ( 9) ) 
(SKIP(2>, A, F ( 9 , 3 ) , A, X(54), A, F(li,7), SKIP, X(95), A, 

F ( 11 , 7) , SKIP, X ( 106) , F(U,7)),. 

/* PRINT STATISTICS */ 

PUT EDIT (ISKPED, 'RASTER SCANS LOST DUE TO MISSING DATA', 
IS12E, 'WORDS IN THIS RASTER SCAN', IBADCT, 

•QUESTIONABLE DATA WORDS IN THIS RASTER SCAN', IMISS, 
'MISSING DATA WORDS IN THIS RASTER SCAN') (SKIP, F(5,0), 

X( 1) , A) , . 

/* PRINT OUT SUBCOMMUTATOR VOLTAGE READINGS */ 

IF ISUB SI = 0 THEN GO TO NOSUBCOMS,. 

PUT EDIT ('VOLTAGES FROM SAIL ANALOG SUBCOMMUTATOR 1 , ISUBSZ, 

• SETS', '05) 6 VOLT MONITOR', '07) 6 VOLT MONITOR', 

•27) NEGATIVE 6 VOLT SUPPLY', '44) 2.4 KV MONITOR* ) 

( SKIP ( 2 ) , X ( 32 ) , A, X ( 30 ) , F(5,0), A, SKIP, X(L3>, A, 

X ( 12 ) , A, X ( 4 ) , A, X( 11 ) , A) 

DO 1 = 1 TG ISUBSZ,. 

VCRAY = ' • , . 

IF NOT I SC05R ( I ) AND 'OOOOOOOl'B THEN VOLTRAY ( 1 ) = SUBSTR 
f I SC05R ( I ) , 9, 8) , . 

IF NOT I SC07R ( I ) AND 'OOOOOOOl'B THEN VOLTRAY ( 2 ) = SUBSTR 
(ISC07RM), 9, 8),. 
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74 

75 

76 

77 

78 

79 


80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

L09 
110 
ill 
1 L2 

113 

114 


IF NOT I SC2 7R ( I ) AND •OOOOOOOL'B THEN VOL TRAY ( 3 ) = SUBSTR 
[ ISC27RI I) , 9, 8 } y • 

IF NOT l SC44R ( I ) AND *00000001*B THEN VOL TRAY ( 4) = SU3STR 


( ISC44RI I), 9, 8),. 

PUT EDIT (VCRAY) (SKIP, 4 {X( 28 J, A(2n>,. 

END,. 

CHECKRASTREND.. IF I ENDN=1 THEN PUT ED IT ( * LAST RASTER BEFORE SUNSET * ) 
(SKIP, A),. 

IF IENDN ^ 2 THEN PUT EDIT 

{ * LAST RASTER BEFORE TIME DISCONTINUITY IN DATA 1 ) (SKIP(2), 
A) » . 

CALL NEWPAGE, - 

PUT EDIT (COLID) (SKIPi'2), X(4), A(80)l f . 

DO M = 1 TO 48,. 

ILINE = 49 - M,. 

ILINEOUT = * 

IF I L I NE= 1 OR ILINE=48 OR ( IL I NE/5 ) *5= IL I NE THEN I LINEPIC 
- ILINE,. 

KAY = • * , . , 

00 J = 1 TO 39 BY 2,. 

ISU8SCR = 48 * I J-U + 49 - M , * 

CALL PROCWORD ( KAY I J ) ) , • 

ISU8SCR = 48 * J + M,. 

CALL PROCWORD ( KAYU+1)) 

END, . 

PUT EDIT (ILINEOUT, KAY, ILINEOUT) (SKIP, A(2), X(2)i 40AI2), 

X( 2) , A( 2) 

END,. 

PUT EDIT (COLID) ISKIP<2), X14), A(80)j,. 

GO TO NEXTRASTER,. 

NOSUBCOMS.. PUT EDIT ( 'NO SU6C0M VOLTAGES AVAILABLE*) ISKIP(?), A),. 

GO TO CHECKRASTREND, . 

ENDF1LERASTERS. . PUT EDIT ( * END OF FILE ON TELEMEIRY TAPE*) ( SK I p ( 2 ) , 

A) 

GO TO WRAPUP,. 

MAXTIMEXCEEDED.. PUT EDIT ( * MAX I MUM REQUESTED TIME EXCEEDED ON TELEMETR 
Y TAPE* ) ( SK IP ( 2 ) , A) f . 

WRAPUP.. CALL RUNTL,. 

IF IFPRNT THEN CALL RUNEPH,. 

IF NOT IFTAPE THEN STOP,. 

CLOSE FILE (OUTAPE),. 

CALL RUNFIL (OUTAPE),. 

MESSAGE = * * , . 

PUT STRING (MESSAGE) EDIT (* M OSO-D RAW RASTERS FROM*, UTLJTAPE, 

* TO*, UTHIGHTAPE, **••) (A, F(10,3), A, F<10,n>,. 

DISPLAY ('LABEL OUTPUT TAPE ON SYS009 AS FOLLOWS* ),. 

DISPLAY (MESSAGE),. 

STOP,. 

ASIN.. PROCEDURE (X),. 

RETURN (ATAN (X / SQRT <U ~ X * XJl),. 

END,. 
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115 

116 


117 


118 

119 

120 
121 
122 
123 
127 
129 

131 

132 

133 

134 

135 


NEWPAGE.. PROCEDURE,. 

PUT EDIT ( • RAW RASTER SCAN STARTING AT UT 1 , MONAME, 

IDATE1, I YEAR, • UT TIME', IHR, MIN, SECS, 

UTIME) (PAGE, A, X(2), A, F(3,0), A, FI5,0L, A, 2F(3,0), 
FI 6, 2 ) , X<10), FI 10 , 3) ) , * 

PUT EDIT ( FNAME, 'LAMBDA 1 QLAM1, 'LAMBDA 2=*, QLAM2, 

•KN = ', FNX , * /EPS I LON KE = ', FEX, '/EPSILON' , ANAME) 

{SKIP, A, X { 6 ) , A, F 1 6 , 2 ) , X{4), A, F<6,2), X(5), A, 

F ( 8 , 2 ) , A, F I 8 ,2 ) t A, SKIP, A},. 

RETURN,. 

END,. 

PROC WORD. • PROCEDURE (OUTSTR),. 

DECLARE N, OUTSTR CHARACTER t 2 ), . 

OUTSTR = • * ,. 

IF ISU6SCR GT ISIZE THEN DO,. OUTSTR = '==»,. RETURN,. END,. 

IF IRASRY (ISUBSCR) AND 'OOOOOOOl'B THEN DO,. OUTSTR= ■ ++ ' , . 
RETURN,. END,. 

N = SUBSTR I IRASRY ( ISUBSCR) , 9, 8),. 

CALL I NTHEX ( N, INT2, OUTSTR),. 

RETURN,. 

END,. 

END,. 
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*/ 


STMT IFVEL NFST 
1 

0 l 

1 I 

A 1 

4 1 

A l 

7 1 

R 1 

9 1 

10 1 

U l 


12 1 

11 l 


14 1 

15 1 

16 l 

17 1 

1 5 l 

1° 1 

20 1 

1 1 

2 1 

21 1 

24 l 

? 5 1 

2 6 1 

27 1 

2 8 ! 

20 1 

10 1 

11 1 

52 l 

52 l 

16 1 1 

15 1 

50 1 

40 1 

41 1 

42 1 

44 1 

44 1 


I* ASF.E nSO-O PROGRAMS // GFT BACKGROUND RASTERS */ 

BKGR : PROCEDURE OPTIONS (MAIN); 

OECIARF ILINFPIC PICTURE • 99 > , IL1NC0UT CHARACTER ( 2) 

DEFINED IL1NEPIC» COLIO CHARACTER I 80 ) ; 

DECLARF ( INDFXf I) FIXED 8 INARY (31); 

DECLARE PREMNS FIXED DEC I MAL I 3 , l ) ; 

DECLARE MESSAGE CHARACTER { 60 ) ; 

DECLARE SiJMTABl 1800) , SUMTA04 11 flCO ) , ( MSUNTAB (1 800 ) . MS JNTAB411 800 ) ) 
BIT18) ; 

DECLARE MTFMP FIXED BINARY(fi); 

DECLARE IF40MIL B f T( 1 ) ; 

DECLARE MRAY11920) FIXFD DEC I MAL 11*0) ; 

DECLARE RFC I D POINTER; 

DECLARE i OUTRECD (5 AS FD ( R EC ID ) , 2IEPS, UTSTRT!, UTSTDP1, XNE ARC AL » 
DFLCALi XN40CAL, DEL40C4L, XN400NE , DEL400NE, 

RAKR AY ( 1920 ) FIXFD DFC IMAL 1 1 , 1 » ) ; 

DECLARE 0UTF11 F TLF RFCORQ SEOUFNT I AL BUFFERED OUTPUT 
ENVIRONMENT (V(3884)); 

DECLARF 01 RASTER EXTERNAL ALIGNED, ? (UT l M E ,UT I M05 ) FLOAT BINAOY153), 
02 (ISiZE, ISUBS7. I BADCT , ISKPFD, IMISS, IENON, IFOF, 
JFWEEL , JAWEEL) FIXED BINARY 111"), 02 1 TRA SR Y ( 1970 I , 
ISC05R110), ISC07P110), ISC27R110), I SC44P { 10 ) ) BIT1I6); 
DECLARE (DATEA, SOATE, YROATE) FLOAT BINARY! 53) ; 

DFCL ARE ISBEGNBITU), ISBEGN1 RIT(l); 

DCCLAPE ( MflNAM 1 , M0NAM2 ) CHARACTER 1 4) , SGN CHARAC TER I 1 ) ; 

DECLARE KKAY140) CHARAC TF p ( 2 1 ; 

DECLARE NBAR FIXED B I NARY ( 8 ) ; 

DECLARE SIJMRAYC 19?0) ; 

DFCL ARE ISGN CHARACTER ( 2 ) ; 

DECLARE 1MHISUNWD, MH1SUNWD4, I40CAL) FIXED 8 I NAP Y I 1 1 , 0 ) ; 

DECLARE ( I YEAR l , 1 YEAR? , MONTH) FIXED BINARY! 17); 

ON ERROR SNAP SYSTFM; 

COLIO = 101 D5 10 16 20 to 

55 40'; , 

CALL NTRAST; 

UTORG = 291.; 

UTO = 19490.; 

ISBFGNl = *0*B; 

I SR EON = *0*R; 

O p FN F r LE (OUTFIL ) ; 

IfcNDN = o; 

GFT LIST 1UTSTRT, UTEND , TD IF ) ; 

IF UTEND <= UTSTRT | TDIF <= 0 THEN DO; PUT FDIT 

(MLLFGAl INPUT PARAMETERS') (SKIP, A); STOP; FNO; 

UTS TOP = UTSTRT + TDIF; 

CALL NITCOUNTS; 

NEXTRASTER: IFNDNB4 = IFNON; 

CALL NXRAST; 

TESTRASTER: TF IEOF -.= 0 THEN Gt) TO FNOOFILE; 

IF UTIME < UTSTRT THEN GO TO NEXTRASTER; 

IF UT I ME > UT S TOP THFN GO TO ENDT I MFI NTFRVAL : ‘ 
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/* ASSE OSO-D PROGRAMS // GET BACKGROUND RASTERS 


STMT 

48 
53 
55 
' 57 

59 

60 
6'+ 
65 

69 

70 

74 

75 

76 

77 
7ft 
80 
ftl 
82 

83 

84 
8* 

89 

90 

91 

92 

94 

95 

97 

98 
1 0 2 

103 

104 

105 

106 
107 
109 

112 , 
1 14 
115 
1 16 
117 
I 18 
119 
7 ?0 
121 
1?7 


*/ 


LEVFL NEST 


1 

l 

1 

1 

1 

1 

1 1 

1 

1 I 

1 

I ' l 

1 

1 

1 

1 l 

1 l 

1 1 

1 1 

1 1 

1 

1 1 

1 
1 
1 
l 
t 

1 l 

1 1 

1 1 

1 ? 

I 1 

1 1 

1 1 

1 
1 
1 

1 

1 

1 

1 

1 1 

1 1 

l l 

1 l 

i 1 

1 1 


IF ^ iSBEGNl THEN DO; ISBEGNl - 'l’B; UTSTR T=UTI ME 5 END; 

IF IFNDNB4 = 1 THEN GO TO NEXTRASTER ; 

IF IENPN = l THFN ISIZE = ISI7F - 90; 

IF JFWEFL -•= 3 THEN GO TO SUNRASTER; 

/* THIS IS A CALIBRATION RASTER */ 

MTOT * MTOT +1; 

IF [SIZE < 960 THEN DO; MSMALL = MSMALL + L; GO 

TO NEXTRASTER; END; 

IF IMTSS > 100 THEN DO; MMISS = MMISS + 1? GO 

TO NEXTRASTER; END; 

IF I8ADCT > 100 THEN DO; MHORR - MHORR + l; GO 

TO NEXTRASTER; END; 

/* DETERMINE AVERAGE COUNT DURING 1ST ?0 LINES */ 

IDFNOM = 0; 

SUM>0 = 0*; 

DO I - 5 TO 960; 

IF IRASRY(I) t *000000 1 1 • B THEN GO TO ENDLOOP; 

NBAR = SUBSTR( IRASRYC II # 9); 

S U* M 2 0 = SUM20 * NBAR; 

IDENOM = IDENOM + 1; 

ENDLOOP: END; 

IF SUM20 / IDENOM > 5 THFN DO; MNOISC = MNO I SC +■ i? GO 

TO NEXTRASTER; END; 

MMSTOT - MMSTOT + imiss; 
m * M + 1; 

MBADCT = MBADCT + IBADCT; 

IF [SIZE > 1920 THEN ISIZE = 1920; 

DO I = l TO ISIZF; 

IF I RASRY (II £ * 0000001 1 1 B THEN GO TO ENDLOOP, 1 ; 

NBAR = SUBSTR t IRA SRY ( l ) * 9); 

IF NBAR >= 20 THFN DO; MHIWD = MHIWO + 1? GO TO ENDL00P1 ; 
END; 

SUM RAY C I ) - SUMRAYU) + NBAR; 

MRAY(I) = MR AY ( I) + l; 

ENDLOOP l ; END; 

GO TO NEXTRASTER; 

SUNRASTER: IF JAWEEL > 4 THEN GO TO NEXTRASTER; 

IE JAWEFL = 1 | JAWEEL = 4 THEN IF40MIL = U*B; ELSE IF40MIL 
= f 0* B ; 

IF ISIZE < 96 THEN GO TO NEXTRASTER; 

MTEMP = o; 

TH1SAV ^ 0; 

DO J = 1 TO 19 BY ?; 

INDEX = 4ft * j; 

CALL TRF ATS UN WORD; 

INDEX * 48 * J + 1; 

CALL TREATSUNW0R9; 

INDEX = >48 ♦ ( J-i) * l; 

CALL treatsunword; 



STMT 

1 ?1 
l?4 
! ? 5 
1?6 
12 S 
1 10 
131 

m 

114 

ns 

11 * 

117 

110 

140 

141 

14? 

144 
] 4* 
141 

145 
1ST 

m 

i s? 

1 S 3 
1 54 
1 SO 
1 *0 
161 
1 4? 

1 61 
1 64 

165 

166 
167 
1 61 
160 
l 70 
171 
t 77 

173 

174 

1 75 

176 
1 77 


/« askf nsn-n programs // get background rasters */ 


LEVEL NEST 


1 l 

1 L 

1 1 

1 
1 
t 
1 
l 
1 
1 
1 
1 
1 
l 
1 

1 

1 

i 

i 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


index = 48 * < J + tl ; 

CALL TREATSUNWORO; 

END: 

IF MTEMP <- 0 THEN GO TO NEXTRASTER? 

[F 1 F 40 M IL THEN GO TO SUNRA STER 4 » 

MIOTOT = M 10 TDT + l; 

If M 10 T 0 T > 1800 THEN SIGNAL ERROR 1 
SUMTAO (M 10 T 0 T ) = THISAV; 

MSUNTAB(MIOTOT) = MTEMP; 

GO TO NEXTRASTER; 

S UN RASTFR 4 : M^OTOT = M 40 T 0 T + l; 

IF M 40 TOT > 1100 THEN SIGNAL ERROR ; 

SUMT AB 4 ( M 40 TDT ) = THISAV; 

MSUNTA 04 ( M 40 TOT ) = MTEMP; 

GO TO NFXTR ASTER; 

/* END nr TIME interval reached */ 

ENDTIMEINTFRVAL: IF ISBEGN THEN GO TO PR I NTSTA T I ST l C S ; 

TF MTOT 0 THEN GO TO PR I NTSTAT I ST I CS ; 

IF ISBEGNl THFN Gn TO FNDGROUP; 

I = U 1 T 1 ME - UTSTRT ) / TDIF; 

UTSTOP = UTSTRT +- ( I * 1 1 * TD t F 

UTSTRT - UTIME; 

ISBEGNl = * t 1 B * 

GO TO GETNEXTGROUP; 

/* PRINT OUT STATISTICS ON BACKGROUND RASTERS */ 

PR I MTSTAT I ST TCS S SDATF * UTSTRT + UTO; 

IP ISBEGN THFN DO; ISBEGN - f l f B; UTREGN = UTS 1 PT; END; 

ISBEGN = 'PR; 

CALL SDAY (SOATE, I YEAR V, YRDATF, MONTH , DATEA, MONA^H ) ; 

I DA f l = DATEA; 

HRS = 24 * * t DATEA-IDAT 1 ) ; 

CALL DEGPMS (HRS, ISGN, IHRl, MINI, SECS 1 ); 

SDATF = IJTSTOP + IJTO; 

CALL SDAY (SOATE, IYEA*?, YRDATF, «ONTH, DATEA, MONAM?); 

I DAT? = DATEA; 

HRS = 24 . 4 c (DATFA-IDAT 2 ) ; 

CALL DEGDMS (HRS, ISGN, I HR? , MIN 2 , SECS 2 ) ; 

CALL NFWPAGE; 

PUT FDIT (MTCT, ’CALIBRATION RASTERS’) (SKIP(L), R ( ST AT FO Q M ) ) ; 
STATFOrm: FORMAT ( SK IP ( l ) * FI 5 , 0 ), X(l),.A>; 

PUT EDIT (MSMALL, 'HADLFSS THAN ?0 LINES’) (R (STATFOPM) ) ; 

PUT FDIT (MMISS, 1 HAD MORF THAN 100 INTERSPERSED MISSING WORDS’ 
) (R(STATFORM) ) ; 

PUT FDIT (MHHRR, ’HAD MORE THAN LOO QUESTIONABLE WORDS’) 

( R ( STATF FIRM ) ) ; 

PUT EDIT (MNniSC, 1 HAD AVERAGF COUNT OF MORE THAN 5 TN 1 ST ?0 *L 
INFS * 1 (MSTATFORM) } ; 

PUT FDIT (M, 1 W F R F USFABLF'M ( R ( ST ATFORM ) ) ; 

PUT FDIT (MMSTOT, ’MISSING WORDS’, MB ADC T , • QUEST I UNABLE WORDS’ 
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/* ASCF OSn-D PROGRAMS // GET .BACKGROUND RASTERS 


*/ 


STMT LEVEL NEST 


1 7 8 l 

17*? 1 

181 1 

187 I 

183 1 

184 1 

IRS 1 

t 86 1 

1 37 1 

138 1 

189 1 

1 9n l 

191 1 

197 1 

193 1 

1 94 1 

195 1 

196 1 

701 1 

704 1 

705 l 

706 l 

707 1 

?03 1 

210 1 

211 1 

717 1 

21** 1 

214 1 

715 1 

216 1 

213 1 

719 1 

720 1 

7 24 1 

775 1 

224 1 

777 1 

7,23 1 

729 1 

730 1 

23! 1 

737 1 


, MHIWOt 1 WORDS REJECTED BECAUSE COUNT GREATER THAN 70M 
( R ( S TATFORM 10 ) ) ; 

ST ATFORM lOr FORM AT { SK T P < 1 ) , F<10,0), XU), A); 

IF M = 0 THEN GO TO ENDGROUP; 

LOCATE OUTRECD FILE(OUTFIL) SET(RECID); 

UTSTRT1 = UTSTRT; 

UTS TOPI = UTSTOP; 

XN40CAL = D. ; 

I40CAL = 0; 

DO J = 1 TO 19 BY 7; 

1 INDEX = 48 * J; 

1 CALL treatbackrword; 

1 INDEX - INDEX + l; 

1 CALL treatbackrword; 

1 INDEX = 48 * (J-i) + i; 

1 CALL TP fatbackrword; 

1 INDEX = 48 * ( J+l) ; 

I CALL TREATBACKRWORD; 

1 END; 

IF I40CAL O THEN 00; DEL40C4L = 0* ; XN40CAL = 0.; END ; 
ELSE DO; DEL40CAL = SQRT ( XN40CAL ) /I40CAL ; XN40CAL 
1 = XN40CAL / I40CAL; END; 

PUT EDIT <»N BAR 40 CAL = f , XN40CAL, f DEL 40 CAL = S DEL 40CAL ) 

( R { NB ARC AL FORM )') ; 

N8ARCALF0RM: FORMAT (SKIP(2), A(I9), F(12,4), XUO), A(19>, Hl?,4)l; 

CALL PROCSUN ( ' 10 ,f , M10T0T, MHISUNWD» SUMTAB, MSUNT4R ) ; 

IF ID6N0M -.= 0 THEN GO TO PR I NTNBAR400NE ; 

CALL PROCSUN { 1 40 •» M40TOT , MHISUNWD4, SUMTAB4, MSUNT AB4) ; 
PRINTN8AR400NE: PUT EDIT (»N BAR 40 ONE = •, XN400NE , f DEL 40 ONE 

DEL400NE ) ( R ( NB ARC AL FORM ) ) ; 

XNQARCAL = XN40CAL - XN400NE ; 

OELCAL = SQRT ( DFL40CAL * 0EL40CAL + DEL400NE*0EL400NE ) ; 

PREMNS = 0. ; 

DO I = T TO 1920; 

1 IF MRAY(I) > 0 THEN BAKRAYU) - SUMRAYII) / MRAY(I) - XN8ARCAL 

1 + PREMNS + .05 /* ROUND TO 0.1 */; ELSE BAKRAY(I) =0; 

I PREMNS = 0* ; 

1 IF BAKRAYM) < 0 THEN DO; PREMNS = B AKRAY (II; B AKR AY < T ) =0 . ; 

2 end; 

1 END; 

ERS = XNBARCAL / 4.9 * EXP U UTS TRT-UTORG ) / 1400 . I ; 

PUT EDIT (*N BAR CAL =• , XNBARCAL, ‘DEL CAL = • , DELCAL) 
(R(NBARCALFORM) ) ; 

PUT EDIT ( 1 EPS IL ON ~ 1 , EPS) ( R { NBARC ALFORM ) ) ; 

CALL NFWPAGE; 

PUT EDIT ( 1 AVERAGE BACKGROUND RASTER 1 , COLID) ( SK IP I 7) , A, 

SK IP t 2 ) , X I 4 ) , A( 80) ) ; 

' DO' 1,-1 TO 48; 

1 ' IL INF =49-1; 
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/* A S&r fi$0-D PROGRAMS // GET BACKGROUND RASTERS 


*/ 


ST«T 

LFVFL 

NFST 


737 

1 

l 

I L l NFOl IT = 1 • ; 

734 

1 

1 

IF ILTNF - 40 | I L I ME = L 1 ULlNE/5>*5 = [LINE THFN IL INEP 
= f L I NG ; 

236 

1 

l 

nn j = i to 39 by 2 ; 

737 

1 

7 

INDEX = 48 * (J-l) +49-1; 

? 39 

1 

? 

CALL CONVERTER ( BAKRAYUNDEX ) , KKAY(J)>; 

7 3^ 

1 

7 

INDEX - 48 *J + I ; 

?40 

1 

2 

CALL CONVERTER (BAKRAY(INDCX) t KKAYCJ+1M; 

?4l 

1 

7 

end; 

?/»*> 

1 ' 

1 

PUT EDIT { I LINEOUTf KKAY, ILtNEOUTI (SKIP, A(2>, X<2), 40 A (7 
X { 7 ) t A{ 7) ) ; 

243 

1 

1 

END; 

744 

1 


PUT EDIT (COLID) (SKIP, X(4), 4(90)); 

745 

1 


LOCATE OUTRECD FILE (OUTFIL) SFT (RGCIO); 

?4A 

l 


FPS - 0.; 

247 

1 


UTSTRT1 ^ 0; 

74S 

1 


UTSTOPI = 0.; 

740 

1 


XNBARCAL = 0.; 

250 

1 


DFLCAL - 0. ; 

?5i 

1 


XN4HCAL = 0* ; 

7S7 

1 


DEL40CAL = 0.; 

753 

1 


XN400NE = 0,; 

? 54 

1 


DEL400NE = 0.5 

755 

1 


BAKRAY = MRAY * 0.1; 

756 

1 


FNnGROUP: UTSTRT = UTSTOP; 

757 

1 


UTSTOP = UTSTRT + TO l F ; 

?5R 

1 


IF IFOF 0 THEN GO TO- WRAPEOF; 

760 

1 


GFTNEXTGROUP: IF UTSTRT >- UTEND THEN GO TO MAX T I MEF XCFF0E9 ; 

267 

1 


IF UTSTOP > UTFNQ THEN UTSTOP = UTEND + *0T; 

764 

1 


CALL NITCOUNTS; 

765 

l 


GO TO TESTRASTEP; 

764 

1 


FNDOFILF: UTSTOP -UTI ME; 

747 

1 


UTFND = UTIMF; 

760 

1 


GO TO PRINTSTATISTICS; 

769 

l 


M AXTT M FFXCEEDFD: PUT EDIT ( * MAX I MUM REQUESTED T I MF OF « f UTEN n , 

1 EXCEEDED 1 ) { SK I P ( 2 ) , A, F(12,3), A); 

770 

1 


GO TO WINDUP; 

?71 

1 


WRADFnp: PUT EDIT ( 1 FND OF FILE FNCOUN T FRED ON INPUT TAPEM (SKIP(2), 

A) ; 

WINDUP: CALI RUNTL ; 

777 

1 


770 

1 


CLOSE FILE ( OUTFIL ) ; 

774 

1 


DISPLAY ('LABEL TAPE OUTFIL AS FOLLOWS ~ 1 ) ; 

275 

1 


PUT STRING (MFSSAGF) EDIT C fl OSO-D BACKGROUND RASTERS FROM * t 
UTBFGNt * TO UTEND, •"•) (A, F(7,7), A f F(7,7)); 

7 76 

1 


DISPLAY (MESSAGE); 

777 

l 


STOP; 

778 

1 


TRFATBACKRWORD: PROCEDURE; 

77 Q 

7 


r F INDEX <= 4 THEN RETURN; 

701 

7 


XN40CAL = XN40CAL + SUMRAY ( I NDEX ) ; 
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/* A S £F OSO-O PROGRAMS // GET BACKGROUND RASTERS 


*/ 


stmt level nest 


78? ? 

?P3 2 

734 7 

785 I 

286 7 

788 2 

2°0 ? 

79 1 ? 

793 7 

?94 ? 

795 7 

796 2 

?9Q ? 

700 ? 

701 1 

30? ? 

303 7 

304 2 

305 ? 

706 7 

307 7 

70S ? 

700 2 

310 7 

311 7* 

312 ? 

713 ? 

314 ? 

715 7 

316 ? 

317 7 

713 ? 

719 1 

7?0 7 


771 2 

32 2 7 

723 1 

374 7 

375 ? 

326 7 

3?7 2 

328 2 

329 7 

330 2 


I40CAL - 140CAL + MRAY( INDEX); 
return; 

END; 

treatsunword: procedure; 

IF -INDEX 0.4 f INDEX > [SIZE THEN RETURN; 

IF 1 RA SRY ( INDEX ) £ »OOOOOOU'B THEN RETURN; 

NBAR - SUBSTR< tRASRYf INDEX) , 9); 

IF NBAR >= 10 THEN GO TO HIGHWORD; 

THISAV = THTSAV + NBAR; 

MTEMP = MTFMP +1; 

RETURN; 

HIGHWORD: IE IF40MIL THEN MHISUNWD4 = MHI SUNWD4 + 1; ELSE MHISUNWD 

=MH!SUNWD + i ; 

RETURN; 

end; 

NTTCOUNTSt PROCEDURE; 

MHORR “ 0; 

MMSTDT = 0; 

MNfllSC = 0; 

MM I S S = 0; 

MSMALL - 0; 

MBAOCT = 0; 

MTOT = o; 

M = 0 ; 

MRAY = 0; 

SUMRAY =0.; • 

MHI WD = 0; 

MIOTOT = 0; 

MHISUNWD = 0; 

M40TOT = 0; 

MHISUNWD4 = 0; 

RETURN? 

END ; 

NEWPAGE: PROCEDURE; 

PUT EDIT { 1 T I ME INTERVAL FROM », MONAMlt I DAT 1 , SS l YEAR 1 , 

I H R T » MINI, SECS 1 , * TO S M0NAM2 , I DAT? , N't I YEAR 7 , 

IHR2, MIN?, SECS2, MS UTSTRT , 1 TOM UTSTOP, •>« ) 

(PAGE, ?A, E ( 3 , 0 ) , A, F(4,0), 2Ff3,0), F(8,4), A, A f E(3,0) 
t A, F { 4 , 0 ) , 2F I 3 , 0 ) , F(8,4), X(16>, 2f A , F ( 8 » 3 )) , A ) ; 

RETURN; 

END; 

GONVERTFR: PROCEDURF (X, KAY); 

DECLARE KAY CHAR ACTER ( 2 ) , X FIXED DEC 1 HAL 13,1), PIC PICTURE # 99V9»; 
OECLARE CTEMP CHARACTER ( 3 ) , IPIC; 

PIC = X; 

CTEMP = PIC; 

IPIC = TRUNC(X) ; 

KAY = SUBSTR (CTEMP, 2, 2)5 
IF IPIC <= 9 THEN RETURN; 


A-10 



/* ASr.F nsn-D PROGRAMS // OFT BACKGROUND RASTERS 


STMT 

33? 

■*33 

**3^ 

340 

341 
34? 
343 
*44 

345 

346 

347 
343 
349 
3 50 
353 
35* 

356 

357 
353 


359 

360 
36* 

368 

369 

370 

371 


*/ 


LEVEL NEST 


2 

7 

7 

2 

? 

1 

? 

7 

2 

7 

7 

2 

7 

7 

7 

2 

? 

2 

7 

2 


1 

I 

? 

1 

! 

t 


2 

7 

7 

2 1 
7 
2 
1 


1 0 r C - IPIC - 9; 

IF IPIC > 76 THEN DO; KAY = '**■; RETURN; END ; 

SURSTR ( KAY , 1 ,1 ) -SUBS TR ( * ABCpET GH I JKL MNDPQR STUVWXYZ * * ■ , I P IC , 1 ) ; 

RETURN; 

END; 

PPOCSUN: PROCEDURE UD, MTOT, MHISUNWO, SIJMTAB, MSUNTAR ) ; 

DECLARE ID CHARACTER*?), SUMT AB ( 600 ) , MSUNTAB ( 600 ) BIT(8); 

DECLARE t, MTFMP FIXED BINARY(R); 

OECLARE CIDENOM, MHISUNWD) FIXED fl t N AR Y ( 3 1 f 0 ) ; 

M10RE JECTED = 0; 

XM400NE = 0*; 

IDFNOM = 0; 

OO 1 = 1 TO MTOT; 

MTFMP = MSUNTABU ) ; 

r E SUMTArMI) / MTFMP > XN40CAL THEN DO; M 1 OR EJ EC TED 
=M10REJFCTFDf 1; GO TO ENDL00P2; END; 

TDENOM = IDENOM + MTEMp; 

XN400NE - XN400NE + SUMTAB(l); 

ENOLOOP2: END; 

PUT EDIT ( MTflT , ID, • MIL RASTERS PROCESSED', M 1 OR EJECTED, ID, 

1 MIL RASTFRS REJECTED 8FCAUSF AVERAGE TO n HIGH*, MHISUNWD, 
ID, 1 MIL RASTER WORDS RFJECTFD BECAUSE HIGHER THAN 10*) 
(SKIP, E{ 10,0) , X ( 1 ) , ? A ) ; 

PUT EDIT (IDENOM, ID, * MIL RASTER WORDS USED*) 

(SKIP, F ( 1 0 , 9 ) , xm , 2 A ) ; 

IF TDENflM <= 0 THEN DO; DEL400NE - 0*; XN400NE = 0.; END; 
FLSE DO; DEL490NF - SORT ( XN400NE ) / I DFNOM ; XN400N n 
- XN400NE / IDFNOM; END; 

RETURN; 

END; 

END; 
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DOS PL/I COMPILER 360N-PL-464 CL3-5 


OSODCORR 


04/24/70 


/* ASCE USO-U PROGRAMS // GET CORRECTED RASTERS */ 


1 

2 

3 

4 


5 

6 „ 


7 

a 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


/* ASCE OSO-D PROGRAMS // GET CORRECTED RASTERS */ 

CRST . . PROCEDURE OPTIONS (MAIN) * « 

DECLARE ILINEPIC PICTURE 1 99 • , 1LINE0UT CHARACTER ( 2 ) 

DEFINED ILINEPIC, COLID CHARACTER! 80 ) , • 

DECLARE RECID POINTER, • 

DECLARE 1 INRECD BASED! REC ID ) ALIGNED, 2(EPS, UTSTRT , UTSTOP, XN0ARCAL, 
DELCAL, XN40CAL, DEL40CAL, XN400NE , DEL4Q0NE, 

BAKRAYU920) -FIXED DECIMAL (3,1)),. 

DECLARE INTAPE FILE RECORD SEQUENTIAL BUFFERED INPUT 

ENVIRONMENT (V(3884) MEOIUM (SYS008, 2400) NOLABEL),. 
DECLARE 01 RASTER EXTERNAL ALIGNED, 2(UTIME, UTIM05) FLOAT BINARY ( 53 ) , 
02 USIZE, rSUBSZ, I8ADCT , ISKPED, IMISS, IENDN, IEOF, 
JFWEEL , JAWEEL) FIXED BINARY (31), 02 ( IRASRY ( 1 920 ) , 
ISC05RU0), l SCO 7R ( 10 ) , ISC27RU0), ISC44R(10)) BIT<16),. 
DECLARE (DATEA, SOATE, YRDATE ) FLOAT B INARY ( 53 } , . 

DECLARE ARRAY { 27 ) , ISUBSCR FIXED B INARY( 31 ) , KAYI40) CHARACTFR ( 2 ) , . 
DECLARE TWOSQRBAK (1920),. 

DECLARE NBAR FLOAT,. 

DECLARE VOLTRAY( 4) PICTURE ' 99 * , VCRAY(4) CHARACTER ( 2 ) 

DEFINED VOLTRAY,. 

DECLARE FNAME CHARAC TER ( 24 ) , ANAME CHARACTER ( 20 L, . 

DECLARE MONAME CHARACT ER ( 4 ) , SGN CHARACTER i 1 ), . 

DECLARE PRT ALLCHAR CHARACTER(l), I FPRT ALL B IT { 1 ) , • 

ON ENDFILE ( INTAPE) GO TO ENDBACKGROUND, . 

COLID = »01 05 10 15 20 25 30 

35 - 40 ■ , . 

CALL NTRAST , . 

UTO = 39490,. 

CALL REWEPH,. 

OPEN FILE (INTAPE),. 

GET EDIT ( UTBEGN , UTEND, PRT ALLCHAR ) (2F(15,5), A(l)),. 

IF PRT ALLCHAR NE 1 # THEN I FPRT ALL = • l ■ B , • ELSE IFPRT ALL- 1 0* B , . 
UTIMF = -100. E40,. 

GO TO FIRSTBACKGROUND,. 

NEXTBACK GROUND . ♦ READ FILE ( INTAPE ) SET (RECID),. 

FIRSTBACKGROUND.. READ FILE I INTAPE ) SET (RECID),. 

IF UTSTOP LT UTBEGN THEN GO TO NEX TB ACKGROUND, . 

IF UTSTOP - UTSTRT LT 0.25 THEN GO TO NEXTBACKGROUND, . 

IF UTSTRT GT UTEND + 0.25 THEN GO TO BACKGROUNDT IMEHI , . 
TWOSQRBAK = 2 . * SQRT (BAKRAY),. 

GO TO TESTRAST I ME , . 

NEXTRASTER.. CALL NXRAST , • 

IF IEOF NE 0 THEN GO TO ENDF i LER ASTERS , . 

TESTRASTIME.. IF UTIME LT UTSTRT - 0.25 THEN GO TO NEXTRASTER,. 

IF UTIME LT UTBEGN THEN GO TO NEXTRASTER,. 

IF UTIME GT UTSTOP THEN GO TO NEXTBACKGROUND,. 

IF UTIME GT UTEND THEN GO TO MAXT I MEXCEEDED , . 

IF JFWEEL = 3 THEN GO TO NEXTRASTER,. 

IF I FPRTALL THEN GO TO GUDRASTER,. 

IGUDCNT = 0,. 

DO I = 5 TO I S I ZE , ♦ 
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IF IRASRY ( I ) AND •OOOOOOOl'B THEN GO TO ENOLQOP * . 

CALL PTWRDC IIRASRY, It N, NBAR),. 

IF NBAR - BAKRAY ( I i LT TWOSQRBAM!) THEN GO TO ENOLQOP,. 

IGUDCNT = IGUOCNT + l t . 

IF IGUDCNT GE 7 THEN GO TO GUDRASTER,. 

ENOLOOP.. END,. 

GO TO NEXTRASTER, . 

GUDRASTER.. SDATE = UTIME + UTO,. 

CALL SDAY (SDATE, IYEAR, YRUATE , MONTH, DAT6A, MONAME),. 

IDATE1 = OATEA, . 

HRS = 24. * (DATEA - IOATE 1 ) , . 

CALL DEGOMS (HRS, ISGN, IHR , MIN, SECS),. 

CALL FLTLUK ( JPWFEL , FNAME , QLAM1, QLAM2, FNX, FFX ) , . 

FN = FNX / EPS,. 

FE -= FEX / EPS,. 

CALL APTLUK (JAWEEL, ANAME > , . 

/* GET ASPECT-EPHEMERIS DATA */ 

CALL GETEPH (UTIME, ARRAY, SPRATE, ASPANG, ISDAY ) , . 

DELT = OSODLP (UTIME, ARRAY),. 

CALL NEWPAGE,. 

SIDTIM - SIDGRN IUTIME + UTO),. 

/* COMPUTE ANO PRINT EPHEMER I S QUANTITIES */ 

CELONG = ATAN ( ARRAY (2) • ARRAY ( I ) ) , • 

TRLONG = (SIDTIM - CELONG) * 57.29578,. 

TRLAT = ARRAY ( 22 ) * 37.29678,. 

CLONGD = CELONG * 57.29578,. 

PUT EDIT ( • EPHEMER IS AT START OF SCAN*, 'CELESTIAL LONGITUDE', 
CLONGD, • DEGREES', 'POSITION VECTOR', ARRAY ID, ' KM') 
(SKIPI2), A, X ( 1 3) , A, F 1 9,2 ) , A, X(15), A, FUO.l), A),. 

PUT EDIT ('TERRESTRIAL LONGITUDE*, TRLONG, ARRAYI2), 

'GEODETIC LATITUDE', TRLAT, ARRAY I 3 ) , 

'GEODETIC ALTITUDE', ARRAY 123)) (SKIP, X I 37 > , A, F(9,2), 

X ( 3 8 ) , F ( 10, 1 ) , SKIP, X ( 4 1 ) , A, F(9,2), X(38), F(10,1),' 
SKIP, X ( 41 ) , A, F { B , 1 ) ) ,. 

RA ■= ATAN (ARRAY ( 20) , ARRAY (19) ) * 57.29678,. 

ROLLA = ARRAY 1 27 ) * 57.29578,. 

DEC = ASIN ( ARRAY (21) ) ♦ 57.29578,. 

ALPHA = ARRAY ( 26 ) * 67.29578,. 

PUT EDIT {'ASPECT AT START OF SCAN', 'PITCH ANGLE' , ALPHA, 

• DEGREES', 'SPIN AXIS', ARRAY ( 19 ) , 'ROLL ANGLE', ROLLA, 
'UNIT VECTOR', ARRAY 1 20 ) , 'RIGHT ASCENSION OF SPIN AXIS', 
RA, ARRAY! 21 ) , ' OECLINAT ION OF SPIN AXIS', DEC) ISKIPI2), A 
, X ( 24 ) , A, F ( 9 , 2 ) , A, X ( 22 ) ♦ A, r<ll,f), SKIP, X(48), A, 

F I 9, 2 ) , X I ?8 ) » A, Fill, 7), SKIP, X ( 30 J , A, F(9,2), XI39), 
Fill, 7), SKIP, X ( 34 ) , A, F(9,2)),. 

SIDTIMO = SIDTIM * 57.29578,. 

PUT EDIT ('SIDEREAL TIME AT START OF SCAN', SIDTIMD, « DEGREES* 
f 'SOLAR', ARRAY I 7‘) , 'UNIT VECTOR', ARRAY ( 8 ) , ARRAY I 9 ) ) 

( SKI P I 2 ) , A, F 1 9 , 3 ) , A, X(54), A, Fill, 7), SKIP, X(95), A, 

F ( 11 , 7) , SKIP, X ( 106 ) , Fill, 7)),. 

PUT EDIT IISKPED, 'RASTER SCANS LOST DUE TO MISSING DATA', 
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ISIZE, * WORDS IN THIS RASTER SCANS I8ADCT, 

'QUESTIONABLE DATA WORDS IN THIS RASTER SCAN', IMiSS, 
'MISSING DATA WORDS IN' THIS RASTER SCAN*) (SKIP, F<5,0), 
X(l) f A) f 

/* PRINT OUT SUBCOMMUTATOR VOLTAGE READINGS */ 

IF ISUBSZ = 0 THEN GO TO NOSUBCOMS,. 

PUT EDIT ( 1 VOLTAGES FROM SAIL ANALOG SUBCOMMUTATOR 1 , ISUBSZ, 

• SETS 1 , * 05 ) 6 VOLT MONITORS *07) 6 VOLT MONITORS 
1 27 > NEGATIVE 6 VOLT SUPPLY* , *44)’ 2*4 KV MONITORS 
( SK l P ( 2 ) , X{ 32 ) , A, X ( 30 ) , F(5,0), A, SKIP, X(13), A, 

X( 12) , A, X ( 4 ) , A, X( 11) , A) 

DO I = l TO ISUBSZ,. 

VCRAY *= ■ 1 , • 

IF NOT I SC05R ( I ) AND *00000001*8 THEN VOLTRAY ( 1) = SUBSTR 
( ISC05RI I J , 9, 8) 

IF NOT I SCO 7R ( I ) AND *00000001*8 THEN VOLTRAY (2) = SUBSTR 
( ISC07RI Il t 9, 8) , i 

IF NOT I SC27R ( I ) AND *00000001*B THEN VOLTRAY ( 3) * SUBSTR 
( I SC2 7R ( I) , 9, 8) 

IF NOT ISC44RIU AND *00000001*8 THEN VOLTRAY (4) = SUBSTR 
( I SC44R ( I), 9, 8)f 

PUT EDIT (VCRAY) ( SK l P, , 4 ( X ( 28 ) » A(2)])r. 

END , « 

CHECKRAS TREND* * IF I ENDN= 1 THEN PUT EDIT! • LAST RASTER BEFORE SUNSET*) 
(SKIP, A ) , • 

IF IENDN = 2 THEN PUT EOIT 

( 1 LAST RASTER BEFORE TIME DISCONTINUITY IN DATA*) (SKIPf?), 

A ) f • 

CALL NEWPAGE,. 

PUT EDIT (COLID) (SKIP(2), X(4), A(80)),. 

DO M = 1 TO 48,. 

ILINE = 49 - M, . 

ILINEOUT - • • 

IF ILINE=l OR I L I NE=48 OR ( I L I NE/5 ) *5= I L I NE THEN ILINEPIC 
= ILINE,. 

KAY = » 

DO J = 1 TO 39 BY ?,. 

ISUBSCR = 48 * IJ-1) + 49 - M,. 

CALL PROCWORD ( KAY.IJ)),. 

I SUBSCft = 48 * J + M,„ 

CALL PROCWORD .( KAY(J+ll),’. 

END,. 

PUT EDIT (ILINEOUT, KAY, ILINEOUT) (SKIP, A<2), X(2I, 40A(2), 

X ( 2 ) , A ( 2 ) ) , . 

END., . 

PUT' EDIT (COLID) <SKIP(2>, X(4), A(80)),. 

GO TO NEXTRASTER,. 

NOSUBCOMS.. PUT EDIT ( *N0 SUBCDH VOLTAGES AVAILABLE*) (SKIP(2), A),. 

GO- TO CHECKRASTREND,.' 

END8ACKGR0UND... PUT EDIT (/END OF FILE ON BACKGROUND RASTER TAPE* I 
(SKIP(2) , A),. 
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GO TO WRAPUP,. 

BACK GROUND! IMEHI ♦ . PUT EDIT ( •MAXIMUM REQUESTED TIME EXCEEDED ON 8ACKGR 
OUND RASTER TAPEM ISKIPI2), A),. 

GU TO WRAPUP,. 

ENDF I LER ASTERS* • PUT EDIT (•END OF FILE ON TELEMETRY TAPEM CSKIP(2), 
A),. 

GO TO WRAPUP,. 

MAXTIMEXCEEDED.. PUT EDIT ( 1 MAX I MUM REQUESTED TIME EXCEEDED ON TELPMETR 
Y TAPE 1 > (SKIPI2J , A) 

WRAPUP.. CALL RUNEPH,. 

CALL RUNTL , . 

CLOSE FILE (INTAPE) 

STOP,. 

ASIN.. PROCEDURE (X),. 

RETURN (ATAN (X / SQRT ( l . - X * X))),. 

END, . 

NEWPAGE.. PROCEDURE,. 

PUT EDIT ('CORRECTED RASTER SCAN STARTING AT UT», MON AME * 

I DAT El , ' , 1 , I YEAR , « UT TIME 1 , l HR , MIN, SECS, 

MM UTIME, MM (PAGE, A, X<2), A, M3, 01, A, F(5,0), A, 
2F( 3, 0 ) , F ( 6 , 2 ) , X(10), A, F{8,3), A),. 

PUT EDIT { FNAME, ' LAMBDA 1 =M QLAMl , 'LAMBDA 2= M QLAM2, 

•KN =M FN, ' KE * • , FE) (SKIP, A, X{6), A, F(6,2), X(4), 

A, F l 6 , 2 ) ♦ X ( 5 ) , A, F ( 8 , 2 ) , X(5), A, F(3,2)),. 

PUT EDIT (ANAME, 'EPSIION *•, EPS, ' N BAR CAL -M XNBARCAL, 
'DELTA =', DELT) (SKIP, A, XllOM A, F(7,2), X(10), A, 

F ( 8 , 2 ) , X( 9) , A, F(8,2) ) 

RETURN,. 

ENO,. 

PROC WORD. . PROCEDURE (OUTSTR),. 

OECLARE OUTSTR CHARACTER ( 2 ) , PIC PICTURE »9M NBAR FLOAT BINARY,. 
DECLARE N, CORRN , QLOG, QLOGTR, QLOGFR , . 

DECLARE ITEMP, IPIC,. 

OUTSTR = * ' , . 

IF ISUBSCR GT ISIZE THEN DO,. OUTSTR = •==',. RETURN,. END,. 

IF IRASRY (ISUBSCR) AND 'OOOOOOOi'B THEN DC,. OUTSTR '++ M . 
RETURN , . END,. 

CALL PTWRDC (IRASRY, ISUBSCR, N , NBAR ) , • 

CORRN = NBAR - B AKRAY ( I SU l3SCR ) , . 

IF NOT I FPRT ALL THEN IF CORRN LT TWDSQRRAK ( ISUBSCR) THEN GD 
TO SETSMALL,. 

IF CORRN LT i.O THEN GO TO SETSMALL,. 

QLOG - LOG? (CORRN),. 

QLOGTR * TRUNC (QLOG),. 

QLOGFR = (QLOG - QLOGTR) * 10.,. 

ITEMP * QLOGFR + 0.5,. 

IPIC = QLOGTR + 1.0,. 

IF ITEMP GT 9 THEN DO,. ITEMP = 0,. IPIC = (PIC + l,. END, . 

PIC = ITEMP,. 

SUBS TR (OUTSTR, 2, l) = PIC,. 

IF IPIC GT 35 THEN DO,. OUTSTR- • M . RETURN,. END,. 
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160 


SUBSTRIQUTSTR, 1»1)=SUBSTR( ' 0123456789A8CDEFGHI JKLHNOPQRSTUVWXY • 



, IPIC, 1),. 

161 


RETURN,. 

162 

SETSMALL 

.. OUTSTR = • — 

163 


RETURN,. 

164 

END* • 


165 

ENO t - 
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/* ASEE OSO-D PROGRAMS // SUM COUNTS IN SU RP AS TFR S--SUR SUM */ 

SUB SUM. . PROCF.QUR F OPTIONS (MAIN),. 

DECLARF ZILCH CHARACTER!!),. 

DECLARE NUMDAY PICTURE , 99.»,. 

DECLARE MONTHNAME CHARACTERS),. 

DECLARE NUMP1C CHARACTER ( 7 ) 

DECLARE NFIRSTNO PICTURF ■ <>9 « , . 

DECLARE UTIM2 FLOAT BINARY (53),. 

DECLARE 1 SUE { 1 000 ) , • 

DECLARE RECIO POINTER 

DECLARE 1 IMRFCO BASFJMR^CID) ALIGNED , 

2 (EPS, UTSTRT, UTSTUP, XN4A«CALt DFLCAL, X.N40C AL , DEL 4 DC AL , 
XN400NF , I1EL400NE, DAK RAY ( 1920} TlXLu DECIMAL (3,1)} , . 

DECLARE INTAPE FILE RECORD SEQUENTIAL BUFFERED INPUT 

ENVIRONMENT (V(38d4) Ft 0 IU M ( SY 5009 , ?40D ) NCLA-^l LEAVE),. 
DEC! ARE l RASTER EXTERNAL ALIGNED, ? (UTIMF, UTI”95) FLOAT BINARY ( c 3 ) , 

2 (IS17E, 1 SUBSZ, 1BADCT, iSKPCD, l MISS, IP NON, I COP , 

JFWEEL, JAWEEL) FIXED QINARY131), 2 { I * A f>P Y ( 1 <>20 ) , ! SCO * * ( 1 0 ) , 

I SC 07R ( 10 ) , I SC 2 7R ( 10) , ISC44<M10)) BITUM,. 

DECLARE (DATE, SDATF, YRDATE) TLOAT H I N A R Y ( 5 3 ) f . 

HFCLARE NBAR FLOAT 

DECLARE COP TYM FLOAT BINARY (63),. 

DECLARE FNAMF CHARACTER ( 74) , AN A ME CHAR ACTCR ( 29 ) , Av*E CM A* ACT ( 4 ) , 
SGN CHARACTER ( 1 ) ,. 

DECLARE RADCO (2,10) 

DECLARE CRDRDR FILF STREAM INPUT 

ENVIRONMENT ( F ( 30 ) MEDItJMf SYSIPT, 1442) ) , . 

ON ENOFILE (CRDRDR) Gr TO FINISH 

ON ENOFILE UNTAPE) GO TO ENPBG , . 

ON TRANSMIT (INTAPF) GO TO BADRG 
N = 1000,. 

LOF V = 11,. 

CALL PLOTS URlJH 1) ,N,LDFV) > . 

CALL PLOT ( 0. 0 P 0,“0. 5E0, LIB ) , . 

ICOUNT ~ 0,. 

INDEC = lb,. 

IRFTA ^ -I 
LINT - 24,. 

I PEN - -3,. 

CALL LGAXS( O.OFOtO.OEO, 'COUNTS/ 140 v$FC* , INDEC , 10.-KT, JO. )M f m.orj f 
0 * 3F0 ) , . 

CALL NTRAS T , . 

UTn = 39490,. 

CALL REWFIL (INTAPF), . 

OPEN FILE ( INTAPE ) , . 

UTIMF = -100. E49,. 

DO IN- = l TO 2 

GET FILF (CRDRDR) EDIT ( L t ( P ADCO ( L , M ) DP M= 1 TO 10)) 

(COLUMN! 1 ) ,F < 2 ) , 10( X(2) ,f < 4, 3))*) 

END,. 

' CALL CRDS 
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CRDS • • PROCEDURE 

GET FILE { CRDRDR ) EDIT ( DA YB , HRB , BM I N f DAYN f HRN , FM I N, l X , I Y ) 
(SKIP, 2 ( X ( 5 ) y F ( 3 ) y X ( 3 J y F(2), X(2), F(2)), 2 ( X ( 5 ) , F ( 2 ) ) ) ,♦ 
UBGN = DAY8 + HRR/24.0 + ( BMl N -l.)/1440.0 
UEND = DAYN + HRN/24.0 + ( EM I N + 1. ) / 1440 . 0 ». 

END f . 

CALL PPINFO 
GO TCI FST8G , . 

NXTBG. • READ FILE (INTAPP) SFT (RECID) 

FS TBG. . READ FILE ( I NT APE) SET (RECID) 

BGTM.. IF UTS TOP LT UBGN THEN GO TO NXTBG , . 

UTOT = UTSTOP -IJTSTR T , . 

IF UTOT LT 0.334 THEN GO TO NXTBG 

IF UTSTRT GT UEND +0.5 THFN GO TO RNGBG 

CALCOM L.340*EXP( ( 300 . O-UTSTRT )/ 1400 . 0 > ) /XNBARCAL y • 

/* PRINT OUT 8G RASTER HEAD */ 

PUT EDIT {'NEW BACKGROUND RASTER FROM UTSTRT, ■ TO ', UTSTOP, 
•NBARCAL = * yXNBARCAL, 1 CALCO* = ',CALCO) ( SK I P ( 3 ) y A , F ( 10 , 5 ) , A t 
F{ 10,5) , SKIP, COLUMN ( 10) f A,F(5y3) ,X(5),A,F(5,3)>,. 

/SPRINT OUT INFO */ 

PRINFO .. PROCEDURE 

PUT EDIT («RING INTEGRATION OF (MX, ,IY, *1 FPQM ', 

DAYB, HRB r ’ HRS BMfN, 1 MINS TO ', DAYN, HRN, 

« HRS FM I N , 1 MINS') i p AGE , A , F ( 3 ) , A , F ( 3 ) , A ,2 F(4),A, 

F(3), A, 2 F ( 4 ) , A, T ( 3 ) , A) 

END , . 

/♦RAW RASTERS */ 

GO TO TESTIM 
NXTR ST * • CALL NXRA ST « . 

IF IEDF NE 0 THEN GO TO NORST ,. 

TESTIM.. IF UTIME LT UBGN THEN GO TO NXTRS T , • 

IF UTIME GT UTSTOP THEN GO TO NXTBG, . 

IF UTIME GT UEND THEN GO TO SWITCH 

IF JFWEEL = 3 THEN GC TO NXTRST , • 

I AWL = JAWEEL 

IF I AWL - 3 THEN IAWL = 2 y . 

IF IAWL * 4 THEN IAWL l 

IF IAWL = 5 OR JFWEEL = 5 THFN GO TO NXTRST , . 

/* PRINTOUT RASTER HEADING HERF */ 

SDATE = UTIME + UTO, . 

CALL SDAY ( SDAT E , I YE AR , YRD ATE , MONTH ,OATE , MONAMF ) , . 

IDATE1-DATE,. 

IYROATE =YRD ATE , . 

HRS ^ 24. ♦ (DATE - I DATED , . 

CALL DEGDMS ( HRS, I SGN, I HR , M I N , SEC $ ) 

CALL FLTLUK ( JFWEEL , FNAME , QL AM 1 , QLAM2 , FNX , FFX ), . 

CALL APTLLJK I JAWEEL t ANAME) 

PUT EDIT ( 1 RASTER SCAN STARTING AT', MONAME , T0ATC1, IYPAP, 

• UT TIME', I HR , MIN, SECS, 'DAY NOTATION', UTIME) 

(SKIP(3) , A , X { 2 ) ,A,F{3yO)yA,F(5fO) ,A,2 F(3,0),F(6,2),X(5),A, 
F(10,5D , « 
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PUT EDIT (FNAME, ANAME) (SKIP, A, X(10),. A) 

IF I S A V NE 1 THFN GO TO CONI , . 

I S AV=0 ,. 

on J = 1 TO 240 ,♦ 

IRASRY(J) = I R A SPY ( J ) OR 'OOOOOniO'K 
END , . 

CONT..IF IFNON = 1 THEN ISAVF - 1 , . 

FLSE ISAV =0 , . 

INOG = IHAllCT + IMISS , . 

IF INOG GE 100 THEN DO, . 

PUT EH IT ( 1 OAD PASr.FP' f UTI<M- ) ( A, MU t 7) ) ,. 
GO TO MXTPST t . 

END, . 

I SUM = 0 ,. 

00 OC = 10 TO 40 1Y 10,. 

00 JR - 1 TO 48,. 

JNDX = 48*JC - JR + l,. 

IF l R ASR Y ( J NOX 1 ANO *00000011*8 THIN GO TO NO IP,, 

CALL PTWPOCI I R ASRY , J NHX , NU M , N8 A [> ) , . 

I SUM = I SUM + NRAR 
NDLP. .END, . 

END , . 

IF I SUM GE 4000 THEN 00,. 

PUT EO I T ( * V AN A R A $T F9‘ * , U T i M F ) ( A , F ( 1 l , 7) ) , . 

END, . 

IF ISIZE LT 500 THEN OH,. 

PUT EDIT I * SHORT R A STEP % UT l MF ) t A , F ( I T, 7 ) 1 , . 

GO TO NXTRST , . 

ENO,. 

IF I SIZE LT 19?0 THEM m J = I SIZE -96 TCI ISI7T 
IRASRY(J) = IRASRYIJ) Ok *00009010*1 ,. 

ENO , . 

IF ISIZE LT 1 9?0 THFN DO, . 

IT = MOD { I X , 2 ) 

IF IT - 1 THEN LOX = 4S*(IX -1H- IY 

ELSE LDX = 48* I X — I Y f L , . 

IF LDX GE ISIZE - 96 THFN Of? 

PUT EDIT ( * I NCOMPL F T C R AST FP *) f Al 
GO TO NXTRST ,. 

END ,. 

END, . 

/* THIS IS THE RING INTEGRATIPN ROUTINE */ 

SUM =0 
SIGMA = 0 
itflg - 10 ,. 

DO N = l TO 10 
LIM = N - 1 ,. 

JXl = IX - LIM 

JX2 = [X + LIM f# 

JYl = IY - LIM 

JY? = IY + LIM , . 
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JDY? = 2* LIM 
DEL = 0, 

NBAD = 0 , . 

RING = 0 t * 

BACK = 0 T . 

IF N = 1 THFN NBOX = 1 , . 

ELSE NBOX = R*LIM 

nn jx = jxi to jx? 

IF JX = JXI OR JX = JX2 THEN JDY = 1 

ELSE JOY = JOY? 

IF JX LE 0 OR JX GT 40 THFN GO TO HA OX 
JT = MOOUX,?) 

00 JY JY1 TO JY2 BY JOY 

IF JY LE 0 OR JY GT 48 THEN GD TO BAOY f . 

IF JT . = i THEN NDX=48*(JX - 1 ) 4- JY 
FLSE NOX = 48* JX - JY +1 
IF NDX GT I S 1 7 F THEN GO TO BAOY , . 

IF IRASRY l NDX ) AND '00000011*8 THFN GO TO BADY 
CALL PTWROC ( I R ASRY , NOX ,NUM , NB AR ) , . 

RING = RING + NBAR 

BACK - BACK + BAKRAY(NOX) , „ 

GO TO NDYLP f . 

BAOY.. NBAO = N3AD+ l , . 

NOVLP.. END. ». 

GO TO NOXLP ♦ . 

BADX.. DO JY 5= JY1 TO JY? BY JDY t# 

NBAO = NBAO ♦ l , . 

END f. * 

NOXLP* • END 

CONTRIB *= RING - BACK , 

IF N GE 3 AND CONTRIB LE SORT t BACK) THEN GO TO ENDI T ,• 
BOX - NBOX 
BAD = NBAD 

IF BOX - BAD LE 0.9 THEN GO TP ENDI T 

GEOM = BOX/tROX -BAD ) fm 

DEL = GEOM * SQRT (RING ♦ BACK) 

SUM = SUM + CONTRIB * GEOM ,• 

SIGMA - SORT! SI GMA*S I GMA + DELUDED , . 

GO TO NDNLP 
FNDI T. . 1TFLG = N-l 

GO TO COMPUT f . 

NDNLP. .END 

COMPUT . . KR = ITFLG *. 

IF KR = 0 THEN KR = 1 f • 

FSIJM = SUM*C ALC Q/R ADCO ( I AWL » KP. ) , . 

ERR = $IGMA*CALCO/RADCOUAWL,KR) 

COR TYM = UTIME + I X*5. 1 2/48/ 1440 t . 

ICOUNT = ICOUNT + l t „ 

IF ICOUNT - 1 THFN DO,. 

UTIMFLAST * CQRTYM,* 

NUTIM = C0RTYM*?4 » . 
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/* ASCF OSH-D PROGRAMS// SUM COUNTS IN SUBRA ST FR S — SUB SUM */ 


ISO 

181 

18 ? 

183 * 

184 

13S 

186 

137 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 
193 
199 
?00 
201 
20 ? 

203 

204 
?05 
?06 
?0 7 
203 

209 

210 
?i i 
212 

213 

214 

215 

216 

217 

218 

219 

220 
? 9 1 
222 

223 

224 

225 

226 

227 

228 

229 

230 


UT I M? - COR T YM t • 

MONTHNAME = MONAME t • 

NU^OAY - I DATE l 9 * 

ZILCH =* 1 * , . 

END,. 

IF F SUM LE 0, 1 E 03 THEN DO,. 

XPAG * tCORTYM * 24 - NUT !M 1 *0*90, . 

CALL PLOT ( XPAG,0-0E0, i IB) l# 

CALL PLOT ( XPAG,0*23E0, 100 ) , • 

GO TQ CDMPTit. 

E-NOt • 

XGPSUM = LOG 10 ( FSUM ) t • 

XGFSUMHIHH ■= L0G1 0( PSMM+CRR ) t • 

XGFSUMLuW = L0G10{ FSIJM— ERR ) * • 

INTEQ = 0 ,. 

IF XGFSUM GT 5 THEN DQ* . 

XGFSU* = XGFSUM— 3 1 • 

XGFSUMHIGH = XGFSUMHIGH-3, . 

XGFSUMLOW - XGFSUMLOW - 3 , . 

INTEO - 6 f 
ENOf » 

YPAG = ( XGF SUM-2 ) * 1 0 /3 f • 

YPAGIOW = ( XGFSUMLOW-?) *10/8, . 

YPAGHIGH = ( XGF SUMH I GH-? ) * 10 / 3 , * 

IF JFWEEL = 2 £ TAWL - 1 THEN INTFQ = INTFO + 1 

IF JFWEEL - 4 £ IAWL = l THEN INTEO = INTEO + 2 f • 

IF JFWEEL = I £ IAWL * 2 THEN INTE'O - INTEO. + 3 , . 

IF JFWFEL - 2 £ IAWL * 2 THEN l NT TO = INT C 0 + 4 

IF JFWEEL - 4 fi IAWL ' * 2 THFN I NT f 9 = TNTCO * 5 , „ 

DIFFER = (CORTYM - UT I MFL AST ) *24 , . 

If DIFFER GF 3 THEN DO,, 

I COUNT =■ 0 , . 

IASTONE.. IOOWN = 

IPEN ^ 3 f , 

AXLEN = ( (UT1MFLAST^24) - l¥UTIMJ*0.90 + 9 , ♦ 

CALL PLOT ( AXLEN t 0# 0£0 t I PCN ) f . 

CALL PLOT! O.OFO, 0*0 EO, IpMWNJ , * 

NUMHRS = UT I MELA ST*24 - NUT I M t • 

NS T ART 1 = NUTIM/24,. 

F I RSTNO * UTIM? - NS T ART 1 f ♦ 

N A N A M = F I RSTNO*? 4 , . 

NF I RSTNO - NANAM , * 

ISYMB = 1 3 , „ 

ICHARNUM = 2,. 

DO J s i TO (NUMHRS + 3) ; 

XPAGHRS = ( J- 1 ) *0. 90 ; 

CALL SYMBOL ( XPAGHRS t i 0* OFO t 0 . 07E0 , I SYMB , 0. 0* 0 1 1 1 00") , * 

END; 

DO J - 1 TO (NUMHRS + 3) t • 

XPAGHRS = (J“l)*0.9O, # 

CALL SYMBOL ( XPAGHRS t 0«0h0,D.07C0t I S YMO , 9 *0f 0 t ID'D*- ), . 
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/* AS£E 

231 
23 2 
231 
234 
?3S 
236 
->31 
238 
2 39 

240 

241 
24? 

243 

244 

245 

246 

247 

248 

249 

250 

251 


252 

253 

254 

255 

256 

257 

258 

? 59 
260 

261 

26 ? 

?63 
?64 
26 5 
266 
?6 7 
268 
?69 

770 

271 

27? 

273 


nso-n PROGRAMS // SU M COUNTS IN SURRASTERS — SUBSUM */ 

CALL SYMBOL ( XPAGHRS ,-0. 14EO»O.IOEO,NFIRSTNO,O.OEO, 1CHARNUM) 
NANAM =: NANAM + l,. 

IF NANAM GE ?4 THEN NANAM = NANAM - 24,. 

NFIRSTNO = NANAM, . 

END,. 

NUMPIC = NUMDAY CAT ZILCH CAT MONTHNAME,. 

NFXGL = 7, . 

CALL SYMBOL (0. OFO ,-0.460,0. 07 FOrNUMP IC, 0 *0E0 , NF XGL ) , . 

IPEN = -3,* 

XP AGE = ( UT l MELAST *24 - NUTIM)*0.90 + 5,. 

IF ICflUNT NE 0 THEN GO TO VERYLAST,* 

CALL PLOT(XPAGEtO*OFO, IPEN) , * 

CALL LGAXSt 0. OEO, 0. OFO, ’COUNTS/ 140 MSECS IMDEC, I 0 *OFO , 90 • OGO , 10.0F1 , 
0.3EO) 

NUT IM - C0RTYM*?4, . 

ENO, . 

LASTONF?.*XPAG = ( C0RTYM*?4 - NUT I M ) *0. 90 , . 

ICOOE = -l f . 

IF YP AG GE 0*0C0 £ YPAG LE lO*OEO THEM 

CALL SYMBOL t XPAG, YPAG, 0.07EO, I NTEQ, 0. OFO , ICDOE ) , . 

1NTEP = 15,. 

XPG = XPAG - 0J017: 

IF’YPAGHIGH GE 0*0F0 £ Y PAGH I GH LE 10. OFO THEN 

CALL SYMBOL ( XPG , YP AGH I GH , 0 * 0 7E0 ? I NT EP , 0 . OEO , I CODE ) , . 

IF YPAGLOW LF 0.0 THEN GO TO TOOLDW,. 

IF YPAGLOW GE 10. OEO THEN GO TO TOOLOW, . 

CALL SYMBOL! XPG, YP AGLOW, 0.07EO, I NTEP,0. OEO, ICOOE), . 

TOOLOW.. IF I SUM GF 4000 THEN DO,. 

NVANALLENR ASTER = 29,. 

GUIDE « YPAGLOW - 0.10,. 

IF GUIDE GE O.OEO THEN 

CALL SYMB0L(XPAG,GUIDE,0.07E0, NVANALL EH RASTER , O.OFO, ICOOF ) t . 
END,. 

UT I M ELAST = CORTYM , • 

/* NOW PRINT THE RESULT */ 

C0MPT1 * . PUT EDIT{ ’INTEGRAL OVER SKR,’ RINGS - • , SUM, ’SIGMA = SIGMA) 
{SKIP , COLUMN {10),A,F(31 , A , F ( 8, U , X ( 3 ) , A , F ( 7, 1 ) ) 

PUT EDIT (’CORRECTED VALUES INTEGRAL = ’ , FSUM, ’ RRAC = SERR) 

(COLUMN! 10), A,F{ 8,1 ),Xt2),A,F(7,l)> ,. 

GO TO NXTRST ,. 

SWITCH. .CALL CRDS ,. 

CALL PR INFO ,. 

GO TO BGTM ,. 

CNDBG. -PUT EDIT {’END OF FILE ON BACKGROUND TAPE’) (PAGE, A) 

GO TO FINISH 

RNGBG. * PUT EDIT ( ’ BG TAPE READS START = U T STAR T C ARDS ENO AT', 
USWCH) { PAGE, A, F { 10, 5) , A, F{10,5)),. 

GO TO FINISH 

BADBG..PUT EDIT (’BAD BG AT’ , UT S TRT ) (SK I P ( 3 ) , A , X( 2 I , F { 10 , 6 1 ) , . 

GO TO NXTBG , • 

NDRS T- - PUT EOT T { 9 ENO OF FILE ON TELEMETRY TAPE’) (PAGE, A) 
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/* AS CE OSO-ll PROGRAMS// SUM COUNTS IN SUBRASTE RS— SUBSUM */ 


274 

275 

276 

277 

278 

279 

280 
281 
23? 
281 

284 

285 

286 
287 
? 88 

289 

290 
29 L 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 
30? 
303 


F IN I SH • .PUT EDIT ! * PROGRAM TERMINATED INTERNALLY*) !SKIP(3),A) 

GO TO LASTONE,. 

VERYLAS T * . DO J = l TO 12,. 

TNTEO = J - 1, . 

YCOOR “ 2*00 + 0.2 5*INTEQ, * 

CALL SYM30L (AXLEN, YCOOR , 0. 07E0 , I NTFO , 0 . 0 EG , I CODE ) , . 

END * . 

NK ^ 2 5 ? • 

SA = AXLEN + 0.3,. 

CALI SYMBOL (SA, 2.00 EG, 0.07E0, * MYL AR— 4AR CM IN 1 ,O.OFO,NK) , 

CALL SYMBOL ( S A , 2 • 2 5F0 » 0 .07 EG , 1 1/2MIL BF-4ARCMIN * f O.OFO t NK) T 

CALL SYMBOL ! S A , 2 . 50EO,O.O7EO, *?MtL BF-4ARCMIN *,O.OEO,NK), 

CALL SYMBOL! $ A, 2. 75F0, 0.07 EG, ‘MYLAR- 1ARCM IN ■ , 0. OFO ,NK ) , 

CALL SYMBOL{SA t 3.O0E0 t 0.07G0,* 1/2MIL BE-1A»C*IN ■,G.OEO,NK), 

CALL SYMBnLC$A,3.25F0,0.07EQ, *2MIL 3E-1ARCMIN *,O.OEO,NK), 

CALL SYMBOL! SA, 3. 5060 ,0.07E0, * MYl AR-4 ARCM l N/ L900 • ,O.OEO.NK ) , 

CALL SYMBOL ! SA , 3 . 75F0 t 0.07E0t * 1/2MIL I3F-4ARC M I N/ l OOO * ,O.OFO,NK ) t 
CALL SYMBOL (SA, 4.0060, 0.07E0, ■ 2MIL B F-4ARCM l N/ 1 OOP *,O.OCO,NK| t 
CALL SYMBOL ($A,4.2560,0.07E0, ‘MYLAR- 1ARCMIN/1000 * ,O.0EO,MK) , 
CALL SYMBOL! SA, 4. 50E 0,0. 07E0, ‘ 1/2 MIL BF-1 ARC M I N/ 1 000 * ,O.OFO,NK) , 
CALL SYMBOL! SA, 4.7560,0. 07E0, 1 ?MIL BE-IARCM IN/10Q0 » ,0.0£0,NK), 

IJKLMN - 24,. 

CALL SYMBOL !SA, 1 .0060,0*07F0, *V MEANS VAN ALLFN RA S TER • , 0 . OE 0 , 
IJKLMN) 

IPEN = 999, . 

CALL PLOTUO. OEO, AXLEN, IPFN) 

PUT FD I T ( * C ALCOMP TAPE MADE * ) I SK I P ( 3 ) , A ) , ♦ 

CALL RUNTL ,. 

CLOSE FILE { INTAKE) 

CALL RUNFIL ( I NT APE ) 

END ,. 
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1 

7 

3 

4 

5 

6 

7 

8 


9 

10 


11 

12 

11 

14 

15 

16 

17 

18 
19 
70 
21 
27 

23 

24 

25 

26 
? 7 
7S 
79 

30 

31 
.32 

33 

34 

35 

36 

37 

38 
*9 


/* A$EF nsn-D PROGRAMS / /SUM COUNTS IN SIJGRASTERS — SUBSUM */ 

subsum*. procfourf options (main),* 

DECLARE I MU p ( 10 00) FIXED BINARY,. 

DECLARE M0NAM1 CHARACTER (4),* 

DECLARE DAYNUM PICTIJRF * 99 * , . 

DEC! ARE CHAR CHARACTER (7),* 

DECLARE ZILCH CHARACTER (1),. 

DECLARE P«=r, to PHI M T r r 

declare i i NR r t.o has rn f recto ) aligned , 

7 (TPS, UT.STRT , UTSTOP 9 XNR ARC AL , DELCAL, XN40CAL, DFL40CAL, 
XN400ME f D EL 4 DON F , BAKPAY(1 970) FIXED DECIMAL (3,1)) 
DECLARE INT \ oq FILE RECORD SEQUENTIAL BUFFERED INPUT 

CNVISONHPNT (V(3HS4) MFDIUMf SYSOO^ , ?400 ) NOLABEL L^AVE),* 
OETLAPC 1 RASTER EXTERNAL ALIGNED, ? (IJTI.MEf UTIM05) FLOAT R l NAR Y ( 53 ) , 

? (ISI7.F, I SUBS 7 , IBADCT, ISKPFD, I^ISS, IENDN, IFUF, 

JEWFPL, JAWFEL) FIXED B I NARY (31), 7 ( I R ASR Y ( 1 9 ^0 1 , tSC05R(I0>» 
I S C 0 7 R ( 10) 9 I SC 7 7R { JO) , I SC44R ( 10 ) ) FUTdM,. 

DFCLAPE (DATE, SO ATE , YROATF) FLOAT B I NARY (5 3) 

DFCLARE NBAR FLOAT , . 

D^Cl ARC FNAMF CHAR ACTFR( 2 a ) , AN AM F CHAR ACT ER ( 20 ) , MONAME CHARACTFR ( 4) , 
SON CHARACTFR( l ) f „ 

DECLARE RADCfl ( °, 10) ,* 

DFCLAPF CRDRDR FILF STREAM INPUT 

ENVIRONMENT ( F ( 80 ) MEDl UM( SYSI PT , 1442 )) , * 

DM ENOFILE (CRDRQP) GO TO FINISH 

ON FNDFILF (INTAPF) 0,0 TO ENDBG, • 

0*' TRANSMIT (INTAPE) GO Tn BA DBG 
N = 1000,. 

LOF V - 11,* 

CALL PLOTS ( IBUF ( 1 ) , N, LDEV) , * 

CALL PLOT (0.0F0,-0*5E0, 11B ) , * 

NCHAR = 17 , # 

CALL SYMRflL (O*75FO,10.OE0,O*15F0, f i/2‘ MIL BC F IL TER • , 0 *0F0 , . 

NCHAR) ,. 

INDFC = 15,* 

I R FN ^ -3 , • 

I3FTA = -l,. 

LINT * 15,* 

CAl L LGAXS (O.OFOtO.OEO, ’COifNTS/140 MSEC 1 , I NDFC , I O. OEO , 90 . DEO , 1 „ OE 7 , 
0.3FO) , « 

IB - 0,* 

CALL NTRAST,. 

UTH := 39490,. 

CALL RFWFIL ( l NT APE ) ,♦ 

OPEN FILF UNTAPF),. 

UT I M E = -100. F40, * 

DD IN = l TO 7 

GET FILF (CRDRDR) EDIT I L , ( P ADCO (L»M) DO M=1 TO 10)) 

(COLUMN (I ),F{?) , 10( X ( 2) , F ( 4«3 ) ) )». 

END, . 

CALL CRDS 
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/* ASGF osn-n PR 1GR4MS//SLJM counts in subpastfrs — SUBSUM */ 


40 

41 

4 9 
4 3 

44 

45 

46 
4 7 
4 R 
4 « 

50 
4 1 
«; ? 

50 

54 


55 

56 


57 

5 4 

59 

^0 

61 

62 

6 “* 
44 
64 
64 
67 
4 5 

69 

70 

71 
7? 
70 

74 

75 

76 

77 


C«OS. . PROCEDURE , - 

GPT PTL c (r«nRO») FO[T ( 0 AYR , HRB * DM I N , QAYN , HP N , FM [ N, I X r l V ) 
(SKIP, 2<Y(5), ClOJf X( 0) , FC7) t X(7), F(7)l, ?( X( 5 ) ,F(?) ) ) ,* 

UBGN = DAYF + HPR/24.0 + ( BM I N -1.1/1440.) 

UFNO = nAYN ¥ HRN/24.0 +(F^TN +1*1/1440.0 ,* 

p ND , • 

CALL PRINPO , * 
on rn fstrg , . 

NXTBO.. 9 FAD F II F { I NT AP F ) SFT (RFC!!)) 

FSTBG.* RFAO FILL UNTAPE) S C T (RECIU) t . 

BGTM*. I P UT STOP l.T UBGN THEN C,n TO NX TUG , . 

U TOT = UTSTOP -UISTRT 

IF UTHT IT 0*334 THEN GO TO NX TOG t. 

IF UTSTRT GT 'JFNO + 0.5 THFN G" TO RN^BG 
CALCGaC 1 *340*EXP( UOO.O-UTSTPT 1/1400. D))/XNB ARC \L 
/* npjMi nUT >JG ° AST FR HFAO */ 

PUT EDIT { 1 N F W BACKGROUND RASTER FROM 1 , UTSTRT , ' TO f f lJTSTn p » 
•NBARCAl = 1 , XNR ARC AL , • CALCO = «»CALCO) { SK IP ( 0 ) , A, f ( 1 0 , 5 ) , A , ' 

F (10, 51 , SKIP, COL U M N( 10l 9 AfF(5f3) v X(5)vA v FC5f^))f 
/SPRINT OUT INFO */ 

PPINFO .. PROCEDURE t . 

PUT F 0 1 T (»PING INTEGRATION OF U,IX, •,* ,1Y, 1 1 FRO* G 
OAVR, HP B t 1 HRS BMIN, 1 MINS TO ** DA YN « H&N, 

» HRS f ,.FMJN, • MINS*) { PAGEt A r F( t A , F{ J J , A F(4),A, 
F { 0 ) f A » 7 F ( 4) , A, F(3) , A) 

FNO 

/*RAW RASTERS */ 

GO rn TFSTtM f . 

NXTRST..CALL NX^AST,. 

IF I EOF NE O THCm GO TO WORST ♦- 
TFSTIM*. IF UTIMF LT IJ3GN THFN GO TO N X T R S T » . 

IF UTIMF GT UT STOP THEN G" TO NXT Bf» , • 

IF liriME GT UEN O THEN GO TO SWITCH 

TF JFWEEL = b TH P N GO TO NX TR $T , . 

I AWL = JAW^fl , . 

IF T AWL - 0 THEN IAWL = 7 

IF I AWL = 4 THEN I AWL -l f . 

IF I AWL = * hp JFWPFL = 5 THFN GO TO NXTrtST , - 

/* PRINTOUT "AST* 7 * HEADING HE RF */ 

SDATF = UTIMF + llTO,* 

CALL SO AY (SOATF, I YEAR , VRO A TF , MONTH f OATF t MPN AM^ ) f . 

T0ATF1=0ATE, . 

I YP D AT P = Y D DATE t * 

HP S = ? 4 . * ( 0 AT F - I Mrcn i. 

CALL DEGOMS (HRS, ISDN, IHR, MIN, srcs ) t . 

CALL FLTI UK ( JFWFEL ,FMA Mp, OLAMl t QLAM7 , FNX t FFX) ♦. 

CALL APTLIIK (JAWF C L, ANAMF ) 

PUT EDIT (»RASTFR SCAN STARTING AT 1 , MONAMF, IDATEl, »,*, !YFA°, 

■ UT T I MF 1 , I HP t MIN, SECS, OIAY NOTATION 1 , IJTIM^) 

( SK IP(O) , A F X ( ? ) ,A,F(3,0) ,A,F (5,0) ,A,7 F(O f O) f F{6 f 7),X(4) t A, 

F ( 1 0 , 5 ) ) 
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/* askf nso-o programs// sin counts in simp astf r s — suhsijm */ 


PUT EO IT (FNAMC.ANAMF) (SKIP, A, X ( 101 . A) 

if isav ne I Th«-n r,n to cont 

I SAV=r» 

nn j = i to pao 

IPASPY(J) = T° ASRY I J J OP, * OOOOOO 1 0 ’ R 
END, . 

CPNT..IF ITNUN = I TIHFN ISAVF = 1 ,. 

ELSE ISAV -n 
IMOr. = IPAOCT + I >MI S S 
IF I NOG OP inf) THFN DO,.- 

, PUT EDIT COAIl PASTER* ,UTIMC)( A,T II l » 7)) 
G n TO NXTRST , . 

F^O,. 

ISOM = 0 

r>C JC = 10 I" 40 pY 10,. 

00 JR = 1 TP AS,. 

JNO* = 4R*JC - JR *• !,. 

IF I RASPY { JNflX ) AND '0000001 l'R THEN GO TO N01P,. 

CALL PTWROC ( I R A S° Y , JNDX , NU M » NH A° ) ,. 

I SUM -I S’JM i- NDA1 

NOl P..r.nn t . 

FNO , . 

IF I S' JM oc /.Con THFN -00, . 

PUT EDIT ( 'VAN A ft ASTFR ' ,1IT I MF ) (A,F(11,7)),. 

FNO , . 

IT ISI7F LT FOO THEN DO,. 

PHI FiJIT ('SHORT PASTFP' ,UTIMF) (A,F(11 T 7)1,. 
on Tn NXTRST,. 

FNO,. 

! c ISI7F LT l'»PO THEN 00 J = 1SI/F -is TO ISIZF ,. 

IRASRY(J) = I ft ASRY ( J ) OR ' 000000 10 ' 0 . 


1 09 

FNO 

t • 

110 

IF ISI/F L T 

19?0 ri^M nn, . 

111 

IT = MHO ( l X t ? i 

11? 

IF IT * 

I THFN LUX = 4S*(IX -U* IY 

11 ’ 


r LSF LOX = 48*IX - I Y ♦ l 

11A 

1 p 1 ox 

OF ISIZF - 14 THEN 00 ,. 

1 1 5 

OUT 

FOir ( • I'lCOMOLFTF p.A$TFP ») 

1 1 6 

on 

TO NXTPST 

n t 


t * 

i \* 

FMO t . 



/* r H T S Is THF Q[>’C y riTFO&ATIQN! pnnTTN^ ! 

11 0 

SUM = 0 

t • 

l?0 

S 1 0 M A - 

0 ,. 

1 ?1 

TTFtn = 

10 , . 

12? 

nn n = 

1 Tfl 10 ,. 

1 

[ = 

f l - l ,. 

124 

IX 1 = 

IX - LI ** ,. 

125 

JX? = 

IX *■ 1 IN ,. 

126 

JY1 = 

17 " LIM ,. 

127 - * 

JV7 = 

IY <• L I M , . 
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/* AS£E PSM-0 P J OGPA*S//SUM COUNTS IN S'1RRASTE*S — SU3SUM */ 


12* 

1? Q 

m 

i** 

134 

13* 

1?-) 

1*7 

13* 

l*o 

140 

141 
14* 
14 * 

144 

145 
14*. 
147 

1 4 i'S 

14^ 

150 

151 
15* 
IS* 
154 
154 
1 5*i 
157 
15* 
150 
140 
14 1 
.14* 

163 

164 

165 
144 
V6 7 
16* 
163 

170 

171 
17* 
17* 

174 

175 
174 
177 
1 7* 


JOY? = ** LIM 
Oh I. “ 0. , . 

N*AO = * t . 

‘ 7 1 MG = * 

3ACK = 0 , • 

IF N - l THEN N3PX = 1 # . 

TLSF N*GY = 3*1.1“ t . 

00 JX = JXi TO JX* . 

IF JX - JXl OR JX * IX? THEN JOY = ] 

ri sf J')Y =: jny? t . 

IF JX LF 0 P.P JX GT 40 TH C N 00 TO BAOX , ♦ 

jt = Mom jx 1 7 ) f # 

on JY = JYL Tn JY* BY JPY 

IF 1 Y LE 0 OR JY GT 4* *T HEN G" yn «AOY 
IT JT - 1 THFN MDX=48*(JX - IH-JY 
FLSF NOX ~ 48# JX - JY f] 

1^ NOX GT ISI/F THFN GO TO RADY , * 

jr TRASPYINOXl AND >0000001 1*3 THEN 00 TO RADY t# 
TALL °TW* , )C( IRASHY t NOX ,NUM,NRAR ) t » 

* l Nil = RT«g + N3AR , • 

RACK = RACK + R AK R AY (NOX ) 1. 

GO TO N**YLP f 
RADY. . MIJAO a NRAIJ4- 1 

N n YL p • » r “0 

GP TO N0X1.P , ♦ 

rapx.. on jy = jy] to jy? *v jpy t „ 

AH AO ^ NR AD + 1 t • 

C ND *. 

00X1 . END 

rnNT*J« - RING - a ACK r . 

IF N GF * AN 3 f.ONTRIR L r SO^T(*ACK) THFN GO TP FN01 T 
OOX ^ N^nx 
mao - N n AO 

IF BOX - RAD LF ).0 THFN GO TO F N P I T 

GFOi v - «nx/(3nx -rad) ,♦ 

OFL = OEOM # SORT{PING + 3ACK ) 

SUM = SUM + r.riMTRrO * GFO« 

SIGMA* = S ) R T t S 1 GM A * S I GM A + PEL *051 ) r . 
r.n to NOMLP ♦ . 

FNOl T.. I TF l G - N-l , . 

G(| TO f. *md|jt 

NDNLP..FNO , # *- 

COMPi JT . * K* - I r C L G , . 

F SUM - SUM*C ALC P/P AO r 0 ( I AW| ? K ,? ) 
n«cj - S IG M A#CAl GP/p AOC0{ TAWLtKP ) 

CUP TY M - UT1MT + IX#5*1?/4R/ 1440 f « 

SA* n I..!F JFMFn ME * THEN GP TO COMPTl*. 

1“ = IR +1 , . 

IF 10 = 1 THEN IP,. 

YON AMI' = -M P N A M c t * 

IDAU - tOATFl 
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/* Asr,r psn-o prugpams//su v counts in subrastfrs — sunsuM */ 


I 79 

lao 

131 
IR7 
1*3 
1 *4 
1 34 
1*6 
1*7 
1*3 
1? ^ 
109 
IOO 

I<J4 
19 4 
194 
197 
1 9 c 
|QO 
709 
?0l 
7Q9 
20* 
704 
709 

206 

707 

7Q9 
?1 9 

711 

712 
71 3 

714 

714 

2 14 
71 7 

713 
719 
223 
7?1 
277 

273 
724 
77 4 
?76 
777 


n AYNUM = T n ATI f • 

7ILCH = • • 

CHAP = MONAMl CAT 7 RCH CAT D4YNIJM,. 

NCHAP. = 7 t . 

CALL SV^Bilt ( O.0Jr0,-9.4F0,n.0 7EO,CHAR,0.0E0,NCHAR ) , . 

I (JR ON = C n R TY^» , 

ENn, . 

IOMEGA = 3, . 

IF I AWL = 7 T MPM IOMEGA - 4,. 

IF fsum I f 9.1 r '03 THEN 00,. 

VfniD - 7.9*(CHRTv\i - PJBGN ) r , 

CALL °L FIT txrOR3,n*oro f 1 1-i) *• CALL PLOT ( XCHR 0 , 0. 70^0 , 1 33 ) , . 
00 Tn Cn^^Tl,* FNO t . 

FSU^O = l.nOlO('-SHM) , . 

IMP =9,; 

if psu*n or *> . 0 thtn on,. 

FSU M n = FSHMH - 3 f . 

I MO = i,. 

CNO, . 

VC90R - (PSLf v n~7}*i 0/3, . 

XC90R = CC^RTYV - TNRGN) *7.0, . 

ichoe = -1 , . 

C A I 1 SYMBOL t XCQnP,YCnOR, O.07E0, IOMEGA , 0 .0 EO , IC90F ) , . 

Jp \ *. n = i TH*N on, . 

VN * YCOOR + 0,7,. 

L t N n E * = 4?, . 

call $ y m q o ( xcn 1 9 , v n 1 0 * 0 7 f 0 , L l N oif x , 0 • 0 r 9 , 1 r no r ) ,. 

IF I GUM or 4030 THFM m T . 


NVN = 70, . 

vcnnt> = ycqop - 0 . 10 ,. 

call SY^f>L(xcnnp,Ycnnn t o.o7F;o,NVN, o.oEO, rcnoc) , . 

FNO,. 

/* N .IK PR I N'T TMF 0 FSULT */ 

C rt, «°Tl..PMT F)f T{ ■ INTEGRAL OVER « , KR , 1 RINGS = # , SUM , » S I G MA =: »,SIG*A) 
' (SKIP , COLUMN! 10 ) , A , F { 3) , A , F ( 8 , I ) , X ( 3 ) ,A,F{7,1 I ) 

PUT rniT (•C^R^FCTFO VALUES 1 NT F ORAL = *,f$UM, f «RAC = ' ,PPP) 

(C n f H'«m< 10) ,A,F{A,1 ) ,X(7),A,F(7,l) ) 

G f i TO NXT»ST 
SWITCH. -CALL C >0S , . 

CALL PR INF(1 , . 

OP TO OGTM 

FN01G..PUT FniT ( 1 FNO np FJLF ON BACKGROUND TAPE*) (PAGE, A) 

GP TO FINISH 

DM GHG# . P'iT 40 1 T ( 1 R G TApf REAns START - UT ST AP T , 1 C AROS r N9 AT«, 

HSWCH) (PAGE, A F ( LO , S ) , A, F(10,S)),. 

GO TO FINISH ,, 

tt AORO. . PU T CHIT ('RAO OG A T f ,'ITSTPT) ( S K I P ( 1 ) , A , X ( 2 ) , F ( 1 0 , 5 ) ) , . 

GH rn NXT3G, . 

NOPST..PUT EOT T ( 1 UNO np FILE ON TELEMETRY TAPE*) (PAGE, A) 

FINISH. .PUT FniT ( ‘PROGRAM TFPMINATFO INTFRNAILYM (SKI<M3)«A) 
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/* as € ^ 'isn-o p»nr,QAv»s//su« counts in surrastf&s — sussu« */ 


77Q 

730 

?^1 
73 ? 
73 s 
?3* 
73 3 
7 yi 

p/1'1 

’/.i 

74 ? 

>4’ 

’44 

?4 '. 
’4fi 
?'t l 
140 
540 
’50 
’51 


AXL C N = 1P.0F0 ,. 

•town = ru'xr.M * i.oron ,. 

C 4 [ 1 AX to (O-O^OtO.nrn, *l)Tl«F IN HAYS' ,-INOrC»TO.orO, 0 . 0”0 
O.S^OO),. 

On JKL = I T'l TO,. 

Minor = JM*1.0cn,. CALI pi ot cfijogc, io.ifo, 1 i»i , . cai l n 
O. fl )rn, 1 o«) , . FMO , . 

5 A = AXLCN + l.o,. T<v«cr,/v - 

CALL SYMTOI.C SA,4 .OF ) , 1.07F0 , IflMFGA, O. oro , ICl’O r 1 , . 

CALL SV.«ont.( SA,4.5Fl,0.07l : 0, l nM FGA,4S.0F0, 1 COOT I , . 

SAS = S* + 0.5,. 

CALL SY«AOL ( S AS , 4 . 0 - 0 , 0 , 07 TO, > 4 ARCM IN ' , O.OFO, \|CHA5 ) , . 

C ALL SYMi'^l ( S A 0,4. 5- 0, 0.07FP, » lAROMN’ , 9. OfO , NCMAA ) , . 

Nl JvfA ■= 4 T , . 

CALL SY“"VIL< SAS.O.OFO.O.O’FO, 
i/\U(?n s'-OI Vf 0 ton T Y 10.10; v — VAN ALL e N •» A 5 TT ’ * ,0 . O^O , NO'-* 
I °FN = o’'),. 

CALL PL"T{ 0 .on, 0 .i)FO, fopM) ,. 

PUT FIMT ( • CALCoMP TAPF MAO c >) ( SKI P ( T ) , A ) , . 

CALL p UNT l 

CLOSF Fit r (INTAPI-I ,. 
tall punriL n nt \ p(-) 
run ,. 


» IG'JN, 
LOTIFIpOF, 
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/* AS£E OSO-D PRHGRAMS//SUM COUNTS IN SUBPASTERS — SUBSUM */ 


1 

2 

3 

4 

5 
4 
7 
rt 


9 

10 


11 

‘12 

1? 

14 

19 

16 

17 

IB 

19 

20 
21 
22 

23 

24 

29 
26 
27 
2 B 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


/* A S £E OSO-D PROGRAMS/ /SUM COUNTS IN SUBRASTERS — SUBSUM */ 

SUBSUM* . PROCEDURE OPTIONS (MAIN),. 

DECLARE IRUr(lOOO) FIXED BINARY,. 

DECLARE MONA Ml CHARACTER (4),. 

DECLARE DAYNU W PICTURE ’99’,. 

DECLARE CHAR CHARACTER (7),. 

DECLARE ZILCH CHARACTER (1),. 

DECLARE RECIU POINTER 

DECLARF 1 INRFCD BASFD(RFCID) ALIGNED , 

2 ( EPS r UTSTRT,' UTSTOP, XNBARCAL, OELCAL, XN40CAL , 0EL40CAL » 
XN400NE, DF.L400NE , BAKRAYI1920) FIXED DECIMAL (3,1)) ,. 

DECLARE INTAPE FILF RECORD SEQUENTIAL BUFFERFO INPUT 

ENVIRONMENT (V(3884) MED I UM ( SYS009 , 2400 ) NOLABEL LEAVE),. 
DECLARE I RASTER EXTERNAL ALIGNED, 2 ( UT I ME, UTIM05) FLOAT BINARY! 53 ) , 

2 (ISIZE, ISIJBSZ, IRADCT, ISKPEn, IMISS, IENDN, 1 EOF , 

JFWEEL, JAWEFL) FIXED B 1 NARY 131), 2 U R ASR Y{ 1920 ) , I$C05P(10), 
r$C07R(10), I SC? 7R ( 10), ISC44PI10)) BITH6),. 

DECLARE (DATE, SO ATE , YROATE ) FLOAT B I NARY ( 53) 

DECLARE N3AR FLOAT 

DECLARE FNAME CHARACTER ( 24 > , ANAMF CH ARA CT ER { 20 ) , MONA MG C HARACTER ( 4 ) , 
SON CHARACTER 1 1 ) 

DECLARE RADCC (2,10) 

DECLARE CRDRDR FILE STREAM INPUT 

ENVIRONMENT (F(80) MFD IUM( SYS IPT , 1442) ), . 

ON ENDFILE (CRDRDR) GO TO FINISH ,. 

ON ENDFILE CINTAPE) GO TO ENOBG, * 

ON TRANSMIT (INTAPE) GO TO RADBG , . 

N = 1000,. 

LDEV = 11,. 

CALL PLOTS ( I BUF ( 1 ) ,N, LDEV) , . 

CALL PLGT<O.OEO,-0.5EO, LIB ) 

NCHAR = 17 

CALL SYMfiPL(0.75F0, 10.0E0,0.15E0, *1/2 MIL BE F I LTER 1 , 0 . OE 0 , 

NCHAR) 

INDEC ^ 1 5 , • 

IPEN = -3,. 

IBETA = -1,. 

LINT = 15,. 

CALL LGAXS ( 0 . OEO , 0 . 0 EO , * C OUN TS/ 140 MSEC 1 , I NOEC , 10 . OEO, 90 . DEO , 1 . OF 2 , 
0.3E0) 

IB - 0,. 

CALL NTRAST,. 

UTO = 39490,. 

CALL RFWF IL ( INTAPE) , . 

OPEN FILF UNTAPfc), . 

UTIME - -100. E40,. 

DO IN = 1 TO 2 

GET FILF (CRDRDR) EDIT ( L , ( R ADCO ( L , M) DO M= 1 TO 10)) 

( COLUMN (1) ,F(2),10(X(?) , F ( 4, 3 ) ) 

END,. 

CALL CRDS 
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/* ASSC OSOH1 PROGRAMS// SUM COUNTS IN SU 3R ASTERS — SUB SUM */ 


40 

41 

4? 

43 

44 

45 

46 

47 
43 

49 

50 
5 L 
5 2 
5 * 

54 


55 

56 


57 

55 

50 

60 

61 

62 

63 

64 

65 

66 
67 
6 H 

69 

70 

71 
7? 

73 

74 

75 

76 

77 


CkDS . . PROCEOURE 

GET FILE (C^URDP) EDIT I D AYH , HR 6, RM I N f DAYN t H»N , EM I N , I X , l Y ) 
(SKIP, 2 I X ( 5 ) , F(3), X(3), F{ 2 ) t X{2), F(?)), 2 I X ( 5 » , F ( 2 ) ) 1 ,♦ 
UBGN = DAYB + HRR/24.0 +{OMIN -l.J/1440.0 ,. 

UEND = DAYN +■ HRN/24.0 +IEMIN +1.)/1440.0 ,♦ 


END 

CALL PR INFO 
GO TO FSTBG 

NXTBG. • READ FILE (INTAKE) SFT (RFCID) 

FST3G.. READ FILE (INTAKES SET (RECID) 

RGTM* - IF UT STOP LT UPGN THEN GO TO NX TBG 
UTOT = UTSTDP -UTSTRT ,* 

IF UTOT LT 0.334 THEN GO TC NXTBG t* 

IF UTSTRT GT UEND + 0.5 THEN GO TO RNGRG i. 

C ALCO= ( 1.340*FXP( (300. O-UTSTR T ) / 1400 . 0 ) } /XNBAPC AL , . 

/* PRINT' OUT BG RASTER HEAD */ 

PUT FOIT ( ' NEW BACKGROUND RASTER FROM * , UTS fPT , 1 TO «,UTSTn;>, 
•NBARCAL = f , XNB ARC AL, 1 CALCO = ',CALCQ) ( SK I P < 3 I , A , F (10, 5 ) , A , 

F ( 10 , 5 ) , SKIP, COLUMN! iO),A,F(5,3),X(5),A,F(5t3>),. 

/SPRINT OUT INFO */ 

PRINFO ♦ . PROCEDURE 

PUT EDIT (’RING INTEGRATION OF (*,IX, 1 , 1 ,IY, M FRO * 

DAYB , HP B , • HRS RMIN, 1 MINS TV DAYN, HRN » 

« HOS EMIN, • MINS’) ( P AGE , A ,F ( 3 ) , A , F ( 3 ) , A ,2 F(4),A» 

F(3), A, 2 F ( 4 ) , A, F ( 33 , A) 

END , . 

/*RAW RASTERS */ 

GO TO TEST l M ,. 

NXTRST • • CALL NXRAST,. 

IF I EOF NE 0 THEN GO TO NDRST , . 

TESTIM.. IF UTIME LT UBGN THEN GO TO NXTRST,. 

IF UTIMF GT UTS TOP THEN GO TO NXT3G , * 

IF UTIMf GT UFNO THEN GO TO SWITCH 

IF JFWEEL = 3 THTN GO TO NXTRST , . 

I AWL = JAWEEL 

IF I AWL = 3 THEN IAWL = 2 

IF I AWL = 4 THEN IAWL = L ,. 

IF I A wL = 5 OR JFWEEL ^ 5 THEN GO TO NXTRST 

/* PRINTOUT RASTER HEADING HERE */ 

SDATE = UTIME + UTO, . . 

CALL SDAY ( SCAT E , I TEAR , YRO ATE , MONTH , DAT E , MON ArIC ), . 

I DAT E 1=DATE , • 

IYRDATE = YRD ATE » • 

HRS = 24. * (DATE - I DATE l ) , • 

CALL DEGDMS ( HR S , I SON , I HR ,M I N, SECS ) 

CALL FLTLUK ( JEW F El. , FN A ME , OL AM 1 , QL AM2 , FNX , F T X ) , . 

CALL APTLUK (JAWCFL, ANAMF > 

PUT EDIT ('RASTER SCAN STARTING AT', MONAME , I DATL } • I YEAR , 

• UT TIME', I HR , MIN, SECS, 'DAY NOTATION', UTIMF) 

(SKIP (3) , A , X ( 2 ) ,A«F(3,0) ,A,F(5,0) ,A,2 F ( 1 , 0 ) , F ( 6 , 2 ) , X ( 6 ) , A , 

F ( 10 , 5 ) ) 
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78 

79 

80 
81 
82 

83 

84 
35 
86 

87 

88 

89 

90 

91 
9? 

93 

94 

95 

96 

97 

98 

9Q 

100 

101 

102 

.103 

104 

105 

106 

107 

108 
109 

no 

in 

li? 

in 

114 
1 15 
116 
117 
lie 

119 

120 
121 
122 

123 

124 

125 
i 26 
127 


PUT EDIT ( FNAME , ANAME) (SKIP, A f X(10), A) 

IF I SAV NE l THEN GO TO CONT ,. 

1 S A V-0 , • 

DO J = 1 TO 240 

IRASRY(J) =1RASRY(J) OR ‘OOOOOOIO'B 
END 9 . 

CONT... IF IFNDN 1 THEN I SAVE = 1 

ELSE 1 SAV “0 

INOG = IBADCT + IMISS ,. 

IF INOG GE 100 THEN DO,. 

PUT EDIT { * B AD RASTER • t UT IME ) ( A, F ( 1 1 f 7)) f . 
GO TO NXTRST ,. 

END, « 

I SUM - 0 

DO JC *= 10 TO 40 BY 10,. 

00 JR = 1 TO 48,. 

JNDX = 48*JC - JR + U, 

IF I R ASR Y ( JNDX ) AND 'OOOOOOll'B THEN GO TO NDLP, . 

CALL PTWRDC U R ASR Y, JNDX, NUM,N8AR) 

ISUM = I SUM + NBAR ,. 

NDLP*. END,. 

END,. 

IF I SUM GE 4000 THEN DO,. 

PUT ED I T ( • VAN A RASTER * , UT I ME ) (A,F(U,7)),. 
END,. 

IF ISI2E LT 500 THEN DO,. . . 

PUT EDIT (•SHORT RA STE R • , UT I ME ) ( A,F{ i l, f I ) ,. 

GO TO NXTRST,. 

END,. 

IF ISIZE LT 1920 THEN DO J = ISIZE -96 TO IS1ZE 
IPASRY(J) = IRASRY(J) OR '00000010*3 
END 

IF ISIZE LT 1920 THEN DO,. 

IT = MOD ( I X , ? ) 

IF IT = l THEN 1DX = 48MIX -1) + IY 

ELSE LDX = 48 * I X - I Y + i f . 

IF LDX GE [SIZE - 96 THEN DO 

PUT EDIT ('INCOMPLETE RASTER ') (A) 

GO TO NXTRST' , . 

, END ,. 

END,. 

/* THIS IS THF RING INTEGRATION ROUTINE */ 

SUM = 0 ,. 

SIGMA = 0 
ITFLG = 10 f . 

DO N * 1 TO in t, 

LIM = N - 1 t . 

JXl - IX - LIM 

JX2 = IX + LIM ,. 

JY1 = IY - LIM f . 

JY2 = IY + LIM 
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/* A $8 F (1Sn-0 PROGRAMS/ / SU M COUNTS IN SUBRASTFRS — SUB SUM */ 


1 ?8 
12^ 
1 30 
131 
13? 

133 

134 
13 5 

136 

137 
1 38 

139 

140 

141 
147 
147 

144 

145 

146 

147 

148 

149 
1 5() 

151 

152 

153 

154 

155 

156 

157 

158 
1 59 
160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 
177 

173 

174 

175 

176 

177 

178 


JDY 2 - 7* HM 
DEL = 0, 

NB AD = 0 f . 

P I NG • = 0 , . 

BACK = 0 t . 

IF N = 1 THEN NBHX = 1 , . 

ELSE NBOX = B*LIM 
DO JX = J X 1 TO JX2 f . 

IF JX « JX1 OR JX = JX2 THEN JOY = L 

ELSE JDY = JDY? ,. 

IF JX LE 9 OR JX GT 40 THEN GO TO BADX 
JT = MOD ( JX 1 2 ) , . 

DO JY = JY1 TC JY2 BY JDY , . 

IF JY LE 0 OR JY GT 48 THEN GU TO BADY 
IF JT = l THEN NHX=48* ( JX - 11 f JY 
ELSE NOX = 48*JX - JY +1 
IF NDX GT I SIZE THEN GO TO RADY , . 

IF IRASRY C NDX ) AND ‘0000001PB THEN Gu TO RADY 
CALL PTWROC ( I RASRY , NDX f NUM , NB AR ) , . 

RING - RING + NR AR , . 

BACK “ BACK +■ BAKRAY(NDX) 

GO TO NOYLP , . 

BADY * . NBAl) - NB A0+ 1 

NDYLP . * END 

GO TO NO XL P , * 

BADX . * DO JY = JYl TO JY2 BY JPY 
NB AO = NO AD + 1 
END , . 

NDXLP-- END 

CQNTRIB = RING - 3ACK 

IF N GE 3 AND CONTRIB LC S GRT ( BACK ) THEN 90 TU LNDIT 
BOX = NHOX , . 

BAD = NB AH T . 

IF ROX - BAD LE 0.9 THEN GO TO INOIT ,♦ 

GEOM - BOX / ( BOX -BAD) 

DEL = GEOM * SQRT ( R I NO +■ BACK) 

SUM = SUM + CONTRIB * GEOM , . 

SIGMA - SQRT< SIGMA*SIGMA + DFL^DEL) , . 

GO TO N,9NL P 
ENOIT.. ITFLG = N-l 

GO TO COMPUT ,♦ 

NDNL D . .END , . 

COMPUT. .KR = ITFLG 

IF KR - 0 THFN KR = 1 ,. 

F SUM - SUM^CALCO/PADCPII AWL,KR) 

FRR - SIGMA*CALCO/RAOCQ( IAWL,KR) ,. 

CORTYM = UTIME + 1X^5.17/48/1440,. 

SAMP I ♦ . 1 F JFWEEL NE ? AND JFWETL NE i THEN GO TO COMPTW' 

IF I AWL NE 1 THEN GO TC COMPTl,. 

IB = IB +1,. 

IF IB = 1 THEN DO,. 
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/* ASEE 

179 

180 
18 L 
L 82 
183 

L 84 

185 

186 

187 

188 

189 

190 

191 

192 

193 
195 

197 

198 

199 

200 
201 
202 
203 
'204 

205 

206 
207 

— 208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

213 

219 

220 
221 
222 

223 

224 
2 25 

226 

227 


USO-D PROGRAMS// SUM COUNTS IN SUBRASTERS — SUBSUM */ 

M0NAM1 = MONAME#. 

1DATI - I OAT El , • 

DAYNUM = I DAT 1 ? . 

ZILCH = 1 * 

CHAR = M0NAM1 CAT ZILCH CAT DAYNUM,. 

NCHAR = 7,. 

CALL SYMBOL! 0. 0 EO ,-0 . 4E0, 0 .07E0 , CHAR ,0.0E0, NCHAR) 

TUBGN = CORTYM, • 

END,. 

IOMEGA - 3 T *. 

IF I AWL = 2 THEN IOMEGA - 4,. 

IF JFWEEL 1 THEN IOMEGA=O t . 

IF F SUM LE 0.1E03 THEN DO,. 

XCORO = 5*t CORTYM - IURGNJ/6; 

CALL PLOT ( XCORD , 0. OEO , 118) CALL PLOT { XCORD , 0 . 20 EO , lOR . 

GO TU C0MPT1 , « * ENO,. 

FSUMO •= LOG 10 ( F SUM ) , « 

IMP = 0, . 

IF F SUMO GT 5.0 THEN DO, • 

FSUMO - FSIJMP - 3,. 

IMP ' It. 

END,. 

YCOOR = ( F S U MO— 2)*10/3,. 

XCOOR - 5* ( CORTYM - IUBGNI/6; 

ICODE = -l,. 

CALL SYMBOL! XCOOR fYCnnR,0.07E0, IOMEGA, O.OEOt ICOOt > ,. 

IF IMP = l THEN D0». 

VN ^ YCOOR + 0.2,. 

L1NDEX = 62 » . 

CALL SYMBOL (XCOOR, VN,O.O7E0 »L INDEX, O.OEO, ICODE) 

END,. 

IF I SUM GE^OOO THEN DO,*. 

NVN - 29,. 

YCOOR = YCOOR - 0.10,. 

CALL SYMBOL (XCOOR, YCOOR, 0.07E0, NVN, O.OFO, ICODE) 

END,. 

/* NOW PRINT' THE RESULT */ 

COMPT 1 . ♦ PUT EDIT! * INTEGRAL OVER 1 , KR , 1 RINGS - SUM, •SIGMA = ', SIGMA) 
{SKIP , COLUMN ( 10 ) , A , F ( 3) ,A f F!8,l),X(3),A f F{7,n) 

PUT EDIT ( f CORRECTED VALUES INTEGRAL = 1 ,FSUM, 1 BRAC = * , ERR ) 

(COLUMN! 10) ,A,F(8,1),X(2) , A, F ( 7, 1) ) 

GC TP NXTRST ,. 

SWITCH. .CALL CRDS ,. 

CALL PRINFf) 

GO TO BGTM ,. 

END8G..PUT EDIT ('END OF FILF ON BACKGROUND TAPE') (PAGE, A) ,. 

GO TO FINISH ,. 

RNGBG..PUT EDIT ( • BG TAPE READS START = *, UTSTAP T ,• CARDS END AT * f t 
* * USWCH) ( PAGE, J A , F { 1 0, 5 ) , A, F(10,5>>,. 

GO TO FINISH - * ■ 

BAOBG..PUT EDIT ('BAD BG AT' , UTSTRT ) { SK I P ! 3 ) , A , X ( 2 ) , F { 10 , 5 ) ) , . 
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228 

??9 

230 

231 
23? 

233 

234 

235 

236 
?37 

238 

239 

240 

241 
24? 

243 

244 

245 

246 
? 47 

248 

249 

250 


GO TO NXTBGt* 

NOP S T • . PUT FO l T ( 1 END OF FILE ON TELEMETRY TAPE 1 ) (PAGE, A) 

FINISH. .PUT EDIT ( 9 PROGRAM TERMINATED INTERNALLY*) CSKIP(3),A) t* 
AXLEN = 25.0FC; 

BONN = IUBGN * l.OEOO , . 

CALL AX I S { O.OEO, O.OEO, 1 * ,-l IB , AXLEN, O.OEO , 3GNN , 1 .2E0 ) ; 

DU JKL = 1 TO 30; 

FUDGE = JKL * 20 EO ; 

CALI PLOT (FUDGE, 1 0. OHO, 1 IB) ; 

CALL PLOT (FUOGF, 9.80E0, LOB); 

END; 

DO JKL = 1 TO 30; 

FUDGE = JKL*0.20E0;- 

CALL PLOT (FUDGE, O.OEO,! IB); 

CALL PL0T(FUDGF,-0.2E0, 103) ; 

END; 

IPEN = 999,. 

CALL PLOT(O.OEO, O.OEO, IPEN),. 

PUT EDIT ( f CALCOMP TAPE MADE 1 ) ( SK I p ( 3 ) ,A),. 

CALL RUNTL 

CLOSE F ILE ( INTAPE) 

CALL RUNFIL (INTAPF) 

END 
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/* ASCE OSO-a PROGRAMS// PLAGE INTEGRALS — PLGINT */ 


1 

2 

3 

4 

5 


6 

7 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


33 

34 

35 

36 

37 

38 


/* A S£E QSO-D PROGRAMS// PLAGE INTEGRALS — PLGINT */ 

PLGINT.. PROCEDURE OPTIONS (MAIN) ,. 

DECLARE CRDRDR FILE STREAM INPUT 

ENVIRONMENT (F(80) MEDI UM< SYS I PT , 1442 ) ) , . , 

DECLARE ICLC lOJtlROl 10) 

DECLARE RECIO POINTER f . 

DECLARE l INRECD BAS EDI REC ID I ALIGNED , 

2 (EPS, UTSTRT, UTSTOP, XNBARCAL , DELCAL, XN40CAL , DEL40CAL t 
XN400NE , DEL400NE, BAKRAYC1920) FIXED DECIMAL (3,i)> ,. 

DECLARE INTAPE FILE RECORD SEQUENTIAL BUFFERED INPUT 

ENVIRONMENT (V(3884) MED IUM( SYS009 , 2400 ) NOLAQFL LEAVE),. 
DECLARE 1 RASTER EXTERNAL ALIGNED, 2 (UTIME, UTIM05) FLOAT BINARY ( 53 ) , 

2 (ISIZE, ISUBSZ, IBADCTt ISKPED, IMISS, IENDN, 1E0F, 

JFWEEL, JAWEEL) FIXED B1NARY(3I>, 2 ( I RASRY { 1920 ) , ISC05RUO), 
I SC07R ( 10 ) t I SC2 7R (10), ISC44R(10U BITC16),. 

DECLARE (DATE, SD ATE, YRDATE) FLOAT B INARY ( 53) ,. 

DECLARE NBAR‘ FLOAT 

DECLARE FNAME CHARACTER ( 24) , ANAME CHARACTER ( 20 ) , MONAME CHAR ACTER ( 4 ) » 
SGN CHARACTER 1 1 ) 

DECLARE RADCO (2,10) 

ON ENDFILE (INTAPE) GO TO ENDBG , . 

ON ENDFILE (CRDRDR) GO TO ENDCD , . 

ON TRANSMIT (INTAPE) GO TO BADBG, ♦ 

CALL NTRAST , • 

UTO - 39490,. 

CALL REWFIL ( INTAPE) ,. 

OPEN FILE (INTAPE),. 

UTIME = -100. E40,. 

DO IN = 1 TO 2 

GET FILE (CRDRDR) EDIT ( L, { RADCOI L, M ) DO M = i TO 10)) 

( COLUMN (1),F(2) ,10( X(2 ) ,F(4, 3 ) ) ) , • 

END, . 

CALL CRDS 
CALL PRINPUT,. 

GO TO FSTBG, • 

NXTBG. .READ FILE (INTAPE) SET (RECID) 

FSTBG. .READ FILE I INTAPE) SET (RECID) 

BGTM • • IF UTSTOP LT UBGN THEN GO TO NXTBG, • 

UTQT * UTSTOP - UTSTRT 

IF UTOT LT 0.334 THEN GO TO NXTBG,. 

CALCO=f 1 .340*EXP I ( 300. 0— UTSTRT ) /I 400.0 ) ) /XNBARCAL , • 

/* PRINT OUT BG RASTER HEAD */ 

PUT EDIT ( • NEW BACKGROUND RASTER FROM f , UTSTRT , 1 TO •, UTSTOP, 
•NBARCAL = •, XNBARCAL, 1 CALCO = ^CALCO) ( SK l P < 3 ) , A , F ( 10 , 5 > , A , 

F ( 10,5) , SKIP, COLUMN! lO) t A,F(5,3),X(5)tA,F{5,3)),. 

GO TO TESTIM 
NXTRST. .CALL NXRAST 

IF IEOF NE 0 THEN GO TO NDRST ,. 

TESTIM.. IF UTIME LT UBGN THEN GO TO NXTRST 
I AWL " JAWEEL 

IF I AWL = 3 THEN IAWL =2 ,. 
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39 

IF I AWL = 4 THEN I AWL =1 


40 

IF l AWL = 5 OR JFWEEL =5 T|H£N QO 

TO NXTRSTf. 

41 

CALL PUTRAS 


42 

DO JPL = 1 TO NPL 


43 

. IX = I CL t JPL ) 


44 

IY [ROIJPLJ 


45 

CALL RINGINT 


46 

END f • 


47 

CALL CRDS ,. 


48 

CALL PR INPUT ,. 


49 

GO TO BGTM 


50 

ENDBG. .PUT EOIT ( • END OF FILE ON BACKGROUND TAPE') {PAGE, A) 

51 

GO TO FINISH 


52 

NDRST..PUT EOI T C ' END OF FILE ON TELEMETRY TAPE') (PAGE, A) 

53 

GO TO FINISH ». 


54 

BADBG..PUT EDIT (•FOUL UP ON BG TAPE* ) 

{ SK IP ( 2 ) t A ) 

55 

GO TO 8GTM , . 

* 

56 

ENDCO • .PUT EDIT ( • END OF FILE ON CARO 

READER * ) (PAGEfA) ,, 

57 

GO TO FINISH i * 


58 

FINISH. .PUT FOIT ( 1 PROGRAM TERMINATED 

INTERNALLY • ) (SKIP(3),M 

59 

CALL RUNTL *. 


60 

CLOSE FILE ( INTAPE) t . 


61 

CALL RUNFIL (INTAPE) 


62 

RINGINT.. PROCEDURE t . 



/* THIS IS THE RING INTEGRATION ROUTINE */ 

63 

SUM - 0 t. 


64 

SIGMA - 0 f . 


65 

ITFLG - 10 9 . 


66 

DO N - 1 TO 10 ,♦ 


67 

LIM = N - 1 ♦ . 


68 

JX1 = IX - LIM 


69 

JX2 = IX + LIM f . 


70 

JYl = IY - LIM 


71 

JY2 = IY + LI M 


72 

JDY2 = 2* LIM f . 


73 

DEL = 0. v . 


74 

NBAD = 0 


75 

RING " 0 f . 


76 

BACK -Ot. 


77 

IF N = I THEN NBO X = 1 


78 

ELSE NBOX = 8*L I M , . 


79 

DO JX = JX1 TO JX2 ♦. 


80 

IF JX = JXI OR JX = JX2 THEN 

JDY = 1 

81 

ELSE 

JOY = JOY? 

82 

IF JX LE 0 OR JX GT 40 THEN 

GO m BADX 

83 

JT = MOD ( JX t 2 ) ». 


84 

DO JY = JYl TO JY2 BY JOY 


85 

IF JY LE 0 OR JY GT 48 THEfN 

1 GO TO RADY 

86 

IF JT = 1 THEN NDX=48*(JX - D+JY 

87 

ELSF NOX = 48*JX 

- JY +1 

88 

IF NDX GT IS17E THEN GO TO 

RADY 
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89 

90 

91 

92 

93 

94 

95 

96 

97 
90 
99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 


IF IRASRY (NDX) AND "00000011" B THEN GO TO BADY 
CALL PTWROC { IRASRYt NDX* NUM , NBAR ) 

RING = RING + NBAR 

BACK = BACK + BAKRAY I NDX ) , . 

GO TO NOYLP , . 

BADY. . NBAD = NBAD* 1 

NDYLP. • END »• 

GO TO NDXLP , . 

B ADX. * DO JY = JY1 TO JY2 BY JDY 

NBAD = NBAD + 1 , • 

END 9 . 

NDXLP* • END 

CONTRIB - RING - BACK 

IF N GE 3 AND CONTRIB L6 SCRT ( BACK ) THEN GO TO ENDiT 
BOX = NBOX 
BAD NBAD ,. 

IF BOX - BAD LE 0.9 THEN GO TO ENDIT 
GEGM - BOX/CBOX -BAD) ,* 

DEL = GEOM * SQRT (RING + BACK) ,. 

SUM = SUM + CONTRIB * GEOM , . 

SIGMA = SQRT( SIGMA*SI GMA + DELUDED *. 

GO TO NDNLP 

ENDIT.* ITFLG = N-l ,* 

GO TO COMPUT t * 

NDNLP. .END 
COMPUT*. KR - ITFLG 

FSUM = SUM*C ALCO/RADCO ( IAWLtKR) , . 

ERR = SIGMA*CALCO/RADCO( IAWL,KR) 

/* NOW PRINT THE RESULT */ 

PUT EDIT ("RING INTEGRATION OF REGION AT (SIX, , IY, ")") 

( SKIP (3) , A, F ( 3 ) , A* F ( 3 ) ,A) ,* 

PUT EDIT ("INTEGRAL OVER ", KR, " RINGS = ", SUM, "SIGMA = ", SIGMA) 
(SKIP , COLUMN I 10) ,AtF(3),A,F{8,l),X(3),A,F( 7, 1) ) , . 

PUT EDIT ( "CORRECTED VALUES INTEGRAL = " , FSUM, 1 BRAC = " , ERR ) 

( COLUMN ( 10),A,F(8,1),X(2),A,F(7,1>) 

END RINGINT 
CRDS. .PROCEDURE ,. 

GET FILE ( CRDRDR ) EDIT ( DAY , HR , XMN, N f ( I CL (M ) , I RO ( M ) DO 1 TO N)) 
(SKIP,F(3),2{X(l),F(2)),F(2),9(X(2),F(2)tX(l),F(2))) 

UBGN = DAY ♦ HR/24. 00+( XMN-1 . 0 ) / 1440. 0 , . 

NPL = N 
END CRDS 

PUTRAS.. PROCEDURE,. 

/* PRINTOUT RASTER HEADING HERE */ 

SPATE = UTIME + UTO,.. 

CALL $DAY (SPATE, l YEAR, YRDATE, MONTH, DATE, MONAME) , . 

I DATE1=DATE , • 

I YRDATE * YRDATE , . 

HRS ^ 24* * (DATE - IDATE1 ) , • 

CALL DEGDMS (HRS , I SGN, I HR , M I N , SECS ) 

CALL FLTLUK ( JFWEEL, FNAME, QLAM1, 0LAM2 ,FNX ,FEX ) , . 
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CALL APTLUK (JAWEEL, ANAME) 

PUT EOIT ( * RASTER SCAN STARTING AT 1 , MONAMF , I DA Tfc 1 , 1 , * t IYCAk 
« UT TIME', 1 HR, MIN, SECS, 'DAY NUTATION', UUMLJ 
(SKIP ( 3 ) , A , X I 2 ) ,A,F{3 f 01tA,F(5f0) tA,2 F(3,0),M6.2)»X(5J»A, 
F I 1 0, 5 1 ) 

PUT EDIT (FNAME, ANAME) (SKIP, A f X(10), A) 

END PUTRAS , . 

PRINPUT** PROCEDURE , ♦ 

PUT EDIT ( ' INTEGRAT ICN OF ',NPL, f PLAGES ON DAY * , DAY , HR , 

* HRS 1 , XMN j 1 MINS UBGN “ 1 , U6GN ) ( PAGE , A , F ( 2 ) , A , F ( T ) , X ( ? ] , 
F(3‘),A,XI21 ,F(3),A,F(I0t5) ),* 

1 AO * PUT EDIT ( 1 ICL« IRQ' ) ( SK I P ( 3 ) , COLUMN ( 9 } , A, COLUMN ( 1 9 ) , A 1 , • 

141 DO M -= 1 TO NPL , • 

14 ? PUT EDIT (ICUMJ , IRO(M) ) ( SK I P , COLUMN ( O ) , F I ? I t COLUMN ( 1 9 ) 

F(3)) t . 

ENC,. 

END PRINPUT,. 

END 


134 

135 


136 

137 

138 

139 


143 

144 

145 
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1 

? 

3 

4 

5 

6 


7 

R 

9 

10 

11 

12 


13 

14 

15 

16 

17 

18 

19 

20 
21 
2? 

23 

24 

25 

26 
27 
?R 

29 

30 

31 

32 

33 

34 

35 


/* A SEE OSO-D PROGRAMS ft GET AVERAGE CORRECTED RASTERS */ 

ARST . . PROCEDURE OPTIONS (MAIN),. 

DECLARE* ILINEPIC PICTURE • 99 * , ILINEOUT CHARACTER ( ? ) 

DEFINED ILINEPIC, COLID CHARACTER f 80 ), . 

DECLARE RECID POINTER, « 

DCCLARF 1 INRECD BASFOtRECID) ALIGNED, 2 (EPS, UTSTRT , UTSTOP, XNBARCAL, 
DELCAL, XN40CAL, DEL40CAL, XN400NE , DEL^OONE, 

0 AKR AY( 1920 ) FIXED DECIMAL (3,1)),. 

DECLARE REC I DO POINTER,. 

DECLARE 01 OUTRECD BASED (RECIDO) ALIGNED, 2(UTAVL0, UTAVHI, ITOTCTAV, 
JAWEELAV, JFWEELAV, I NCOHPLETECT AV, LARGM I SSINGCTAV, 

L ARGQUESTCTAV, I DECL AR EDBAOCTAV, NRASTAV, 

I FR STI NDAYAV , AVR AY ( 1920 ) ) , . 

DECLARE' OUTAPF FILE RECORD SEQUENTIAL OUTPUT ENVIRONMENT (V(7732) 

M F 0 1 UM ( $ Y S 0 0 8‘, 2400) NOLABEL LEAVE),. 

DECLARE MESSAGE CHAR ACTER ( 60 ) , INDICAV BITC1),. 

DECLARE DUMMY CHAR AC TER ( 1) , . 

DECLARE SYSIN FILE STREAM INPUT ENVIRONMENT {Ft 801 
MFD I UM ( S YS l PT » 1442)),. 

DECLARE INTAPE FILE RECORD SEQUENTIAL BUFFERED INPUT 

ENVIRONMENT (VC3884) MEDIUM ( SYS009 , 2400 ) NOLABEL LEAVF ) , • 

DECLARE 01 RASTER EXTERNAL ALIGNED, 2(UTIME, UTIM05) FLOAT BINARY ( 53 ) , 
02 (ISIZE, ISUBSZ, I8ADCT , ISKPED, IMISS, IENDN, IEOF, 
JFWEEL, JAWEEL) FIXED BINARY (311, 02 ( [RASRY 1 1920 ) , 
TSC05RU0), I SC07R ( 10 ) , ISC27RU0), ISC44RU0)) BIT(16),. 

DFCLARF UTO FLOAT B INARY ( 53) , • 

OECLARE NRAR FLOAT,. 

DECLARE KEYUSED PICTURE M8)9», NRAST(3,2), I RECRY( 3 ,2 ) , . 

DECLARE 01 TEMPRECD, 02 CORRAY(320), 02 NRAY(320),. 

DECLARE FIRSTRASTERIJT (3,2), RASTERUTLAST ( 3 , 2 ) , . 

OFCLARF TDTFIL FILE RECORD DIRECT KEYED UPDATE ENVIRONMENT (F12560) 
MFDIUM ( SYSOOO, 2311) RFG ION AL ( 1 ) ) , . 

DECLARE INCOMPLETECT (3,2) , LARGM I SS I NGCT ( 3 , 2 ) , LARGQUESTCT (3,2) , 

I DECL AREDBADCT ( 3 ,2 ) , ITOTCT(3,2), I FRST I NDAY ( 3, 2) * . 

DFCLARE DTM AX ( 2 ) , I RASTRWANTED ( 2 ) , . 

ON ENDFILE (INTAPE) GO TO ENDBACKGROUND, . 

ON ENDFUE (SYSIN ) GO TO ENOSYSIN, . 

COLID * '01 05 10 15 20 25 30 

35 40',. 

CALL NTRAST , • 

UTO = 39490,. 

CALL REHF1L (INTAPE),. 

CALL REWFIL (OUTAPE),. 

OPEN FILE f TOTE IL ) , FILE (INTAPE), FILE (OUTAPE),. 

NRAST = 0,. 

INDICAV * f 0 * 8 , . 

ITOTCT = 0,. 

IFRSTINDAY = 0,. 

IENDN =0,. 

INCOMPLETECT * 0,. 

L ARGMI SS INGCT ^ 0,. 
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3 A 

3 7 

3R 

39 

40 

41 
4“» 
43 
^4 

45 

46 

47 
49 
51 

5? 

53 

56 

57 

58 

59 
69 
61 
4? 

63 

64 

65 

66 

67 

68 

69 

70 

71 
7? 

73 

74 


75 

7^ 

75 

70 

80 

51 

57 

83 

84 
86 
88 
80 
on 


L ARGQUFSTC T = Of 
I DEC L ARED3 ADCT = Of 
CORRAY - O.t. 

NRAY = Of 
IRECO = Of 

nn i = i to 3m 
on j - i to 
I R E C R Y { I, J| = IRFCD, . 
on jj = i to 6 f 
KEYUSEO = IRFCD,. 

WRITE FILE (TOTFIL) FROM (TEMPRECO) KEYFROM (KEYUSEO),. 

IRFCD = IRECD ♦ 1 END,. 

END,. END t# 

OFT T I LE C SYS I N ) FD IT ( UTBEGM, UTENO v DTMAX, l RASTP WANT ED, DUMMY ) 
(?PIl5 t 5), 4F(lo f 0), AUOMf 

TP OTMAX(l) LT 0. OR DTMAM2) LT 0. OR I R A SIR W ANTE D ( I) t E 0 
OR I P ASTRW ANTFDf ? ) LF 0 OR UTREGN GF UTENO THFN pn f . - 
PUT EDIT {'INPUT ERROR') (SKIP, A),. G0“ TO W°AP0P,. END”,. 
FLAPurnr = -i . f+60, . 

FLAOUTtn = 

UTIMF = — 100. F40 , . 

GO TO F T R S TB AC K GROUND, . 

NEXT8ACKGROUNO.. READ FILE I INTAPE) SET (RECIO),. 

FIRSTBACKGROUNO.. READ FILF ( INTAPE) SET (RECIO),. 

IF UTSinp LT UTBFGN THFN GO TO N FX TO ACKGROUNO 

IT UTSTOP - UTSTRT LT 0.3333 THEN GO TO NFX TO ACK GROUND , . 

IF UTSTRT GT UTENO * 0.333333 THEN GO TO 3 ACKGRHUNpT ( MEM I , • 

GO Tn TFSTRASTIME,. 

NEXTRASTER.. IENDNB4 = IFNON,. 

CALL NXRAST,. 

IF I EOF NF 0 THFN GO TO FNDF I LFR A S TFR S ♦ . 

TFSTRASTIME.. IF UTIMF LT UTSTRT THFN GO TO NF XTR AS TFR 
IF .UTIME LT , UTBFGN THEN GO TO NEXTRASTER,. 

IF UTIMF GT UTENO THEN GO TO M AXT I MR XCFFDEO ,. 

DO ! - 1 TO 3,. 

DO J = 1 TO 

IF NRAST (I, J) NE 0 THEN IF UTIME - F IR S TR AS TF R UT (I, J)' 

GT DTMAX (J) THFN CALL PR I NT R ASTFR (I, J) f . 

FND,. FNn,. 

IF UTIME GT UTSTOP THFN GO TO NEXT8ACKGR0UN0,. 

IF JFWEEL - 3 OR JFWEEL - 5 OR JAWEEL = 5 THFN GO T i) NEXTR^STf-R, 
IF 1ENDN = l THEN ISIZE = ISIZE - RO,. 

IF JFWFFL = 4 THFN JFWFEL = 3 t . 

IF JAWEFL = 3 THFN JAWEEL =? , . 

IF JAWFFL - 4 THFN JAWFFL = !,. 

TFSTFLARES.. IF UTIME LT FLARUTLO THFN GO TO NOFLARES , . 

IF UT l MR GT FLARUTHI THEN DO,. GET FILE (SYSIN) m(T (FLARUTLO 
FLARUTHI, DUMMY) (ZF(15,5), A(50)),.G0 TO TE ST F( AR F S t. FND , 
I DE C LARF OR ADCT { JFW EE L, JAWEFL ) “I DECLARE DBA OC I {JFWFEL , JAWFFL >♦) , 
GO TO NEXTRASTFR,. 

NOFLARES.. I TOTC T ( JFWFEL , JAWEFL) = ITOTCT I J FW FEL , J AWfFL ) ♦ l t . 
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91 

93 

95 

97 

9R 

99 

101 * 

103 

104 

105 * 

,106 

107 

.108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
12 ? 
128 

124 

125 

126 

127 

128 

129 ‘ 

130 

131 ‘ 

132 

133 1 

134 ' 

135 

136 

137 


IF 1SIZE LT 1920 THEN DO*. I NCOMPLET ECU JFWEEL * JAWEEL) - 
INCOMPLETECT I JFWEEL* JAWEEL) * 1 , . 

GO TO NEXTRASTER * • END* . 

IF IMISS GT 100 THEN DO*. LARGMl SS INGCT i JFWEEL* JAWEEL) » 
LARGMI SS INGCT (JFWEEL* JAWEEL) ♦ 1** GO TO NEXTRASTER*. 
END*. 

IF IBADCT GT 100 THEN DO, • LARGQUESTCT ( JFWEEL * JAWEE L ) = 

LARGQUESTCTC JFWEEL* JAWEEL ) + l * • GO TO NEXTRASTER*. FND, . 

/* THIS RASTER 1$ OK* SO ADD IT TO TOTALS */ 

IF I ENDNB4=l THEN IFRST INDAY ( JFWEEL , JAWEEL ) = IFRST INDAY 
'(JFWEEL* JAWEEL) + 1*. 

IF NRAST (JFWEEL* JAWEEL) * 0 THEN F IRSTRASTERUT ( JFWEEL* 

JAWEFL ) = UTIHE * . 

RASTERUTLASTf JFWEEL, JAWEEL) = UTIME*. 

IRECD = IRECRY (JFWEEL, JAWFEL),. 

NR A SJ 1 JFWEEL* JAWFEL) * NRAST (JFWEEL* JAWEEL) + 1*. 

I.NDEX1 - 0*. 

00 H = 1 TO 6». ' 

KEYUSED = IRECD*. 

READ FILE (T0TF1L) INTO (TEMPRFCD) KEY (KEYUSED)*. 

DO J = I TO 320*. 

INOEXl = INDEX! +1*. 

IF* IENDNB4 s 1 AND INDEXl LE 240 THEN GO TO SKIPWRD*. 

‘ IF IRASRY ( INDEX! ) AND 'OOOOOOIUB THEN GO TO SKIPWRD*. 

CALL PTWPDC (IRASRY* INDEXl'* N* NBAR ) , . 

CORRAY(J) = CORRAY(J) + NBAR*. 

NR AY( J ) = NRAY(J) + !,. 

SKIPWRD.. END*. 

REWRITE FILE (TOTFIL) FROM ( TEMPRECD) KEY (KEYUSED)*. 

IRECD a IRECD ♦ 1*. 

END* • ‘ 

IF NRAST ( JFWEEL * JAWEEL ) GE IRAS TRW ANTED! JAWEEL ) THEN CALL 
PRINTRASTER (JFWEEL, JAWEEL)*. 

GO TO NEXTRASTER*. 

ENDRACKGROUND.. PUT EDIT ("END OF FILE ON BACKGROUND RASTER T APE f ) 
(SKIPC2), A)*. * 

% ITOTCT * - 0*. 

'go to wrapup*. 

BACKGROUNDTTMEHI • • PUT EDIT (’MAXIMUM REQUESTED TIME EXCEEDED ON 8ACKGR 
OUNO RASTER TAPE’) (SKIPC2), A),. 

'GO TO WRAPUP,. 

ENDFILERASTERS.. PUT EDIT ( 1 END OF FILE ON TELEMETRY TAPE 1 ) (SKIP(2), 

A ) * . 

GO TO WRAP UP*. H 

MAXTIMEXCEEDED.. PUT EDIT (•MAXIMUM REQUESTED TIME EXCEEDED ON TELEMETR 
Y TAPE 1 ) (SKIP(2), A), . ' 

WRAPUP.. CALL RUNTL * • 

- ” CLOSE FILE (INTAPE)*. 

CALL RUNFIL (INTAPE),.- ' 

DO JFWEEL = 1 TO 3*. 

DO JAWEEL = V TO 2*. 


A-42 



nns nt/I COMPILER 360N-PI-464 CL1-5 


AVRCRAST 


05/11770 


/* asse nsn-o programs // get average corrected rasters */ 


1 13 
119 
141 
14^ 
143 
1 44 


145 

146 

147 
143 

149 

150 

151 
IS? 


153 
1 54 
1 55 
I 56 
157 
1 53 
1 59 
1 60 
161 
16? 
161 

164 

165 

166 

167 
163 
169 
1 71 
17? 
1 73 
1 74 
175 
t 76 
177 


CALL PRINTRASTER [ JF WELL » JAWEEL>»* 

' FNO * • END,* 

CLOSE PILE CT0TFIL1, FILE (OUTAPF)«. 

CALL RUNFIL ( HUT A PE ) t * 

DISPLAY ( # L A BEL OUTPUT TAPE ON SYSOOfl AS FOLLOWS — •>,. 

PUT STRING (MESSAGF) EDIT ( 1 11 AVERAGE CORRECTED RASTERS FROM* ,/ 
UTAVLOLn f « TO UTAVHIHIt (A, F(7 t l) f A, F(7 f l‘), A 

j ? • 

DISPLAY (MESSAGE),. 

STOP,.. 

FNOSYSTN.. PL ARUTLO = l.F*-60,. 

FLARUTHT - l.F+60,. 

r.n rn noflares,. 

ASIN.. PROCEDURE (X),. 

RFTIJRN (ATAN (X / SORT (1. - X * X)]),. 

END, . 

/* THIS PROCEDURE PRINTS OUT AN AVERAGE RASTER 
AND ASSOCIATED STATISTICS */ 

PRI NTR ASTFR . . PROCEDURE ( JFWEEL , JAWFEL) , . 

OFCLARF RIGR AY( 1D?0) CHARACTERS) 

DECLARE FNAME CHARACTER ( 24) , ANA«F CHAR ACTER { 20) , . 

DECLARE MONA Ml CHARACTER (-4) , M0NAM2 CHARACTER ( 4) , SGN CHAR AC TER ( 11 , . 
OECLARE KA Y ( 40 ) CHAR ACTFR ( 7 ) , PIC PtCTURE •<?•, IHRl, MINI, SECS1,. 
DECLARE I YRt , ( YRDATCl , YR0ATE2, DATEA) FLOAT BINARY ( 5R ) » FKNX , • EKE X , . 
DECLARE QLAML, OL AM? , [HR?, MIN2, SECS?, FN, FE, QLOG, QLOGTR , . 

OFCLARE QLDGFR, ITFMP, TF«P, M, J, INDEX1, I DATE I , IDATE? , t 0 ECO,IPIC,: 
OFCIARE I YR 7 , • ' . 

TF ITOTCT (JFWEEL, JAWEEL) LE 0 THEN RFTURN, . . • 

CALL SDAY ( FIRSTRASTFRUT( JFWEEL, JAWEEL > * UTO, IYRI, YPOAfF 1 , 
MONTH , DATFA, MONAMt ) , . 

I OAT E I’ = DATEA, . 

CALL DEGDMS ( 2 4. * ( OAT EA- 1 OA TE L ) , t SGN , IHRl, MINI, SKS1I,. 

CALL SDAY ( RA S TFRUTL AST ( J FWEFL , JAWFEL) ♦ UTO, IYR2, YR DATE? » 
MONTH, DATFA, M0NAM2 1 , . 

IDATE? = DATEA,. 

CALL DEGDMS { 24. * ( DA TEA- I DA T6? ) , I SGN, IHR?, «IN?, SFCS?),'. 

IF JFWEFL = T THEN M = 4,. ELSE M = JFWEEL,. 

CALL FLTLUK (M, FNAME, Ql AM 1 , QLAM2, FKNX, FKEX ) , . 

CALL APTLUK (JAWTEL, ANAME1,. 

FN = FKNX / FPS,. 

FE = FKEX / EPS,. 

CALL NEWPAGF , . 

PUT EDIT ( • • ) (SKIPC2) , A) ,. 

PUT FDIT UTOTCTt JFWEEL , JAWFEL), 'RASTERS RCAD' , 

INCOMPLETECT( JFWFEL, JAWEEL), 'WERF I NCOMPLFT.F ' , 
LARGMISSINGCT( JFWFFL, JAWFFL), 'HAD MORF THAN 100 INT FRSPFR 
SED MISSING WORDS', 

LARGOUESTCTI JFWEEL, JAWEEL), 'HAD MORE THAN 100 QUESTIONABL 

E WORDS' , 

IDECLAREDRADCTI JFWFFL, JAWEEL), 'RASTERS WERE IN TIME I NT F R 
VALS DFCLAPFD AS BAD', 
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17* 
179 
1 80 
1 33 
134 
185 
18* 
187 
133 

189 

190 

191 
19 ? 

193 

194 

195 
194 
197 
193 

199 

200 
201 
20 ? 
? 03 ~ 

207 

208 
209 
21 2 

213 

214 
715 
216 
217 
718 
27 ? 

223 

224 
227 
728 

229 - 

230 

231 . 

232 
73t3 
234 
735 


NR4ST( JFWEEL, JAWEEL), ♦.RASTERS WERE USED'* 

I'FR ST INDAY ( J FWEEL • JA WEED, 'RASTERS WERE 1ST SINCE SUNRISE 

♦ ) 

(SKIP, F ( 6*01,' X(l) , A1 , • 

* LOCATE OUTRECD FILE (OUTAPEI SET (RECIDO),. 

UTAVLO = FIRSTRASTERUTI JFWEEL, JAWEEL 1 ,. 

IF NOT IN01CAV THEN DO* . INDICAV^' l ♦ 8 * • UTAVLOLO = UTAVLO*. 

FND * • 

UTAVHI = RASTERUTLAST (JFWEEL, JAWEEL 1 , « 

UTAVHIHI =» UTAVHI 

I TOTCTAV = ITDTCT (JFWEEL, JAWEEL)*. 

JFWEEL AV - H f . 

JAWEELAV = JAWEEL*. 

rNCDMPLETECTAV = INCOMPLETECT ( JFWEEL , JAWEEL),. 

LARGM I SSI NGCTAV = LARGM I SS I NGCT (JFWEEL* JAWEEL)*. 

LARGQUESTCTAV = LARGOUESTCT (JFWEEL* JAWEEL)*. 
inECLAREDBADCTAV - 'I DECLAREDBADCT (JFWEEL, JAWEEL)*. 

NRASTAV = NRAST (JFWEEL, JAWEEL),. 

I FR ST INDAYAV * I FRST INDAY (JFWEEL, JAWEEL),. 

AVRAY « 0.,. * 

IRECD * IRFCRYt JFWEEL, JAWEEL),. 

INDEX1 - Ot. 

DO M a I TO 6,. 

KEYUSED * IRECD,. . 

READ FILE (TOTFIL) INTO (TEMPRECO) KEY (KEYUSF01*. 

DO J = l TO 320,. 

INDEX1 = INDEX1 + l*. 

rFNRAY(J) LE 0 THEN DO* . BIGRAY( INDEX! )*♦♦♦*,. GO TO ENDJ, • END, . 
AVRAY (INDEXU = CORRAY(J) / NRAY(J) - BAKR AY(t NDEX1 1 * . 

TEMP ^ AVRAY ( I N0EX1 ) + 1.0,. 

IF TEMP LT 1.0 THEN DO*. 81 GR AY( I NDEX1 ) = GO TO ENDJ*. 

END,.. . 

QLOG = LOG? ( TEMP) *. 

QLOGTR = TRUNC (QLOG) , . 

QLOGFR = (QLOG - QLOGTR) * 10.,. 

TTEMP = QLOGFR + 0.5,. 
iPIC, =: QLOGTR. + 1 . Q * • . 

IF ITEMP GT 9. THEN, DO,. ITEMP = 0,. IPIC = TP1C + I*. FND,. 

PIC = ITEMP*. 

SU8STR < BIGRAY ( INDEX1 ), 2, 1) - PIC*. 

IF IPIC GT 35 THEN DO*. BIGRAY ( I NOEX 1 ) « •*♦«,* GO TO FND J * . 
END*. 

SUB STR ( BIGRAY ( I NDEXl ) , 1,1) ^ SU8STR ( f 0 1 23456789ABCDEFGH I JKLMN 
OPQRSTUVWXYZ', IPIC, ]>,* 

- ENDJ..- END* • 

CORRAY * 0.,. 

NRAY - 0*. 

REWRITE FILE (TOTFIL) FROM (TEMPRECO) KEY (KEYUSED),. 

IRECD = IRECD + I*.‘ 

END,. 

INCOMPLETECT (JFWEEL, JAWEEL) = 0*. 
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?36 

737 

738 
7 39 
240 
?41 
74? 
24“* 

H4 

745 

246 

747 

?48 
74Q 
750 
? 5 1 
75? 
253 
7 54 
255 

7 56 
757 
258 
759 
260 


761 


7 67 


263 

764 

765 
?66 


LARGMISSINGCTt JFWEEJ. , JAWEEU = Of 
ITnTCT (JFWEEL, JAWEEL) = 0,. 

LARGQUFSTCTC JFWEFL, JAWEEL) = n t . 

NR A ST (JFWFELt JAWEEL) = 0,* 

IFRST1NDAY ( JFWEFL , JAWEEL) = 0,. 

I DECLARED RADCT (JFWEEL, JAWFEL ) = 0,. 

CALL NFWPAGFf 

PUT EDIT {'LOG BASF 2 UP (AVERAGE CORRECTED COUNT RATE * 1)0 
COLID) { SK I P ( 2 ) v X ( 70 ) t A, SKIP, X(4), ACBO)),. 

DO M = 1 TO 48 1 * 

ILINF = 49 - M f . 

ILINFOUT = * 

IF l L I NE" 1 OR 1 L I NF=48 OR ( I L INE/5 ) *6-1 L I NF THEN ILlNfEPJC 
= ILINF,. 

KAY - » ',. 

00 J = 1 TO 39 BY ?,. 

I SUBSCR = 48 * (J-l) + 49 - M,. 

KAY(J) = B I GR AY ( I SUB SCR ) , . 

I SUBSCR = 48 * J + M t . 

KAY ( J+l ) = BIGRAY( ISUBSCP ) ♦. 

END,* 

PUT EDIT (ILINEOtJT, KAY, ILINEOUT) (SKtP, A{2), X(2), 40A(?>, 

X ( 7 ) , AC 2 ) ) , . 

END, . 

PUT EDIT { COL I 0 ) (SKIP, X C 4 ) , A(80)>,. 

RETURN,. 

NFWP AGE. . PROCEDURE,. 

PUT EDIT ( 1 AVERAGE CORRECTED RASTER FM ', IOATEl, MBNAM l , ivpi, 
I HR I , MINI, SECS1, 1 TO «, IDATE2, M0N4M? , IYR7, I HR? , 

MIN2, SECS2, U', FIRSTRASTERUTI JFWEFL, JAWEEL), 1 TO’, 
RASTFRUTLASTC JFWEEL, JAWEFL) , • ) f ) 

(PAGE, A, F ( 3 , 0 ) , X(l), A, ?F(5,Q), F(3,D), F(5,l), A, 

F ( 3 f 0 ) , X(l), A, 7 P { 5 f 0 ) , F C 3,0 ) ? F(5,l)* X(6), ?(A t F(9,7) 

) , A ) , . 

PUT EDIT (FNAME, 1 L AM BOA l QLAMl, • LAMBDA 0L4M7, 

■KN = ', FN, • KE FE) (SKIP, A, X(6), A, F(6,7), X(4>, 

A, F ( 6 , 7 ) , X ( 5 ) , A, F ( 8 , 2 ) , XI5), A, F ( 3 , 7 ) } , . 

PUT EDIT ( SUBS TR (ANAMC, 4, 17), •FPSILON EPS, 

•N BAR CAL =■♦ XNBARCALt •DELTA DEL T ) 

(SKIP, X ( 3 ) , A, XC10), A, F ( 7 , 7 ) , X(IO), A, F(3,2), XC9),A, 
F { 8 , 7 ) ) 

RFTURN,. 

END,. 

END,. 

^ FNO , . 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
;22 

23 

24 

25 

26 
27 
'28 

29 

30 

31 

32 

33 

34 


35 

36 

37 

. *. > -‘38 1 < 

39 

40 

41 

42 

43 

44 

45 

46 

47 


/* AS LE OSO-D PROGRAMS // ASPECT— EPHEMER IS SUMMARY LISTING */ 

OSSM: PROCEDURE OPTIONS (MAIN); 

DECLARE (REWFILt RUNFIL, FLCN93, OEGDMS) ENTRY; 

DECLARE UTIME FLOAT BINARY (53 ) \ 

DECLARE TAPEIN FILE INPUT RECORD BUFFERED SEQUENTIAL ENVIRONMENT 

(U ( 3 1 20 ) LEAVE MED I UMI SYS007 « 2400) BUFFERSt 2) NOLABEL) ; 
OECLARE 01 RECORDIN ALIGNED, 02 WGRD_7090 (520) BITI48); 

ON RECORD ( TAPEIN )T 

CALL REHFIL (TAPEIN) ; 
l PAGE = o; 

LINES. = 80; 

NEWSFILE: OPEN FILE (TAPEIN); 

ON ENDFILE(TAPEIN) GO TO ENO_RUN; 

DO l = 1 TO 4; 

- READ FILE (TAPEIN) INTO (RECORDIN); 

end; 

READ FILE (TAPEIN) INTO (RECORDIN); 

IORBIT - FLCN93 ( WORD_7090 ( 54 )) ; 

SECS = FLCN93 ( W0RD_709G{ 56J ) / 1QG0.; 

DAYNO = FLCN93 ( WORD_7090 ( 55 ) ) ; 

UTIME = DAYNO + <*490, ; 

UTIME UTIME + SECS / 86400,; 

UT_D AY-START = FLCN93 ( W0RD_7090( 63 1) ; 

UT_DAYEXIT = FLCN93 { WORD_7090 (6 1 ) ) ; 

PITCH = FLCN93 { WORD_7090 (28)) 5 
ROLL = FLCN93 ( WOR0_7090 (29)); 

SPINX = FLCN93 ( W0RD_7090< 42) ) ; 

TEMP = FLCN93(W0RD_7090(43) ); 

TEMPI = FLCN93(W0RD_7090(44) ) ; 

' SPINY = TEMP * 0.91-74532 + TEMPI * 0.3978437; 

SPINZ = - TEMP * 0.3978437 + TEMPI * 0.9174532; 

SUNX FLCN93 ( WORD__7090( 16) ) ; 

TEMP = .FLCN93(WOR0l7090( 17)) ; 

TEMPI = FLCN93 ( W0RD_7090 ( 181); 

SUNY = TEMP * 0.9174532 + TEMPI * 0.3978437; 

SUNZ = - TEMP * 0.3978437 + TEMPI * 0.9174532; 

*/* ASCE 0S0-4 SUBROUTINES // COMPUTE ANGLE BETWEEN SPIN AXIS AND THE 
. SUN* S NORTH POLE */ 

S PD = SORT ( S P I NX*S P I NX * SPINY*SPINY + SP INZ*SP 1NZ ) , . 

SUND = SQRT (SUNX*$UNX + SUNY*SUNY + SUNZ*SUNZ)t. 

SPINX = SPINX / SPD, • 

SPINY = SPINY / SPD,. 

SPINZ = SPINZ / SPD,. 

SUNX = SUNX / SUNOf. 

SUNY = SUNY / SUND,. 

SUNZ = SUNZ / SUND, . 

FO = ASIN (SPINX * SUNX ♦ SPINY * SUNY ♦.SPINZ * SUNZ),. 

XL = ATAN (SUNY, SUNX),. 

R =* 1.307899 + 0.6146465E-06 * (UTIME); 

BO = ASIN (SIN(XL-R) * 0.126199),. 

ALPHA = ACOS (0.992005 / C0S(80>); 
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48 

49 

50 
5 l 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
6? 

63 

64 

65 

66 

67 

68 

69 

70 

71 
7? 

73 

74 

75 

76 


77 

78 

79 
30 
81 
82 

83 

84 

85 

86 

87 

88 
89 


90 

9t 


R_l = R + 0 . 614665 E- 06 ; 

= AS I N CSINIXL - R_ 1 ) * 0 . 126199 ); 

U- IF 8 __ l - 80 <= 0 THEN ALPHA = - ALPHA; 

RX = SIN(R) * 0 * 126199 ,. * 

RY = - COS ( R ) * 0 . 126199 ,. 

COST 2 - SPINX * RX + SPINY * RY + 0.992005 * SPINZ,. 

V = ACOS ( C 0 ST 2 ) f * 

S - 0*5 * ( 3.141593 + FO + V - 3 n>,» 

RP - ((SIN ( 0 . 5 * 1 V— FO— BO ) ) * SINCO . 5 *( 3 . 141693 +FO-V-BO) ) 

* SIN (O. 5 *(F 0 +V+B 0 ) ) ) / SIN(S)),. 

RP = SORT (ABS (RP) ) t . 

AR G 4 ' ” SIN (S-Vi 

PR 001 = SPINZ * SUNY' — SPINY * SUNZ,. 

PR 00 2 = SPINX * SUNZ - SPINZ *'SUNX,. 

PR OO 3 = SPINY * SUNX - SPINX * SUNY,. 

STEST - PROD 1 * RX * PRQD 2 ^ PY f PRHD 3 * 0 . 992004 ,. 

ARG 3 = RP * SIGN (STEST * ARG 4 ) • • 

DELTA = - 2.0 * ATAN ( ARG 3 , ARG 4 ) + 1 . 570796 ; 

HRS - SECS / 3600 .; 

CALL DFGOMS (HRS, ISGN, I HR , IMIN, SFC); 

HRS = UT^DAYST ART / 3600000 .; 

CALL 0 EGOMS (HRS, ISGN, I HR 1 , I 31 N 1 , SEC 1 ); 

HRS = UT_DAYEXIT / 3600000 .; 

CALL OEGOMS (HRS, ISGN, 1 HR 2 , IMIN 2 , SEC 2 ) ; 

ORB I T_END = FLCN 93 ( wnPD_ 7090 ( 33 ) > / 3600000 .; . 

CALL OEGOMS '(CRBIT^FND, ISGN, IHR 3 , IMIN 3 , SFC 3 ) ; 

C ALPHA = - ALPHA; 

^PITCH * - PITCH; 

CALL NEW_PAGE_TFST ; 

PUT EDIT (IORBIT, DAYNU, 1 HR, IMIN, SFC* I HR 3 , I*IN 3 , SEC 3 , 
IHR 1 , IMIN 1 , SFC 1 , IHR 2 , IMIN 2 , SFC 2 f P I TCH* 5 ( . 29 577 ^ 5 , 
ROLL* 57 . 2957796 , ALPHA* 47 . 2957796 , DC L TA*S 7 . 295 7 795 ) 
(SKIP( 1 ), 2 F ( 5 ) , X( 5 ), 3 F ( ^ | , X( 3 ) T 3 F( 3 ), X{ 7 ), 3 T { 3 1 , 
X( 2 ) , 3 F( 3 ) , 4 F { 10 , 2 M ; 

ON END P l LE( T APEI N ) GO TO END_OF_F ILE ; 

NEXT_RECORO: READ FILE (TAPEIN) INTO (RECOROIN); 

GO TO NEXT_REC 0 RO; 

FND_OF_F ILE; CLOSE FILE (TAPEIN); 

GO TO NEW_F ILF? 

END^RUN: CALL RLJNFIL (TAPEIN); 

stop; 

NFW_PAGE_TEST: PROCEOURE ; 

LINES = LINES + 1 ? 

IF LINES < 50 THEN RETURN ; 

LINFS = o; 

IPAGE = IPAGE + 1 ? 

PUT EDIT ( f OSn— D SUMMARY LISTING FROM ATTITUDE-EPHEMERIS* , 
f ORB I r T DAY ORBIT START ORBIT END SUN ENTRANCE SUN EXIT' 

TCH ROLL ALPHA DELTA 1 ) (PAGE, A ? SKiP(l), A); 

RETURN; 

END; 
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92 

ASIN.. 

PROCEDURE ( XI , . 


93 


RETURN (ATAN (X / 

SORT (1. - X 

94 

END,. 



95 

ACOS • . 

PROCEDURE IX),. 


96 


ANGLE = ATAN (SORT 

I 1.0 - X*X) , 

97 

* 

RETURN (ANGLE),. 


98 

END,. 



99 

END, . 




*/ 

* XI)),. 
XI ; 



